R语言代码请教

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我手上有两组数据，想用y=ax^b来拟合，在使用R运行的时候，一直报错。求大佬们指点一下，怎么改代码才能算出估计值a 、b

x<-c(2000,2001,2002,2003,2004,2005,2006,2007,2008,2009,2010,2011,2012,2013,2014,2015,2016,2017,2018,2019)
y<-c(8.6,8.9,9.5,11,15.2,20,26.2,33.27,43.41,62.06,101.72,131.15,170.73,127.67,161.55,271.36,287.47,301.04,378.97,411.75)
cor(x,y)
test=nls(y~a*x^b,start=(list(a=1,b=1)))

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

您这本来是时间序列，x是年份，y是观测值。
2013，2014两年的数据偏离明显，
改成
test=nls(y~a*(x-1999)^b,start=(list(a=1,b=1)))
再试试。

建议你x,y取对数后线性回归

llb_321 发表于 2020-4-7 16:10
您这本来是时间序列，x是年份，y是观测值。
2013，2014两年的数据偏离明显，
改成

问题解决了 谢谢大佬