生存分析Stata实现

第十四章生存分析的Stata实现

本章使用的STATA命令：

结局变量为1表示失效事件发生

例14-2 McKelvey et al(1976)收集了3期的某型淋巴瘤患者的生存时间（天）。分别是6，19，32，42，42，43+，94，126+，169+，207，211+，227+，253，255+，270+，310+，316+，335+，346+。现用Kaplan-Meier法计算生存率。

Stata命令为：

stset time,failure(d) sts list

sts graph

结果为：

例14－3 下面是来自于Berkson & Gage(1950)的一个研究队列。为了叙述方便，把原来的出院后的生存时间改称为某恶性肿瘤术后生存时间。共有374名患者进入研究队列。

表14－3 寿命表法计算生存率的计算用表

(1) (2) (3) (4) (5) (6) (7) (8) (9)

序号术后生存

年数

期初观察

例数

期内

死亡

期内截尾

人数

校正期初

人数

死亡

概率

生存

概率

生存率t n d c n c=n-c/2 q=d/n c p=1-q S(t)

1 0～374 90 0 374 0.2406 0.7594 0.7594

2 1～284 76 0 284 0.2676 0.7324 0.5561

3 2～208 51 0 208 0.2452 0.7548 0.4198

4 3～157 2

5 12 151 0.165

6 0.8344 0.3503

5 4～120 20 5 117.5 0.1702 0.8298 0.2907

6 5～95

7 9 90.5 0.0773 0.9227 0.2682

7 6～79 4 9 74.5 0.0537 0.9463 0.2538

8 7～66 1 3 64.5 0.0155 0.9845 0.2498

9 8～62 3 5 59.5 0.0504 0.9496 0.2372

10 9～54 2 5 51.5 0.0388 0.9612 0.2280

11 10＋47 21 26 34 0.6176 0.3824 0.0872 Stata数据格式为：

Stata命令为：

stset time [weight=w],failure(d=1) ltable time d [weight=w],graph

结果为：

例14－4 在例14－3中介绍了McKelvey et al收集了19名3期淋巴瘤病人的生存资料，该研究者还同时收集了4期淋巴瘤病人的生存数据共61例：4，6，10，11，11，11，13，17，20，20，21，22，24，24，29，30，30，31，33，34，35，39，40，41+，43+，45，46，50，56，61+，61+，63，68，82，85，88，89，90，93，104，110，134，137，160+，169，171，173，175，184，201，222，235+，247+，260+，284+，290+，291+，302+，304+，341+，345+。要比较3期淋巴瘤患者和4期淋巴瘤患者的生存率是否有差异。

Stata命令为：

stset time,failure(d=1) sts test group

sts graph,by(group) 结果：