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

 

Logo

魔乐社区(Modelers.cn) 是一个中立、公益的人工智能社区,提供人工智能工具、模型、数据的托管、展示与应用协同服务,为人工智能开发及爱好者搭建开放的学习交流平台。社区通过理事会方式运作,由全产业链共同建设、共同运营、共同享有,推动国产AI生态繁荣发展。

更多推荐