// 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