R中的bootstrapping代码不起作用?

时间:2016-04-19 18:47:27

标签: r bootstrapping

所以我准备跳楼了,因为我无法弄清楚这个代码是否适合我的生活。每次我运行它,我得到的标准误差为0(偏差= 0)所以显然我做错了,因为当我想要1000次迭代时它使用相同的样本。

我想要所有这些样本的模型的AIC。我创建了一个数据框并将其命名为“d。”

这是我的代码:

bootaic<-function(d, idx){
d2<-d[idx,]
fit<-coxph(Surv(time,infection) ~ X48FL + X55FL, method="efron", data=d)
return(AIC(fit))
}

bootf<-boot(data=d,statistic=bootaic,R=1000)

什么是错的。我感谢任何和所有的帮助,我是自我引导的新手,我很擅长编码。

0 个答案:

没有答案