这篇博客介绍了如何在STATA中进行C指数的假设检验，以比较不同预测模型（如TNM分期）的性能。首先，通过安装'somersd'包，然后设置生存时间和失效变量，接着进行Cox回归调整其他变量。通过'estat concordance'命令计算C指数，并预测风险比。使用'somersd'函数计算C指数并进行转换。最后，运用'lincom'命令进行C指数间的统计比较，评估模型间的差异。
// STATA的 STATA的 STATA的
//C指数间差异的假设检验
//我自己看的文章，常见的就是一个指标和TNM分期的性能作比较
//安装
ssc instal somersd
//定义模型结果和时间变量
stset time, failure(dead ==1)
//进行cox回归 调整其他变量
stcox xxx agegroup i.bmi i.TNM
estat concordance
predict hr1
generate invhr1=1/hr1
generate censind1=1-_d if _st==1
somersd _t invhr1 if _st==1 , cenind(censind1) tdist transf(c)
stcox i.TNM
estat concordance
predict hr2
generate invhr2=1/hr2
generate censind2=1-_d if _st==1
somersd _t invhr2 if _st==1 , cenind(censind2) tdist transf(c)
somersd _t invhr1 invhr2 if _st==1 , cenind(censind1) tdist transf(c)
//统计比较
lincom invhr1-invhr2
7991
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
