[問題] 非線性回歸, nlme
[問題類型]:
給定方程式形式，
非線性回歸。
[軟體熟悉度]:
今天第一次用...
但是時間非常緊迫，
查了一個下午沒結果，
被逼得只好上來問。
[問題敘述]:
給定方程式形式，要決定係數。
方程式形式如下：
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
→
03/17 18:04,
6年前
, 1F
03/17 18:04, 1F
→
03/17 21:27,
6年前
, 2F
03/17 21:27, 2F
先謝謝回覆
不過我有寄給paper原著，他說他當初是用nlme package裡面的nlme function跑的沒錯。但是我去翻instruction完全看不懂input的內容 QQ
推
03/20 09:56,
6年前
, 3F
03/20 09:56, 3F
→
03/20 09:56,
6年前
, 4F
03/20 09:56, 4F
→
03/20 09:56,
6年前
, 5F
03/20 09:56, 5F
→
03/20 09:56,
6年前
, 6F
03/20 09:56, 6F
謝謝回覆
可是lm不是用來跑linear model的嗎？我的方程式是nonlinear了，這樣做是可行的嗎？(我自己有用lme跑過，被系統擋下來，上網查原因似乎是因為起始值給的不好)
我的參數其實都已經匯入R準備好了，現在就是只差在如何"使用"nlme這個function (paper原著也有提到另一個叫lme4的新方法也適用)
→
03/20 09:56,
6年前
, 7F
03/20 09:56, 7F
※ 編輯: Anchorites (180.217.137.188), 03/20/2019 15:55:03
※ 編輯: Anchorites (180.217.137.188), 03/20/2019 15:58:21
R_Language 近期熱門文章
PTT數位生活區 即時熱門文章
