关于《R语言实战》中交叉验证内容的疑问

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

R语言实战一书使用bootstrap包中的crossval()进行交叉验证
代码为
shrinkage<-function(fit,k=10){
  require(bootstrap)

  theta.fit<-function(x,y){lsfit(x,y)}
  theta.predict<-function(fit.x){cbind(1,x)%*%fit$coef}

  x<-fit$model[,2:ncol(fit$model)]
  y<-fit$model[,1]

  results<-crossval(x,y,theta.fit,theta.predict,ngroup = k)
  r2<-cor(y,fit$fitted.values)^2
  r2cv<-cor(y,results$cv.fit)^2
  cat("Original R-square=",r2,"\n")
  cat(k,"Fold Cross-Validated R-square=",r2cv,"\n")
  cat("Change=",r2-r2cv,"\n")
}



states<-as.data.frame(state.x77[,c("Murder","Population",
                                   "Illiteracy","Income","Frost")])
fit<-lm(Murder~Population+Illiteracy+Income+Frost,data=states)
shrinkage(fit)


我对照数次都觉得我并没有输入错误，但是R却显示如下信息
Error in theta.predict(u, x[groups[[j]], ]) :
  unused argument (x[groups[[j]], ])

求解答TAT谢谢大大！！

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：R语言实战 交叉验证 R语言 shrinkage Bootstrap

我也是新手。。不明白为什么你的不对，我刚才对着你的一遍，没有找出来，但是我的可以运行，给你看看我的代码吧。


shrinkage<-function(fit,k=10){
require(bootstrap)
theta.fit<-function(x,y){lsfit(x,y)}
theta.predict<-function(fit,x){cbind(1,x)%*%fit$coef}
x<-fit$model[,2:ncol(fit$model)]
y<-fit$model[,1]
results<-crossval(x,y,theta.fit,theta.predict,ngroup=k)
r2<-cor(y,fit$fitted.values)^2
r2cv<-cor(y,results$cv.fit)^2
cat("Original R-square=", r2,"\n")
cat(k,"Fold Cross-Validated R-square =",r2cv,"\n")
cat("Change =",r2-r2cv,"\n")
}
fit<-lm(Murder ~ Population + Income + Illiteracy +Frost,
data=states)
shrinkage(fit)

少女AAAAA 发表于 2017-3-9 15:57
我也是新手。。不明白为什么你的不对，我刚才对着你的一遍，没有找出来，但是我的可以运行，给你看看我的代 ...

谢谢谢谢！！