【R语言新手】数据挖掘非线性回归

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

fdat=read.table(file="lll.txt",header = T)
YM=fdat$Y/fdat$M
LM=fdat$L/fdat$M  
KM=fdat$K/fdat$M  
summary(YM)
summary(LM)
summary(KM)
nl.f=nls(YM~A*(LM^a)*(KM^(1-a)),start=list(A=0.5,a=1),trace=T)
summary(nl.f)

library("maxLik")
loglik=function(para){
  N=length(YM)
  e=YM-para[1]*LM^para[2]*KM^(1-para[2])
  ll=-0.5*N*(1+log(2*pi)-log(N))-0.5*N*log(sum(e^2))
  return(ll)
}
res=maxLik(loglik,start=c(0.5,1),method="NR")
summary(res)

plot(YM~LM)
lines(LM, fitted(nl.f))
plot(YM~KM)
lines(KM, fitted(nl.f))
Rsq<-1-sum(resid(nl.f)^2)/(sum((YM-mean(YM))^2))
Rsq

adjRsq<-1-(length(YM)-1)/(length(YM)-2)*(1-Rsq)
adjRsq

anova(nl.f,lm.f)    ！！！object 'lm.f' not found

Q = -2*(logLik(nl.f)-logLik(lm.f))
df.Q=df.residual (nl.f)-df.residual(lm.f)
1-pchisq(Q,df.Q)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：非线性回归 数据挖掘 线性回归 R语言 非线性

倒数第四行   lm.f  not found
什么意思呀  怎么改进

帮楼主顶贴~~·