方法一: 方法二: library(survival) cindex <-concordance.index(predict(fit),surv.time =result$os,surv.event = result$event1,method = 'noether') > cindex$c.index; cindex$lower; cindex$upper[1] [1] 0.3010949[1] 0.4181514-concordance.index(predict(fit),surv.time =result$os,surv.event = result$event1,method = 'noether') > 方法三,为了验证,我换一种办法: > survConcordance(Surv(result$os, result$Event1) ~ predict(fit, result))Call:survConcordance(formula = Surv(result$os, result$Event1) ~ predict(fit, result)) n= 449 concordant discordant tied.risk tied.time std(c-d) 21313.000 6483.000 13485.000 7.000 1925.018 方法四: library(survival) > library(Hmisc)> fp <- predict(fit)#模型预测> cindex=1-rcorr.cens(fp,Surv(result$os, result$Event1)) [[1]]# > cindex[1] 0.6796226- predict(fit)#模型预测> 方法一,三四都是0.68.而方法二是0.30.说实在的,真不知道 方法二会有这么大差异,为何会不一致。你会认为方法二不对吗? |
|