生存分析 shihy 医学课件.ppt

另一部分是参数部分，可以通过样本的实际观察值进行估计。 多因素分析中的样本含量一般应在自变量个数的10倍以上。 在医学研究中，观察对象的生存时间的长短往往与多种因素有关系。如宫颈癌患者术后的生存期、白血病患者化疗后的缓解期等，除了跟治疗方案有关之外，还可能与患者年龄、体质、病情轻重及营养状况等有关。 医学上将这些因素称为预后因素，统计学上称为协变量。 由于生存时间资料常存在截尾值，不满足正态分布和方差齐性的要求?不适宜用多元线性回归 如果不考虑生存时间仅考虑事件结局，或将生存时间转换为二分类变量，显然会损失信息。?不适宜用logistic回归 Cox模型：由英国统计学家D.R.Cox于1972年提出的比例危险率回归模型。 1、Cox模型的用途 用于分析带有协变量的生存时间资料； 医学上主要用于探讨恶性肿瘤和其它慢性病的预后因素，也可用于临床疗效评价和病因探索。 Cox模型的基本思想 是在风险函数与研究因子之间建立类似于广义线性模型的关联，这样就可以直接考察研究因子对风险函数的影响效果。 2、Cox模型的结构 属于半参数模型 基准风险函数：表示当所有预后因素处于0时的风险函数。函数形式无任何限定，具有非参数的特点。 ?j的流行病学含义：预后因素每改变一个单位时所引起的相对风险度改变量的自然对数值。 3、参数估计和假设检验 参数估计值b的计算：极大似然函数法 参数估计值b的假设检验： Score检验 Wald卡方检验 似然比检验 案 例 胰腺癌手术中接受放射性治疗是否会延长病人生存时间的研究数据。（data：pancer.sav） 因素 变量名 赋值说明 患者编号 caseno 生存时间（月） time 删失 censor 死亡=1，删失=0 手术时的年龄 Age 处理组别 trt 无术中放疗=0，有术中放疗=1 性别 sex 男=0，女=1 占位处 bui 胰脏头部=0，头部以外=1 TNM分类 stage III期=3，IV期=4 胰胆管浸润程度 ch 1-4 有无腹膜转移 p 无转移=0，有转移=1 胰腺癌术中放疗效果的研究中变量赋值情况 由于该研究是一项未经随机化的观察性研究，要正确估计术中接受放射性治疗提高患者生存时间的效果，还需要考虑对其它因子的效果进行评价。 在其它因子固定不变的情况下，平均来说，在任何一个时间点上，手术中接受放射治疗的患者的死亡风险是未接受患者的0.441倍。 注意事项 以科学的设计、有代表性的抽样为前提 必须满足PH条件：死亡的风险与其基础风险在所有生存时间点上都保持一个恒定的比例。当生存曲线有明显交叉时，说明可能有混杂因素，不宜使用。 样本含量：一般在40例以上；自变量的5-20倍 先进行单因素分析，选择单因素分析有意义的变量进行cox回归。 多重共线性 局限性：在估计参数时不是利用精确的生存时间，而是利用生存时间的顺序统计量，损失了一定的样本信息。 提问小结 什么是生存分析？ 什么是生存时间？ 生存资料的特点 反应变量？ 描述平均生存时间应选用什么指标？ 截尾值 提问小结 生存率的估计----乘积极限法和寿命表法分别适用于什么情况？ 生存率的标准误表示？ 生存曲线的纵横轴分别是？ 乘积极限法中，绘制的生存曲线为阶梯型曲线。 （ ） 生存曲线的比较方法是？其统计量是？ 提问小结 含有删失值的生存时间比较可以用t检验，也可以用log-rank检验。 （ ） log rank检验、乘积极限法都属于非参数法。 （ ） Cox回归模型的因变量是？与多重线性回归、logistic回归有何不同？ Cox回归模型属于半参数模型。 （ ） 某医生随访搜集了200例食管癌术后病人的生存时间随访资料，欲分析多个预后因素对生存时间的影响，宜选用（ ）方法 案例讨论 案 例1 22例Ⅲ期非小细胞肺癌患者在不同日期经随机化分配到放疗组和放化疗联合组，从缓解出院日开始随访，随访时间(月)如下： 放疗组 1，2，3，5，6，9+，11，13，16，26，37+ 放化疗联合组 10，11+，14，18，22，22，26，32，38，40+，42+ 案例2：56例某种癌症病人在不同时期经随机化分配到A、B两治疗组，并继续进行随访至1974年5月 31日结束。资料如下表： 问： 判断上述随访时间哪些属截尾值，写出相应的病人号。 要比较A、B疗法对该种癌症病人的疗效，宜选用何种统计检验方法？ 一般用何指标表示A、B治疗组生存时间的长短？ 案 例3 应用组织芯片研究人肾癌组织中Survivin、EGFR和PCNA表达与预后的关系 Data：kidney.sav EGFR表达阳性组与阴性组比较 The end ！ 博客地址：/statistics Logrank检验：属于单因素分析方法 同一条起跑线 右