作者Anchorites (隱客士)
看板R_Language
標題[問題] 非線性回歸, nlme
時間Sun Mar 17 18:02:48 2019
[問題類型]:
給定方程式形式，
非線性回歸。
[軟體熟悉度]:
今天第一次用...
但是時間非常緊迫，
查了一個下午沒結果，
被逼得只好上來問。
[問題敘述]:
給定方程式形式，要決定係數。
方程式形式如下：
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
-----
Sent from JPTT on my OPPO CPH1725.
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 1.200.222.25
※ 文章網址: https://webptt.com/m.aspx?n=bbs/R_Language/M.1552816970.A.72E.html
1F：→ Anchorites: 補充，Dsr75也是已知數據 03/17 18:04
2F：→ andrew43: nlme主要是用來做mixed model的 03/17 21:27
先謝謝回覆
不過我有寄給paper原著，他說他當初是用nlme package裡面的nlme function跑的沒錯。但是我去翻instruction完全看不懂input的內容 QQ
3F：推 bboybighead2: 看起來用lm跑就好了，只是你要先對變數做處理再丟進 03/20 09:56
4F：→ bboybighead2: 去。Dsr取自然對數log(Dsr)、Mw不處理、log(sqrt( 03/20 09:56
5F：→ bboybighead2: Rrup^2+hi^2))、log(sqrt(Rrup^2+hi^2))*Mw、log( 03/20 09:56
6F：→ bboybighead2: Vs30)、Ztor不處理。這些變數全丟進lm就可以得到你 03/20 09:56
謝謝回覆
可是lm不是用來跑linear model的嗎？我的方程式是nonlinear了，這樣做是可行的嗎？(我自己有用lme跑過，被系統擋下來，上網查原因似乎是因為起始值給的不好)
我的參數其實都已經匯入R準備好了，現在就是只差在如何"使用"nlme這個function (paper原著也有提到另一個叫lme4的新方法也適用)
7F：→ bboybighead2: 想要的迴歸係數了 03/20 09:56
※ 編輯: Anchorites (180.217.137.188), 03/20/2019 15:55:03
※ 編輯: Anchorites (180.217.137.188), 03/20/2019 15:58:21