9 生存分析.pptVIP

* * 主编 周仁郁 9.1 生存时间 9.1.1 生存函数 生存分析，是把事件的结果和出现这一结果所经历的时间，结合起来分析的一种统计方法． 例1 6名肝癌病人在 1~12月进入观察，其中A，D，E在研究结束前死亡，C由于某种原因失访，B、F到研究结束仍未获得确切的生存时间 生存 180 12 7 F 死亡 120 9 5 E 死亡 120 8 4 D 失访 120 7 3 C 生存 330 12 2 B 死亡 150 5 1 A 结局 生存天数 终止月 开始月 病人 可以看出，A、D、E?在研究结束前死亡，C由于某种原因失访，B、F到研究结束仍未获得确切的生存时间 生存时间是指从某事件开始发生到出现某种反应所经历的时间，可分为完全数据与截尾数据两种类型完全数据是指从起点到死亡所经历的时间，截尾数据，是指由于病人发生迁移或死于其他疾病造成失访、或者改变方案、结束时间尚未发生等截尾情况 生存时间有其特定的分布，描述生存时间的分布规律，可以使用死亡密度函数、生存函数、危险函数等3种方法 死亡密度函数f?(t)表示患者在时间t的瞬间死亡概率 其中，F(t) 称为死亡分布函数。以 t 为横轴、f(t)为纵轴绘制的曲线，称为死亡密度曲线 生存函数S(t)，表示患者生存时间大于t的累积概率 S(t)＝P(T＞t)＝1－F(t) S(t)为t右边死亡密度曲线下面积，也称为生存率。以t为横轴、S(t)为纵轴绘制的曲线，称为生存率曲线或Kaplan－Meier曲线 危险函数h(t)又称风险函数，表示t的瞬间死亡概率 以?t?为横轴、h(t)为纵轴绘制的曲线称危险函数曲线 9.1.2 生存率计算 生存率的计算，有非参数方法与参数方法。参数方法虽然效率较高，但对生存时间的分布要求符合某种类型，而且计算复杂。非参数方法，主要有寿命表与概率乘积两种方法，计算比较容易，应用广泛 概率乘积法也称为乘积极限法，对不分组资料用概率乘法估计生存率，1958年由Kaplan－Meier提出 把生存时间t从小到大排序，重复数据只列一次。分别计算各时段的初始例数n及该时段的死亡例数d，得到死亡概率q及生存概率p q＝d/n，p＝1－q 生存率S(tk)的估计式 S(tk)＝P(T≥t)＝p1?p2?…?pk 生存率S(tk)的标准误为 例2 中药组与对照组两种疗法治疗白血病，估计两组的生存率 40＋ 4＋ 27 27 24 16 10 9 9 4 2 对照组 40＋ 40＋ 36＋ 28＋ 39 37 37 31 22 17 5 中药组 对照组按生存时间t排序，并分段进行计算 在第1时段，t1＝2，存活时间大于等于2例数为11，死亡时间为2的例数为1，得初始例数n1＝11，死亡例数d1＝1，死亡概率q1＝1/11，生存概率p1＝10/11 生存率S(t1)＝p1＝10/11 在第2时段，死亡概率q2=1/10，生存概率p2=9/10 1 1 40 8 0.0969 0.1023 1/3 2/3 2 3 27 7 0.1475 0.3068 3/4 1/4 1 4 24 6 0.1572 0.4091 4/5 1/5 1 5 16 5 0.1599 0.5114 5/6 1/6 1 6 10 4 0.1557 0.6136 6/8 2/8 2 8 9 3 0.1163 0.8182 9/10 1/10 1 1 10 4 2 0.0867 0.9091 10/11 1/11 1 11 2 1 标准误 生存率 生存概率 死亡概率 截尾数 死亡数 初始数 生存月 序号 类似计算中药组的生存率 9.1.3 时序检验 时序检验也称对数秩检验(?Log-Rank Test), 可用于两个或多个样本生存时间的比较，条件是生存率曲线不能交叉。有交叉时采用分层处理或多因素分析 若k（≥2）组的生存时间相同，则在ti点有相同的死亡率di?/?ni，由此算出第j组在ti时点理论死亡例数Eji Eji＝nji?di?/?ni 第j组理论数Ej与实际数Oj构成卡方统计量 理论数有小于5者取修正值 在假设检验有统计学意义时，可以从生存率曲线目测判断、或半数生存期比较、或相对危险度比较几方面来考察效果的好坏。第i组与第j组相对危险度 例3 对例2数据用时序检验比较两组的生存时间 H0：两组的生存时间相等，H1：两组生存时间不等 两组数据混合排序，k＝2，分别计算两组及合计组的n、d，并分段进行计算 在第一时段，n1＝22、d1＝1，中药组n11＝11、E11＝11×1/22＝0.5，对照组E21＝0.5 5.175 10.825 9 7 合计 0 3 1 1 2 2 40 16 0.2500 0.7500