随访资料的生存分析一、生存分析的意义原理及方法在其他领域的应用利用生存分析的原理和方法，把“出生”、“死亡”的涵义稍加变通，在医学研究中可以有广泛的应用。应用举例：观察某“手术或化疗”对恶性肿瘤病人的疗效，可将“手术或化疗”代替“出生”，用生存分析来分析接受该处理患者的生存时间以及活过某时点的概率。以第一次心肌梗死代替“出生”，以第二次心肌梗死代替“死亡”，用生存分析可以预计在多少年或月内发生第二次心肌梗死的概率。二、生存分析的基本概念（一）生存分析的基本术语1、‘‘死亡’’事件或称失败事件(failureevent)：表示观察到随访对象出现了我们所规定的结局，是反映处理因素失败或失效的特征。如乳腺癌病人手术后复发、肾移植病人肾功能衰竭、白血病患者化疗后的复发等。失败事件的认定是生存分析的基石，必须绝对准确。注意：失效事件应当由研究目的而定，并非一定是死亡，而死亡也并非一定是失败事件。2、截尾值(censoredvalue)：删失值。有的观察对象终止随访不是由于失败事件发生，而是由于中途失访、死于其它原因、随访截止。由于不知道这些观察对象发生失败事件的时间，他们的资料不能提供完全的信息，这些对象的观察值称为截尾值，常用符号“+”表示。如140+天。3、生存时间(survivaltime)即随访观察持续的实足时间，按失败事件发生或失访前最后一次的随访时间记录。按天、周、月、年等时间单位记录，常用符号t表示。一般情况下较细的时间单位准确性较高，应尽量以个体为单位采用较细的时间单位来记录。但在许多大型的随访中，不可能做到按个体记录，常见的是按固定时间段（如一月一次等）记录有多少人失访及多少人发生失败事件，此为分组生存资料。如某病人1990年2月1日进入随访，1992年4月间发生失败事件，他的生存时间为t=26月；又如某白血病患者化疗3月后失去联系，他的随访结果为一截尾值，生存时间记为t=3+月；安放心脏起搏器患者术后2年因意外事故死亡，他的随访结果也为一截尾值，生存时间记为t=2+年。某医院对100例原发性肝癌生存情况随访结果4、生存率(survivalrate)实际应当是生存概率，指某个观察对象活过t时刻的概率，常用P(X>t)表示。如P(X>10)表示某对象活过10天(或10月、10年)的概率。根据不同随访资料的失败事件，生存率可以是缓解率、有效率等。（二）随访内容随访结局和终止随访时间记录影响生存的有关因素随访资料数据的特点：（1）应变量有两个，即生存时间（天数）和结局（死亡与否）。（2）生存时间存在观察不完全的数据。（三）随访方式“×”——表示“死亡”；“o”——表示失访退出研究或死于与本处理无关的其它原因。（四）生存分析的主要内容和基本方法2、比较生存过程两组或多组生存曲线比较。常用方法：对数秩检验、Gehan比分检验、Breslow检验。3、生存过程的影响因素分析比较不同亚人群的生存状况，进行两组或多组生存率比较，以了解哪些因素会影响目标人群的生存过程，这是生存分析方法最重要的研究内容，在临床医学中应用非常广泛。例如分析影响乳腺癌病人手术后预后的因素，可以是病人的年龄、病程、术前健康状况、有无淋巴结转移、术后有无感染、辅助治疗措施、营养等。Cox比例风险回归模型。三、生存资料的统计描述和生存率的区间估计0.103（1）死亡数（dx）：与生存时间对应的发生死亡时间的人数。（2）期初病例数（nx）：恰好在X时刻以前的病例数。（3）条件死亡概率（qx）：恰好在X月前的观察对象在X月时点死亡的概率。（4）条件生存概率：恰好在X月前的观察对象活过X月时点的概率。（5）生存率P(X＞t)：即观察对象活过某时点的概率。某时点的生存率为≤t时刻的各时点条件生存率的乘积。2、生存率曲线生存率曲线阶梯曲线（Kaplan-Meier曲线）由图可直观地比较各样本的生存率曲线，也可对某一病例任意时刻的生存率作出估计，反之亦可由任意生存率估计生存时间。本例中位生存时间为11个月。中位生存时间常用于比较随访资料。3、总体生存率可信区间的估计（二）分组资料的生存分析分组资料的生存率估计使用寿命表(lifetable)法，这也是一种非参数统计方法。例某医院对100例原发性肝癌患者确诊后进行随访，得资料见下表，根据此资料进行分组资料生存率和标准误的计算。生存率曲线总体生存率可信区间的估计四、两样本生存率曲线的比较2、两样本生存曲线比较——对数秩检验(log-ranktest)其基本思想是如果无效假设成立，即两总体生存曲线无差别时，则两种处理在各个时期的实际死亡数与理论死亡数不会相差太大，否则应认为无效假设不可能成立，两条生存率曲线差异有统计学意义。对数秩检验(log-ranktest)例22例Ⅲ期非