生存分析与Cox回归.pptx

《生存分析与Cox回归.pptx》由会员分享，可在线阅读，更多相关《生存分析与Cox回归.pptx（81页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、2023/3/27生存分析与Cox回归1生存分析与Cox回归生存资料概述生存分析的基本概念生存资料的统计描述生存曲线的比较Cox回归第1页/共81页2023/3/27生存分析与Cox回归2一、生存资料概述 举例 某医师分别用中药、西药各治疗急性肝炎病人40例，结果如下表，试问：哪种药物的治疗效果好？第2页/共81页2023/3/27生存分析与Cox回归3一、生存资料概述随访研究(follow up)是医学研究中常用的设计方法随访研究不仅要考虑观察对象的结局，还要考虑出现结局所经历的时间这类(既要考虑结局又要考虑结局出现的时间)资料称为生存资料(survival data)生存分析(surviv

2、al analysis)是将观察的结局和出现结局所经历的的时间结合起来进行分析的统计方法。第3页/共81页2023/3/27生存分析与Cox回归4一、生存资料概述医学随访研究一般有两种所有研究对象同时进入研究(观察起始时间相同)例如，队列研究、动物的随访观察被研究对象逐个进入研究(观察起始时间不同)例如，临床随访研究由于受经费和时间的限制，最终观察时间不能无限延长第4页/共81页第5页/共81页第6页/共81页2023/3/27生存分析与Cox回归71.基本概念 起始事件与终点事件起始事件反映研究对象生存过程的起始特征事件。终点事件(outcome event)又称失效事件(failure e

3、vent)或死亡事件(death event)终点事件研究者所关心的研究对象的特定结局，可以标志某种处理措施失败或失效的特征事件起始事件和终点事件是由研究目的决定的，在设计时就明确规定，并在研究期间严格遵守不能随意改变第7页/共81页2023/3/27生存分析与Cox回归81.基本概念 举例 说明下列研究的起始事件与终点事件急性白血病患者进行骨髓移植后以是否复发来评价骨髓移植效果职业性铅中毒的危险因素(开始职业性接触至出现铅中毒症状)冠心病患者两次发病的时间间隔大肠癌患者手术后存活情况(手术、死亡)接受健康教育对青少年戒烟到复吸的影响因素分析接受某种保险方式后的中途退保分析第8页/共81页20

4、23/3/27生存分析与Cox回归91.基本概念 生存时间(survival time)或失效时间(failure time)生存时间指观察到的存活时间常用符号 t 表示生存时间是生存分析中的重要信息，必须准确明确规定起始事件、终点事件时间的测度单位(年、月、日)第9页/共81页2023/3/27生存分析与Cox回归101.基本概念 整个研究的观察时间研究开始到研究结束的时间因为有起始事件发生时间、终点事件发生时间、观察开始时间、观察结束时间，生存资料数据分为完全数据(complete data)和截尾数据(censored data)第10页/共81页2023/3/27生存分析与Cox回归1

5、11.基本概念完全数据：指从观察起点到发生死亡事件所经历的时间，生存时间是完整确切的。截尾数据或截尾值(censored value)，又称删失值或终检值。生存时间观察过程的截止不是由于死亡/终点事件，而是由于其他原因引起的，称为截尾(censored)从观察起点到截尾时点所经历的生存时间称为截尾数据，习惯上在生存时间右上角标注“”表示第11页/共81页2023/3/27生存分析与Cox回归121.基本概念 截尾的主要原因：失访(withdrawal):失去联系，如信访无回音、电话采访不应答、上门采访找不到人、搬迁没留地址等 退出:死于非研究因素或非处理因素而退出研究，如死于意外或其他疾病 终

6、止:设计时规定的研究时限已到而终止观察，但研究对象仍然存活第12页/共81页第13页/共81页第14页/共81页2023/3/27生存分析与Cox回归151.基本概念完全数据提供了观察对象确切的生存时间，是生存分析的主要依据；截尾数据仅提供了部分信息，研究者并不知道观察对象确切的生存时间。截尾数据太多会影响生存分析的效果那么截尾数据能不能删除？第15页/共81页2023/3/27生存分析与Cox回归161.基本概念 对生存资料的两种错误分析抛弃截尾数据，只考虑确切数据 损失样本含量、损失了信息，截尾数据提供部分信息，说明在某时刻之前仍存活将截尾数据当作确切数据处理 低估了生存时间的平均水平，截

7、尾数据中存在生存时间较长的数据，如果损失，会使结果产生偏性第16页/共81页2023/3/27生存分析与Cox回归171.基本概念对截尾数据进行分析是生存分析的重要特点，在生存分析中，结局变量常以1表示出现结局，0表示截尾。应变量有2个:生存时间 t 和结局变量(0-1)第17页/共81页2023/3/27生存分析与Cox回归182.生存资料的特点 蕴涵有结局和时间两个方面的信息；结局为两分类互斥事件；一般是通过随访收集得到，随访观察往往是从某统一时间点(如确诊、入院或实施手术等某种处理措施后)开始，观察到某规定时间点截止；常因失访等原因造成某些研究对象的生存时间数据不完整；分布类型复杂，需用

8、生存分析第18页/共81页2023/3/27生存分析与Cox回归193.生存资料的数据形式第19页/共81页2023/3/27生存分析与Cox回归204.生存分析的基本内容 描述生存过程：研究生存时间的分布特点，估计生存率及其标准误、绘制生存曲线等。例如，根据乳腺癌患者手术后的生存资料，可以估计不同时间点的生存率及其标准误，如1年生存率、3年生存率、5年生存率等，还可以绘制生存曲线，观察乳腺癌患者手术后的生存过程。常用方法有乘积极限法和寿命表法。第20页/共81页2023/3/27生存分析与Cox回归214.生存分析的基本内容 比较生存过程：获得生存率及其标准误的估计值后，可进行两组或多组生存

9、曲线(生存过程)的比较。例如，比较两种不同治疗措施治疗恶性肿瘤患者的生存曲线，可了解哪种治疗措施较优，从而为临床决策提供依据。常用方法有对数秩检验。第21页/共81页2023/3/27生存分析与Cox回归224.生存分析的基本内容 生存过程的影响因素分析 例如，为了改善鼻咽癌患者的预后，应先了解可能影响患者预后的因素，如年龄、病程、病情、术前健康状况、有无淋巴结转移、术后有无感染、辅助治疗措施、营养状况等，通过随访收集患者术后的生存时间和上述因素的资料，然后采用多因素生存分析方法确定影响患者预后的主要因素，从而为在手术前后进行预防或干预提供参考依据。常用的多因素生存分析方法：Cox比例风险回归

10、模型第22页/共81页2023/3/27生存分析与Cox回归23生存分析的基本方法统计描述统计指标:生存率、中位生存时间列表和绘图:生存曲线统计推断参数法：指数分布、威布尔分布等非参数法：log-rank检验(单因素分析)半参数法：Cox回归(多因素分析)第23页/共81页2023/3/27生存分析与Cox回归24二、生存资料的统计描述生存率及其标准误中位生存期生存曲线第24页/共81页2023/3/27生存分析与Cox回归25基 本 概 念第25页/共81页2023/3/27生存分析与Cox回归26基 本 概 念第26页/共81页2023/3/27生存分析与Cox回归27基 本 概 念第27

11、页/共81页2023/3/27生存分析与Cox回归28基 本 概 念第28页/共81页2023/3/27生存分析与Cox回归29基 本 概 念第29页/共81页2023/3/27生存分析与Cox回归30基 本 概 念第30页/共81页2023/3/27生存分析与Cox回归31常 用 方 法乘积极限法(product-limit method)也称为-法(Kaplan-Meier法)小样本资料或大样本未分段(未按时间分组)资料寿命表法(life table method)大样本分段(按时间分组)资料第31页/共81页2023/3/27生存分析与Cox回归321.未分组资料(乘积极限法/Kapla

12、n-Meier法)第32页/共81页第33页/共81页第34页/共81页2023/3/27生存分析与Cox回归352.分组资料(寿命表法 life-table method)第35页/共81页第36页/共81页中位生存期为2.3第37页/共81页生存资料的统计描述SPSS实现例1第38页/共81页第39页/共81页第40页/共81页生存资料的统计描述SPSS实现例2第41页/共81页第42页/共81页第43页/共81页2023/3/27生存分析与Cox回归44三、生存曲线比较(单因素分析)第44页/共81页2023/3/27生存分析与Cox回归45 对上例进行分析时两种错误的做法：错误1：采用

13、平均生存时间而不是中位生存时间来表示生存时间的平均水平。错误2：采用常规t检验或方差分析进行组间比较。第45页/共81页第46页/共81页2023/3/27生存分析与Cox回归47三、生存曲线比较(单因素分析)第47页/共81页2023/3/27生存分析与Cox回归48三、生存曲线比较(单因素分析)第48页/共81页第49页/共81页第50页/共81页2023/3/27生存分析与Cox回归51 生存资料分析的基本要求样本应由随机抽样得到，要保证一定的样本含量死亡例数不宜太少截尾例数不宜太多生存时间应尽可能精确第51页/共81页2023/3/27生存分析与Cox回归52对数秩检验(log-ran

14、k)注意事项属于单因素分析方法，适用于两组及多组间的比较可用于时间未分组的资料，也可用于时间分组资料各组间生存时间的比较根据各组生存曲线的高低及中位生存时间判断需满足生存资料的基本要求，且各样本生存曲线不能交叉生存曲线若出现交叉，则提示可能存在混杂因素，应采用分层对数秩检验或Cox比例风险回归模型进行分析第52页/共81页生存曲线比较SPSS实现例3第53页/共81页第54页/共81页2023/3/27生存分析与Cox回归55四、Cox回归对生存时间资料进行单因素分析(log-rank检验)时，要求各对比组在非处理因素方面均衡可比，而实际工作中却很难做到；很多因素会对生存时间产生影响，需用多因

15、素分析方法。Cox回归也称比例风险模型(Proportional hazards model)是对生存资料进行多因素分析常用的统计方法第55页/共81页2023/3/27生存分析与Cox回归56四、Cox回归Cox回归模型的一般形式回归系数的解释Cox回归应用条件回归系数的估计及假设检验Cox回归分析应注意的问题第56页/共81页2023/3/27生存分析与Cox回归571.Cox回归模型的一般形式第57页/共81页第58页/共81页2023/3/27生存分析与Cox回归591.Cox回归模型的一般形式第59页/共81页2023/3/27生存分析与Cox回归602.回归系数的解释相对危险度 R

16、R:两个风险函数(率)之比 (风险比)当Xi为有无某危险因素时(0-1变量)第60页/共81页2023/3/27生存分析与Cox回归612.回归系数的解释在其他自变量固定不变时，自变量Xi每改变一个单位，得到的相对危险度RR的对数值。RR=exp()0,RR 1 危险因素0,RR 1 无作用 0,RR 1 保护因素RR的解释与自变量的编码有关第61页/共81页第62页/共81页第63页/共81页2023/3/27生存分析与Cox回归643.Cox回归应用条件独立性等比例风险风险比与时间无关，为常数自变量的取值及作用大小不随时间变化而变化 例：在研究10中，有糖尿病的发生心脏病的风险始终是无糖尿

17、病的3倍。对数线性：自变量与对数风险比呈线性关系 第64页/共81页第65页/共81页第66页/共81页第67页/共81页2023/3/27生存分析与Cox回归684.回归系数估计及假设检验回归系数的估计极大似然法(偏似然函数)假设检验方法似然比检验Wald卡方检验Score比分检验借助统计软件估计回归系数，进行假设检验对 SPSS软件：Analyze/Survival/Cox Regression第68页/共81页2023/3/27生存分析与Cox回归695.Cox回归分析应注意的问题等比例风险的假定影响因素(协变量)各水平的生存曲线无交叉协变量与时间交互项无统计学意义回归系数解释、变量赋值

18、与编码、假设检验、建模策略等与logistic回归类似样本含量估计注意时依协变量Cox模型第69页/共81页Cox回归分析SPSS实现例8.8第70页/共81页Cox回归分析SPSS实现例8.7第71页/共81页第72页/共81页2023/3/27生存分析与Cox回归73上海市青浦区2002 2008 年新发胃癌患者队列生存分析第73页/共81页第74页/共81页第75页/共81页牙齿缺失对食管癌的生存分析第76页/共81页第77页/共81页第78页/共81页第79页/共81页2023/3/27生存分析与Cox回归80小 结生存资料的特点生存资料统计描述方法非参数的假设检验：单因素分析log-rank检验多因素分析Cox回归第80页/共81页2023/3/27生存分析与Cox回归81感谢您的观看！第81页/共81页