[問題] 非線性回歸, nlme

作者: Anchorites (隱客士)   2019-03-17 18:02:48
[問題類型]:
給定方程式形式,
非線性回歸。
[軟體熟悉度]:
今天第一次用...
但是時間非常緊迫,
查了一個下午沒結果,
被逼得只好上來問。
[問題敘述]:
給定方程式形式,要決定係數。
方程式形式如下:
http://i.imgur.com/MiI8cjm.jpg
其中,Mw, Rrup, Vs30, Ztor這4個是已知的參數,要用600筆有以上4個參數的數據,回歸出c0, m1, r1...等7個係數。
結果大概會像這樣
http://i.imgur.com/42SXTPf.jpg
包含回歸出的係數值及其95%信賴區間。
目前已經知道nlme package似乎可以解決我的問題,也有查過youtube影片。
但nlme好像能處理的非線性,
只有類似像,
y = a*x1*x2 + b*x1^2 + c
(每個未定係數都是某項的領導係數)
然後R可以決定a, b, c。
但是我的function form太醜了,
在網上找不到類似的範例可以依樣畫葫蘆。
跪求大神解惑
萬分感謝
[程式範例]:
(略)
[環境敘述]:
macOS 10.14.3
[關鍵字]:
特定方程式形式
非線性回歸
nlme
作者: Anchorites (隱客士)   2019-03-17 18:04:00
補充,Dsr75也是已知數據
作者: andrew43 (討厭有好心推文後刪文者)   2019-03-17 21:27:00
nlme主要是用來做mixed model的
作者: bboybighead2 (腳毛會唱歌)   2019-03-20 09:56:00
看起來用lm跑就好了,只是你要先對變數做處理再丟進去。Dsr取自然對數log(Dsr)、Mw不處理、log(sqrt(Rrup^2+hi^2))、log(sqrt(Rrup^2+hi^2))*Mw、log(Vs30)、Ztor不處理。這些變數全丟進lm就可以得到你想要的迴歸係數了

Links booklink

Contact Us: admin [ a t ] ucptt.com