[問題] 非線性回歸, nlme

[Anchorites (隱客士)]

[問題類型]:
給定方程式形式，
非線性回歸。


[軟體熟悉度]:
今天第一次用...
但是時間非常緊迫，
查了一個下午沒結果，
被逼得只好上來問。

[問題敘述]:
給定方程式形式，要決定係數。
方程式形式如下：
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://www.ptt.cc/bbs/R_Language/M.1552816970.A.72E.html
→ Anchorites: 補充，Dsr75也是已知數據
→ andrew43: nlme主要是用來做mixed model的
先謝謝回覆
不過我有寄給paper原著，他說他當初是用nlme package裡面的nlme function跑的沒錯。但是我去翻instruction完全看不懂input的內容 QQ
推 bboybighead2: 看起來用lm跑就好了，只是你要先對變數做處理再丟進
→ bboybighead2: 去。Dsr取自然對數log(Dsr)、Mw不處理、log(sqrt(
→ bboybighead2: Rrup^2+hi^2))、log(sqrt(Rrup^2+hi^2))*Mw、log(
→ bboybighead2: Vs30)、Ztor不處理。這些變數全丟進lm就可以得到你
謝謝回覆
可是lm不是用來跑linear model的嗎？我的方程式是nonlinear了，這樣做是可行的嗎？(我自己有用lme跑過，被系統擋下來，上網查原因似乎是因為起始值給的不好)
我的參數其實都已經匯入R準備好了，現在就是只差在如何"使用"nlme這個function (paper原著也有提到另一個叫lme4的新方法也適用)
→ bboybighead2: 想要的迴歸係數了
※ 編輯: Anchorites (180.217.137.188), 03/20/2019 15:55:03
※ 編輯: Anchorites (180.217.137.188), 03/20/2019 15:58:21