关于将计算cross-validation写入for循环体的疑问

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我想做的是用一个循环体遍历所有的2变量模型，并计算CV值。用一个list vector（logi）储存glm得出的模型，再将logi[[i]]作为参数传递到cv.glm中，但是这样操作会报错。
下面是我的代码
library(boot)
des=read.table('XXXXX')
k=1
logi=vector('list',15)
cv=rep(0,15)
cost <- function(r, pi = 0) mean(abs(r-pi) > 0.5)
for(i in 5:9){
  for(j in (i+1):10){
    logi[[k]]=glm(V1~des[,i]+des[,j],data = des,family = binomial)
    cv[k]=cv.glm(des,logi[[k]],cost=cost,K=10)$delta[1]
    k=k+1
  }
}
错误信息见图，有大神解答一下吗，谢谢，想了好久

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Validation ATION Valid for循环 Cross library family 模型 信息

错误处换成用数据集变量名试试

nuomin 发表于 2016-9-13 09:47
错误处换成用数据集变量名试试

错误处是：
cv[k]=cv.glm(des,logi[[k]],cost=cost,K=10)$delta[1]。
这个函数会用到上面那个模型V1~des[,i]+des[,j],data = des,family = binomial
我试过把V1~des[,i]+des[,j]换成V1~V2+V3确实不会报错，但是就不能用循环体了，只能一个一个算

feizheyao 发表于 2016-9-13 10:03
错误处是：
cv[k]=cv.glm(des,logi[[k]],cost=cost,K=10)$delta[1]。
这个函数会用到上面那个模型V1~de ...

用names（）【i】函数把变量名提出来

      就本程序来说：大双环中：
     xvars<-names(des)
     xvars<-xvars[!xvars%in%c("V1")]
  for...
      for....
        fm<-paste(xvars[c(i,j)],collapse="+")
        fm<-as.formula("V1~",fm,sep="")
        logi[[k]]<-glm(fm,data=des,family=binomial)

fm<-as.formula("V1~",fm,sep=""),改错为
  fm<-as.formula(paste("V1~",fm,sep=""))
#

nuomin 发表于 2016-9-14 08:59
用names（）【i】函数把变量名提出来

程序已改对了，谢谢！

jgchen1966 发表于 2016-9-14 15:06
fm

程序已改对了，谢谢！