cox单因素独立预后分析提示错误

cox单因素独立预后分析提示：
Error in coxph(Surv(futime, fustat) ~ rt[, i], data = rt) :
No (non-missing) observations

代码如下：
定义独立预后分析函数
indep=function(riskFile=null, cliFile=null, project=null){
risk=read.table(riskFile, header=T, sep="\t", check.names=F, row.names=1) #读取风险文件
cli=read.table(cliFile, header=T, sep="\t", check.names=F, row.names=1) #读取临床文件

#数据合并
sameSample=intersect(row.names(cli),row.names(risk))
risk=risk[sameSample,]
cli=cli[sameSample,]
rt=cbind(futime=risk[,1], fustat=risk[,2], cli, riskScore=risk[,(ncol(risk)-1)])

#单因素独立预后分析
uniCoxFile=paste0(project,".uniCox.txt")
uniCoxPdf=paste0(project,".uniCox.pdf")
uniTab=data.frame()
for(i in colnames(rt[,3:ncol(rt)])){
     cox <- coxph(Surv(futime, fustat) ~ rt[,i], data = rt)
     coxSummary = summary(cox)
     uniTab=rbind(uniTab,
                  cbind(id=i,
                  HR=coxSummary$conf.int[,"exp(coef)"],
                  HR.95L=coxSummary$conf.int[,"lower .95"],
                  HR.95H=coxSummary$conf.int[,"upper .95"],
                  pvalue=coxSummary$coefficients[,"Pr(>|z|)"])
                  )
}
write.table(uniTab,file=uniCoxFile,sep="\t",row.names=F,quote=F)
bioForest(coxFile=uniCoxFile, forestFile=uniCoxPdf, forestCol="green")


#多因素独立预后分析
multiCoxFile=paste0(project,".multiCox.txt")
multiCoxPdf=paste0(project,".multiCox.pdf")
uniTab=uniTab[as.numeric(uniTab[,"pvalue"])<1,]
rt1=rt[,c("futime","fustat",as.vector(uniTab[,"id"]))]
multiCox=coxph(Surv(futime, fustat) ~ LINC01094+AC147067.2+HAGLR+AC016394.2+AL590705.3+AC023511.1    , data = rt1)
multiCoxSum=summary(multiCox)
multiTab=data.frame()
multiTab=cbind(
             HR=multiCoxSum$conf.int[,"exp(coef)"],
             HR.95L=multiCoxSum$conf.int[,"lower .95"],
             HR.95H=multiCoxSum$conf.int[,"upper .95"],
             pvalue=multiCoxSum$coefficients[,"Pr(>|z|)"])
multiTab=cbind(id=row.names(multiTab),multiTab)
write.table(multiTab, file=multiCoxFile, sep="\t", row.names=F, quote=F)
bioForest(coxFile=multiCoxFile, forestFile=multiCoxPdf, forestCol="red")

}
定义独立预后分析函数

#调用函数,进行独立预后分析
indep(riskFile="risk.all.txt", cliFile="clinical.txt", project="all")

结果错误如下：

#调用函数,进行独立预后分析
indep(riskFile="risk.all.txt", cliFile="clinical.txt", project="all")
Error in coxph(Surv(futime, fustat) ~ rt[, i], data = rt) :
No (non-missing) observations
此外: Warning message:
In max(event[who2]) : max里所有的参数都不存在；回覆-Inf

这个怎么改呀？谢谢了！

数据里没这个变量，楼主自己检查一下数据吧