Bootstrap内部验证代码，R自带数据集，保证出结果

Bootstrap内部验证代码，R自带数据集，保证出结果

library(rms)

head(lung)

 

用age+sex+ph.ecog+ph.karno这几个指标来预测结局

m <- lrm(status~ age+sex+ph.ecog+ph.karno, x = T,y = T,data = lung)

m

 

res.val <- validate(m, B=1000)

res.val

                                                                                      

c_index <-(res.val[1,1]+1)*0.5

c_index

 

adjust_c_index <-(res.val[1,5]+1)*0.5

adjust_c_index

#根据Dxy计算c index  C=0.5(Dxy+1)

#index.orig是原始Dxy

#index.corrected是校正后的Dxy