卫生统计学第十五章生存分析.ppt

生存分析 Survival Analysis 一、什么是生存分析？ 在医学研究中，常常用追踪(follow up)的方式来研究事物发展的规律。 了解某药物的疗效 了解手术后的存活时间 了解某医疗仪器设备的使用寿命 这种研究的特点是追踪研究的现象都要经过一段时间，统计学上将这段时间称为生存时间。 生存分析就是用来研究生存时间的分布规律以及生存时间和相关因素之间关系的一种统计分析方法。 二、生存分析的基本概念 生存时间 (survival time)： 从某起始事件起到某终止事件止所经历的时间跨度称为生存时间。例如， 在临床研究中，急性白血病患者从骨髓移植治疗开始到复发为止之间的时间间隔。 冠心病患者在两次发作之间的时间间隔。 在流行病学研究中，从开始接触某危险因素到发病所经历的时间。 在动物实验研究中，从开始给药到发生死亡所经历的时间。 生存时间数据的类型 完全数据(complete data)：它提供的关于生存时间的信息是完整确切的，也就是说它准确地度量了观察对象实际生存的时间。 截尾数据(censored data)：它提供的关于生存时间的信息是不完整不确切的，也就是说它没有准确地度量观察对象实际生存的时间。例如， 在随访过程中某些观察对象失访； 或死于其它原因； 或在规定的研究过程结束时观察对象的终止事件还未发生。 生存时间数据的特点： 所有观察值取值非负 存在截尾数据 生存时间函数： 描述生存时间分布规律的函数，例如， 生存函数、死亡函数、死亡密度函数、风险函数 生存分析的主要任务： 描述生存过程 比较生存过程 分析危险因素 建立数学模型 生存函数S(t)：观察对象的生存时间T大于某时刻 t 的概率称为生存函数(survival function)。S(t)满足条件：S(0)=1， S(∞)=0，且0≤S(t)≤1。 三、生存时间函数的估计 死亡函数F(t)：观察对象的生存时间T不大于某时刻 t 的概率称为死亡函数(failure function)。F(t)满足条件：F(0)=0， F(∞)=1，且0≤F(t)≤1。 死亡密度函数f(t)：观察对象在某时刻t 的瞬时死亡率称为死亡密度函数(failure density function)。 风险函数h(t)：生存到时刻 t 的观察对象在时刻t 的瞬时死亡率称为风险函数(hazard function)。 四、生存分析的基本方法 描述法：根据样本观察值提供的信息，直接用公式计算出每一时间点或每一个时间区间上的生存函数、死亡函数、风险函数等，并采用列表或绘图的形式来显示生存时间的分布规律。 优点：方法简单且对数据的分布无要求。 缺点： 1)不能比较两组或多组生存时间分布函数的区别； 2)不能分析危险因素对生存时间的影响； 3)不能建立生存时间与危险因素之间的关系模型。 Survival Function Estimates 1.0 + ** | ** SDF | **N | L*N | L*NN 0.5 + LLNNN | LL NNNN | LLL NNN | L--L N-N--N | | N-------------N---N 0.0 + L N ----+----+----+----+----+----+----+----+----+----+----+----- 0 10 20 30 40 50 60 70 80 90 100 Time 无淋巴结浸润患者(N)比有淋巴结浸润患者(L)的平均生存时间长。 例如，两组急性淋巴细胞性白血病患者治疗后的随访资料: 非参数法：估计生存函数时对生存时间的分布没有要求，并且检验危险因素对生存时间的影响时采用的是非参数检验方法。 优点： 1)可以估计生存函数； 2)可以比较两组或多组生存布函数； 3)可以分析危险因素对生存时间的影响； 4)对生存时间的分布没有要求。 缺点：不能建立生存时间与危险因素之间依存关系的数学模型。 参数法：根据样本观察值来估计假定的分布模型中的参数，获得生存时间的概率分布模型。 优点： 1)可以估计生存函数； 2)可以比较两组或多组生存布函数； 3)可以分析危险因素对生存时间的影响； 4)可以建立生存时间与危险因素之间依存关系的模型。 缺点：需要事先知道生存时间的