所以我准备跳楼了,因为我无法弄清楚这个代码是否适合我的生活。每次我运行它,我得到的标准误差为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)
什么是错的。我感谢任何和所有的帮助,我是自我引导的新手,我很擅长编码。