全部版块 我的主页
› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛
R中多元非线性回归
iamwilling
11772 10
2011-12-21
是对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)

二维码

扫码加我 拉你入群

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

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

全部回复
qoiqpwqr
2011-12-21 10:59:53
t是什么？
二维码

扫码加我 拉你入群

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

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

iamwilling
2011-12-21 11:06:14
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。崩溃掉，这样改后还不行
二维码

扫码加我 拉你入群

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

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

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

扫码加我 拉你入群

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

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

iamwilling
2011-12-21 12:21:59
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  
二维码

扫码加我 拉你入群

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

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

qoiqpwqr
2011-12-21 12:31:19
iamwilling 发表于 2011-12-21 12:21
谢谢，那请问 m1=nls(y1~a*x+b+K/(1+exp(c-v*x)),start=list( )这里边初始值应该选什么，选初值有什么技巧 ...
如果你的模型里面的参数有些实际意义的话，可以根据实际情况来选初值。不然的话，就大概试几个吧。
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

微信QQ空间QQ微博
扫码加好友，拉您进群
各岗位、行业、专业交流群