生存分析 (卫管 2011)

《生存分析 (卫管 2011)》由会员分享，可在线阅读，更多相关《生存分析 (卫管 2011)（64页珍藏版）》请在金锄头文库上搜索。

1、2018/10/18,1,黄志刚 广东医学院公共卫生学院,第十四章 生存分析,2018/10/18,2,在临床医学实践中, 对病人疗效考查： 治疗结局？ 生存时间？,一、引言,2018/10/18,3,【例14-1】某临床医生将22例肺癌患者随机分成两组，分别采用化疗和放化疗联合治疗，从缓解出院日开始随访，随访时间（月）如下（带“+”号的数据表示患者至少活了多少个月）。 化疗组：1，2，3，5，6，9+，11， 13，16，26，37+ 放化疗联合组：10，11+，14，18， 22， 22，26 ，32，38，40+，42+ 欲比较两种治疗措施的疗效应采用何种统计分析方法？,2018/10/

2、18,4,化疗组：1，2，3，5，6，9 ，11， 13，16，26，37 （月） 放化疗联合组：10，11 ，14，18， 22， 22，26 ，32，38，40 ，42 （月）,+,+,+,+,+,2018/10/18,5,二、生存资料的特点,蕴涵有结局和时间两个方面的信息 结局为两分类互斥事件 一般是通过随访收集得到，随访观察往往是从某统一时间点开始，观察到某规定时间点截止 常因失访等原因造成某些研究对象的时间数据不完整，分布类型复杂，一般不能采用t test、ANOVA、卡方检验等进行统计分析,2018/10/18,6,随访（follow-up）资料的记录,生存资料一般通过随访收集 记

3、录的项目： （起始与终止）随访事件 生存时间（开始观察日期与终止观察日期）（年、月、天、小时、分、秒等） 分组变量（处理方法） 和其它变量或因素（性别、年龄、职业、文化程度等）,2018/10/18,7,三、生存分析的概念,生存分析(survival analysis)是将事件的结局和出现这一结局所经历的时间，结合起来分析的一大类统计方法 近30多年来发展起来的一个新的数理统计学分支，是一个相对独立的理论体系 70年代中期以来应用广泛，目前在工业、农业、国防和医学等领域均有广泛的应用,2018/10/18,8,四、生存分析的应用,描述生存过程 研究生存时间的分布特点，估计生存率及绘制生存率曲线

4、 方法：乘积极限法(kaplan-meier法) 寿命表法 (life-table method),2018/10/18,9,比较生存过程 进行两组或多组生存曲线（生存过程）比较 方法:对数秩检验（log-rank test） 生存过程的影响因素分析 采用多因素生存分析方法寻找影响预后的主要因素，达到在手术前后加以预防和干预的目的 方法：Cox比例风险模型,2018/10/18,10,五、生存分析中的几个基本概念,死亡事件(death event)：又称失效事件或终点事件，不单是指通常意义下的生物体死亡，而是泛指标志某种处理措施失败、失效的特征事件，或者出现了预期的结局事件。 一般是在设计阶段

5、根据研究目的确定。,2018/10/18,11,生存时间(survival time)：指观察到的“存活”时间。不单是指通常意义下生物体的存活时间，而是泛指研究者所关心的某现象的持续时间。 可用天、周、月、年等时间单位记录，常用符号t 表示 根据其特点，可分为两种类型： 完全数据 截尾数据,2018/10/18,12,完全数据(complete data)：是指从观察起点到发生死亡事件所经历的时间 截尾数据(censored data)：简称截尾值 (censored value)，又称删失值或终检值。观察过程的终止不是由于死亡事件，而是由于其他原因。常用t +表示,失访退出（死于其它原因）终

6、止（截止随访结束时尚存活）,2018/10/18,13,临床随访研究中的完全数据与截尾数据,观察起点,观察终点, 起点事件 终点事件失访退出终止,2018/10/18,14,起始事件 生存时间 终点事件（死亡事件） 疾病确诊 死亡 疾病确诊 痊愈 治疗开始 死亡 治疗开始 痊愈 症状缓解 疾病恶化 接触毒物 出现毒性反映 接触危险因素 发病,2018/10/18,15,生存时间资料的特点,2个效应变量 （1）生存时间(天、月、年等) （2）结局(死亡与否、是否复发等) 截尾数据（病人未观察到底，不知他们究竟能活多长时间）产生结尾原因：(1)迁移 (2)死于其他原因 (3)因其他客观原因中途退出

7、(4)预定终止结果迟迟不发生 分布类型复杂：生存时间分布常呈正偏态分布,2018/10/18,16,两种错误的做法： 错误1：丢弃截尾数据，只考虑确切数据（损失了信息） 错误2：将截尾数据当作确切数据处理（低估了生存时间的平均水平）,截尾数据 (censored data),2018/10/18,17,两种错误的做法： 错误1：采用平均生存时间而不是采用中位生存时间来表示生存时间的平均水平。 错误2：采用常规 t 检验或方差分析进行组间比较。（应采用log-rank检验比较几组生存时间 ）,正偏态（positive skewness）数据,2018/10/18,18,生存分析的统计描述指标,1

8、、死亡概率、生存概率 2、生存率及其标准误 3、生存曲线 4、半数生存期（中位数） 及四分位数间距,2018/10/18,19,死亡概率（mortality probability）：指在某单位时段开始时存活的个体在该时段内死亡的可能性大小，记为q 年死亡概率的计算公式为若年内有截尾，则分母用校正人口数,2018/10/18,20,生存概率(survival probability)：表示在某单位时段开始时存活的个体，到该时段结束时仍存活的可能性大小，记为p 某年生存概率的计算公式为,2018/10/18,21,生存率(survival rate) 记为 ，指观察对象经历 t 个单位时段后仍存

9、活的可能性，即活过 tk时刻的概率。 若无截尾数据，则,2018/10/18,22,如果含有截尾数据，采用概率乘法原理估计 假定观察对象在各个时段的生存事件独立，生存概率分别为 ，则根据概率乘法原理得到估计生存率的计算公式为,实质上是累积生存概率(cumulative probability of survival),2018/10/18,23,生存率 ( 累积生存概率)指病人经历t个单位时间后仍存活的概率生存概率指单个时段的概率 生存率指从0t多个时段的累积概率,2018/10/18,24,例 手术治疗50例肺癌病人，术后1，2，3年的 死亡数分别为10，10，10例，无截尾数据。 试求各年

10、的生存概率和3年生存率。各年生存概率： p1 = ( 50 10 ) / 50=40/50p2 = ( 40 10 ) / 40 =30/40 p3 = ( 30 10 ) / 30=20/303 年生存率:直接法： S(3) = ( 50 30 ) / 50 = 0. 4概率乘法定理:或 S(3) = p1 p2 p3 = 40/5030/4020/30=0.4,2018/10/18,25,生存率的标准误(Greenwood估计),nj 表示时刻 tj 的期初观察人数 dj 表示 tj 时刻的死亡人数,2018/10/18,26,生存率的可信区间估计(正态近似法),100(1-)%可信区间为

11、：,2018/10/18,27,生存曲线(survival curve)：以观察（随访）时间为横轴，以生存率为纵轴，将各个时间点所对应的生存率连接在一起的曲线图，用以描述生存过程。,2018/10/18,28,生存曲线（K-M法）,阶梯状：每一级阶梯代表一个死亡时间点 (在截尾时间点无阶梯)；如果最大时间点是截尾则生存曲线不与横轴相交（见下图），否则与横轴相交。,Death,Censored,2018/10/18,29,半数生存期及四分位数间距半数生存期也称中位生存期（median survival time) 即生存时间中位数，表示50%的个体可存活的时间，即生存率为50%时对应的生存时间（

12、集中趋势指标） 生存期的四分位数间距 2575，反映离散程度大小,2018/10/18,30,Median Survival（中位生存期）,1.8 Years,2018/10/18,31,描述生存过程 乘积极限法（kaplan-meier法）适用于小样本的未分组生存时间资料 寿命表法 （life-table method ）适用于大样本的分组生存时间资料 比较生存过程：对数秩检验（log-rank test） 生存时间的影响因素分析：Cox 回归分析 生存率的预测： Cox 回归模型,生存分析的基本研究方法,2018/10/18,32,未分组资料是指由每个观察单位的原始测得值组成的资料，在生存

13、分析中即是各观察单位的具体生存时间资料，常用乘积极限法(product-limit estimate) 乘积极限法是一种非参数法，由Kaplan-Meier于1958年首先提出，故又称Kaplan-Meier法(简记为KM法),六、未分组资料的生存分析,2018/10/18,33,【例14-2】 某医师采用手术疗法治疗12例宫颈癌患者，随访时间（月）记录如下：1，2，4，5，7，8+，11，15，18，33+，36，38+。试估计各时点生存率及其标准误、各时点总体生存率的95%可信区间、中位生存时间，并绘制生存曲线。,2018/10/18,34,Kaplan-Meier法基本思想： 将生存时间

14、由小到大依次排列，在每个死亡点上，计算其期初人数、死亡人数、死亡概率、生存概率和生存率。,2018/10/18,35,2月的死亡人数,在2月 时点以前尚存活的病例数,2月时点以前尚存活的患者在2月死亡的概率,2018/10/18,36,2月时点以前尚存活的患者在2月存活的概率（单位时段）,2月时点以前尚存活的患者活过2月概率,2018/10/18,37,2018/10/18,38,生存率的标准误(Greenwood估计),nj 表示时刻 tj 的期初观察人数 dj 表示 tj 时刻的死亡人数,2018/10/18,39,生存率的可信区间估计(正态近似法),100(1-)%可信区间为：,2018

15、/10/18,40,中位生存时间的计算（内插法）,2018/10/18,41,生存曲线 未分组资料的生存曲线也称Kaplan-Meier曲线。它是以生存时间t为横轴，生存率为纵轴，绘制而成的连续的阶梯形曲线，用以说明生存时间与生存率之间的关系 生存曲线图可直观地比较各组观察对象的生存过程，也可对任意时刻的生存率作出粗略估计,2018/10/18,42,38+,36,18,15,11,7,5,4,2,1,33+,8+,2018/10/18,43,七、分组资料的生存分析,随访结果只有某年或某月的观察人数、发生死亡事件人数和截尾人数，而没有每个观察对象确切的生存时间，只能获得按随访时间分组的资料 当

16、样本含量较大（如n50）时，采用乘积极限法估计其生存率及其标准误较为繁琐，往往也是先将原始资料按照生存时间分组后再进行分析 分组资料的生存分析常用采用寿命表法(lifetable method),2018/10/18,44,【例14-3】 某医师对110例原发性肺癌患者确诊后进行随访，得生存资料见表14-3第(2)栏至第(4)栏。试估计各时点生存率及其标准误、各时点总体生存率的95%可信区间、中位生存时间，并绘制生存曲线。,2018/10/18,45,2018/10/18,46,2018/10/18,47,2018/10/18,48,2018/10/18,49,2018/10/18,50,2018/10/18,51,