我现在搞了一个新的某肿瘤的生存预测模型。生存分析用的C-index做出来了。然后需要继续做Bootstrap。但是无论如何做Bootstrap都提示报错。
Bootstrap用的是以下R语言代码：bootit=200
for(i in 1:bootit)
{
case=noNA[group=="long",]
control=noNA[group=="<24",]
bootindex.case=sample(1:nrow(case),replace=T)
boot.case.data=case[bootindex.case,]
bootindex.control=sample(1:nrow(control),replace=T)
boot.control.data=control[bootindex.control,]
boot.data=rbind(boot.case.data,boot.control.data)
dstr.boot=svydesign(id=~1, prob=~inv_weight, fpc=~ssize, data=boot.data)
boot.fit=svycoxph(Surv(mysurC$Time,mysurC$Statuso) ~AAT1, data=boot.data, x=TRUE,design=dstr.boot)
library(Hmisc)
library(survival)
cindex.train=1-rcorr.cens(lp.boot,Surv(boot.data$time, boot.data$Statuso))
cindex.test=1-rcorr.cens(lp_=.test,Surv(boot.data$time,boot.data$Statuso))
bias=rep(1,bootit)
bias=abs(cindex.train-cindex.test)
}报错原因总是：Error in noNA : object 'noNA' not found。请求会的帮忙解答下。
