[问题] 非线性回归, 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