cox回归生存分析在stata中实现

为了评价各种指标对透析患者生存时间的影响，资料如下，为避免受过大值的影响，实际拟合模型中，ntprobnp数据取对数，记为ntprobnp1 gen ntprobnp1 = log(ntprobnp) 资料如下：

time

die

crea

bun

alb

age

therapy

ntprobnp

ntprobnp1 12 1 1032 34 26 1 1324 46 38 1 1520 49 45 1 870 28 99 1 990 34 120 0 785 23 120 0 456 19 120 0 570 23 120 0 1020 33 120 0 780 29 120 0 670 22 120 0 932 28 120 0 689 27 35 1 670 33 67

1

1210

34

1. 预后因素筛选： logrank time die, by(crea) 检测肌酐是否影响生存率

chi2(13) = 29.07 Pr>chi2 = 0.0064 P<0.05，说明肌酐对预后影响大 logrank time die, by(bun) logrank time die, by(alb) logrank time die, by(age) logrank time die, by(therapy) logrank time die, by(ntprobnp1)

由于bun P>0.05，在COX模型中去除改因素

2. 用COX比例风险模型分析

cox time crea alb age therapy ntprobnp1, dead(die)

Cox regression -- no ties

Entry time 0 ③ Log likelihood = -4.1883032

28 45 26 66 31 56 32 62 39 35 42 29 43 35 39 27 41 42 40 23 39 19 35 23 44 56 28 66 34

72

Number of obs = ① LR chi2(5) = ② Prob > chi2 = ④ Pseudo R2 = 0 20 1 34 0 56 1 112 0 1020 0 1920 0 45 0 88 0 211 0 455 0 2100 0 1320 0 44 1 66 1

77

15 26.21 0.0001 0.7578 2.995732 3.526361 4.025352 4.718499 6.927558 7.560081 3.806663 4.477337 5.351858 6.120297 7.649693 7.185387 3.78419 4.189655 4.343805

------------------------------------------------------------------------------

time | ⑤ Coef. ⑥Std. Err. ⑦z ⑧P>|z| ⑨ [95% Conf. Interval] -------------+----------------------------------------------------------------

crea | -.0032338 .0031113 -1.04 0.299 -.0093319 .0028643 alb | -1.260323 .8536905 -1.48 0.140 -2.933525 .4128797 age | .2428877 .2597176 0.94 0.350 -.2661495 .7519248 therapy | -7.685876 7.077233 -1.09 0.277 -21.557 6.185245 ntprobnp1 | -.2293586 1.013952 -0.23 0.821 -2.216669 1.757952

① 为模型无效假设(即：所有协变量的回归系数为 0) 所对应的似然比检验( 自由度为协变 量个数的卡方); ② 模型无效假设检验对应的 p值；③ 对数似然比；④ 伪决定系数；⑤ 回归系数；⑥ 回归系数的标准误；⑦ 单个回归系数检验 (Ho:该回归系数为0)的Z统计量；⑧ 单个回归系数验的p值；⑨ 回归系数的 95% 可信限。

本例结果表明：在平衡诸多混杂因素后，模型中的因素与死亡风险无明显相关性。

3. 作图

stset time, failure(die) 数据库定义为生存资料数据库；sts list为计算命令。 sts list

sts graph 生存概率曲线

sts graph, by(therapy) 生存概率曲线，不同治疗措施对生存率的影响