第十七章 生存分析(Survival Analysis)

ti时刻存活的人数 n年生存率 s(ti ) 观察的总人数 90 80 s(1) 0.9 s(2) 0.80 100 100
22
生存时间数据分析时整理示意图
•
• •
失访
死亡 失访 死亡
死亡 死亡 0 2 4 6 8 10 12 14 16 18 20 年
t
23
“t”表示从研究起点到结局出现时间
时间(月) 死亡 ti 人数 期初 人数 死亡 概率 生存 概率 生存率
p( x t )
0.957 0.914 0.783 0.652 0. 609 0.565 0.478 0.478 29
1 3 5 6 7 8 10 14+
1 1 3 3 1 1 2 0
23 22 21 18 15 14 13 11
32
第三节 生存曲线的统计检验
比较不同方法的生存率，常进行生存率曲 线间的比较。 方法：时序检验（Log-Rank test），可对两组 或多组生存率曲线做比较. 检验假设：H0：两总体的生存率曲线相同 H1：两总体的生存率曲线不同 α=0.05，如P≤α，拒绝H0
33
Log-rank检验
检验统计量：
注：生存时间（t）是正偏态分布。
25
Survival Functions
1.2
1.0
甲手术
.8
.6
乙手术
.4
.2
0.0
-.2 0 10 20 30 40 50 60
月
图17-2 两种手术治疗方式术后病人生存曲线的比较 26
（二）生存率估计的统计方法 （非参数方法）
1.小样本数据生存率计算*。 用 kaplan-Meier的乘积极限法（productlimit method，PL法） 方法： 1）将生存时间t由小到大排列。截尾值排在完 全数据后，例：20，20+ 2）列出t时刻死亡数（d） 3）生存率估计用概率乘法原理 例：17-1和表17-1
35
两组生存率曲线的检验 H0：s（t1）= s（t2） 2 2 ( A甲 T甲） （A乙 T乙） 2 T甲 T乙
(16 23.809) (20 12.191) 7.56 23.809 12.191
2 2 2
6.63 p＜0.01 ν=组数-1=2-1， 7.56 结论:两生存率曲线有统计差别, 甲手术方法后生存 率高于乙法.
第十七章 生存分析 (Survival Analysis)
随访研究及统计分析
第二军医大学卫生统计学教研室 孟 虹
1
本章内容
第一节 第二节 第三节 生存分析的基本概念 生存率的估计与生存曲线* 生存曲线的Log-rank检验
第四节
第五节
COX比例风险回归模型*
寿命表(不讲)
2
*要求掌握概念、方法、用途。
（一）随访研究设计 1.明确研究对象的起始事件时间，如手术日期等。 2.明确结局事件：如死亡或复发。 3.明确研究跨度时间：如2000年至2005年结束。 4.记录个体影响结果（y）的其他自变量。
11
例：收集生存数据和影响预后的因素 。
某病不同药后随访记录（天）
预后因素
病例 号 1 2 3 4
随访记录
27
例17-1：某手术方法（甲法）治疗23例肾 上腺肿瘤病人后生存情况（讲义358页）
生存时间（t，月），其中“+”者为截尾数据 1，3，5（3），6（3），7，8，10 （2） ，14+，17，19+，20+，22+，26+， 31+，34，34+，44，59 计算生存率s（t）和生存曲线
28
表17-1 甲种手术后病人生存率的计算方法
（17-1）
生存概率（pi=1-F）：
指某时间段开始存活的个体到该时间段结束 时仍存活的概率。
某时间段存活人数 p 某时间段初的观察数
20
生存率（survival rate）
称为生存函数 记为S（t） S（t） ：指观察对象从起始事件（如手术时间
为0点）开始，到t时刻仍存活的概率。常用n年 生存率表示。
是否 死亡 1 1 0 生存 时间 140 160 99 36
12
性别 处理 开始 终止 组 日期 日期 1 A药 98/07/12 98/11/29 2 B药 98/07/01 98/12/29 1 A药 98/08/22 98/11/29 2 B药
98/10/20 98/11/25 车祸死亡 0
31
寿命表法与PL的区别
1.计算在 (ti 1,ti ) 时间段的生存率。 如0-1年、1-2年，时间段组距相等。 2.寿命表方法计算死亡概率，用校正观察人数计 算。假定有截尾事件的人在各时间组内平均生 存为1/2时间。 （校正观察人数=期初观察人数－截尾人数/2） 死亡概率=某时间组内死亡人数/校正观察人数
（二）随访的方式
1.全部观察对象同时接受不同处理（起点相同） 随访方式：多见于动物实验（见图17-1，a）
2.观察对象在不同时间接受处理因素（起点不同） 随访方式：临床试验研究（见图17-1，b）
13
动物实验随访数据（图17-1，a）
×为死亡
× O O × × 0
O 为截尾
起始事件时间
如给药
t 研究结 束时间
概 述
• 临床上疗效、预后的评价常用疾病 的结局指标：如有效率、治愈率、 死亡率比较。对于短期内能明确治 疗效果的疾病是适用的。但对于远 期疗效，上述指标的评价不全面。
3
例
某病的疗效比较
甲药 乙药 治愈率(%) 80 81 平均治愈时间(月) 20 12
疗效除了应评价“结局”的好坏，结局所 经历时间长短也是评价疗效重要的指标。
16
（三）生存分析主要研究的内容
1.统计描述：计算不同时间点（t）的生存率， 描述生存过程。 2.统计推断：检验不同处理方式的生存过程 有无统计差别. 3.自变量（x）对生存时间（t）的关系：影 响生存时间的危险因素分析.
17
第二节 生存率的估计与生存曲线
（一）描述生存资料的几个指标
1.不同时间点生存率 2.生存曲线
( Ai Ti ) Ti
2
2
该χ2服从自由度=比较组数－1 Ai 为某组各时点实际死亡频数合计. Ti 为某组各时点期望死亡频数合计 i 表示比较组，i=1，2，…k组
34
Log-rank检验的基本思想
表17-4部分数据，365页 时间 甲法手术组 乙法手术组
t 1 2
n1i
d1i
1 0
3.中位生存时间
18
几个率概念
死亡率、死亡概率、生存概率、生存率
死亡率：表示在单位时间（年）内死亡发生 的频率（年平均死亡水平）。
某年内死亡数 某年死亡率 1000% 某年的平均人数
19
死亡概率、生存概率
死亡概率（F）：在某时间段(t)开始存活的个体， 死于(t+△t)该时段内的可能性。
某时间段内的死亡数 d F 某时间段初的观察数 n
2.大样本资料的生存分析方法—寿命表法 （Life-table method）
例17-3 表17-3 2418例男性心绞痛病人生存率情况 术后 死亡 截尾 期初 校正 生存 生存率 年数 人数 人数 人数 人数 概率 （t+1） 0456 0 2418 2418 0.8114 0.8114 1226 39 1962 1942.5 0.8837 0.717 2- 152 22 1697 1686 0.9098 0.6524 校正人数=1962-39/2=1942.5
ti时刻存活的人数( X ) n年生存率 s(t ) 观察的总人数
时间ti ，i=1，2，3…n
假设数据是完全数据，计算见例
21
例:某病病人术后生存率
期初 (ti 1,ti ) 人数 [0，1] 100 [1，2] 90 [2，3] 80 生存 死亡 生存 死亡 生存 生存率 s(ti ) 人数 人数 概率 概率 10 90 0.1 0.9 0.90 10 80 0.11 0.89 0.80 20 60 0.25 0.75 0.60
4
例 2：
两种方法对疾病的疗效 方法 治疗人数 生存人数 生存率% 甲方法 100 20 20 乙方法 100 50 50 经χ2检验 p＜0.05，乙法预后优于甲法。 假定： 1.观察期间疾病的死亡率不随时间变化。 2.研究对象观察时间长度相等。
5
随访研究
随访研究(follow-up study)是医学中常用前 瞻性研究. 例:两种方法肾移植病人术后肾的生存时间 和结局(生存率)比较. 例:不同方法对某病人(癌症、反复发作疾病） 生存时间(缓解时间)与结局(生存率)比较. 该类数据通过随访得到，称为随访资料。
9
3.生存数据的特点
1）完全数据：研究对象在规定研究期间提供确 切的“时间和结局”。
2）截尾数据(t+) ：截尾数据虽然提供的信息不 完全，但提供了部分信息，如 t=10+年＞9年。 3）生存数据的结果变量（Y ）有两个： 时间（t）值 ，t＞0 结局状态（y ）=“ 如死亡或截尾值”
10
二 资料的收集
1.分析的变量（y） 1） 结局事件：指结局出现的特征，如疾病的死 亡、复发、发生( y=1或0) 。 2）时间间隔变量 记为（t）
t=结局事件出现日期 － 事件的起始日期
（起始日期可规定:如诊断、用药、手术日期 8 等）, t的单位：可用年、月、周表示
2.截尾数据
观察过程中个体因其他原因未观察到明确的 结局， 称为截尾或删失数据（ censored data）。 截尾原因有： ①失访，②退出研究，如其他原因死亡。③研究 时间结束，未出现结局事件。 截尾值（censored value）： 时间（t）=截尾事件日期 －起始事件日期 记为t+。（例：10+月）