R中多元非线性回归

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

是对y1=a*e^b*x+c*sin(w*t+u)回归，要求a，b，c，w，u； R程序我是这样编的
m1=nls(y1~a*exp(b*x)+c*sin(w*t+u),start=list(a=1,b=0.1,c=100,w=0.01,u=0.1),trace=T)
但是一直提示奇异梯度，或产生缺省值，请问是不是list(a=1,b=0.1,c=100,w=0.01,u=0.1)的问题的，
要怎样改呢？谢谢
x=c(1:12)y1=c(7098.00,7498.00,7848.00,8254.00,8761.00,8801.12,8951.32,9325.03,9683.90,10200.00,11000.00,12360.74)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：多元非线性回归 非线性回归 线性回归 非线性 Start start 程序

t是什么？

qoiqpwqr 发表于 2011-12-21 10:59
t是什么？

m1=nls(y1~a*exp(b*x)+c*sin(w*x+u),start=list(a=1,b=0.1,c=100,w=0.01,u=0.1),trace=T)
，额，是x。崩溃掉，这样改后还不行

试一下(a=1,b=0.8,c=100,w=0.01,u=0.1)，就可以了
初始值的选取很重要

qoiqpwqr 发表于 2011-12-21 11:50
试一下(a=1,b=0.8,c=100,w=0.01,u=0.1)，就可以了
初始值的选取很重要

谢谢，那请问 m1=nls(y1~a*x+b+K/(1+exp(c-v*x)),start=list( )这里边初始值应该选什么，选初值有什么技巧吗？   这里是对y1=a*x+b+K/(1+e^(c-v*x))   的回归，求a,b,K ,c,v  

iamwilling 发表于 2011-12-21 12:21
谢谢，那请问 m1=nls(y1~a*x+b+K/(1+exp(c-v*x)),start=list( )这里边初始值应该选什么，选初值有什么技巧 ...

如果你的模型里面的参数有些实际意义的话，可以根据实际情况来选初值。不然的话，就大概试几个吧。

qoiqpwqr 发表于 2011-12-21 12:31
如果你的模型里面的参数有些实际意义的话，可以根据实际情况来选初值。不然的话，就大概试几个吧。

  m1=nls(y1~a*x+b+K/(1+exp(c-v*x)),start=list(a=1,b=100,K=100,c=100,v=1),trace=T)
前面的回归语言有错码，我毫无头绪地试了几个，无果。。。

先编一个最优函数，然后在用optim（）函数也可以得到参数估计

qoiqpwqr 发表于 2011-12-21 12:31
如果你的模型里面的参数有些实际意义的话，可以根据实际情况来选初值。不然的话，就大概试几个吧。

m1=nls(y1~a*x+b+K/(1+exp(c-v*x)),start=list(a=1,b=100,K=100,c=100,v=1),trace=T)
这个初值我试了很多个还不行，可以帮忙试一下吗？非常感谢你~~~