生存分析讲课.ppt

生存分析SurvivalAnalysis,第四军医大学西京医院神经内科刘学东,第一节概述第二节生存时间资料的特点第三节小样本生存率的Kaplan-Meier估计第四节大样本生存率的寿命表法估计第五节生存曲线比较的假设检验第六节Cox比例风险回归模型第七节实例分析,第一节概述,一基本概念生存分析：是将事件的结果（终点事件）和出现这一结果所经历的时间结合起来分析的一种统计分析方法。二特点含有截尾数据生存时间一般不呈正态分布,三历史最早可上溯至十九世纪的死亡寿命表现代的生存分析则开始于二十世纪三十年代工业科学中的相关应用二次世界大战：武器装备的可靠性研究，这一研究兴趣延续到战后。此时生存分析都集中在参数模型。,二十世纪六七十年代：医学研究中大量临床试验的出现，要求方法学有新的突破，生存分析研究开始转向非参数方法。D.R.Cox在72年提出的比例风险模型为此做出了划时代的贡献。现在，生存分析方法的在医学领域得到了广泛的应用，而通过医学研究要求的不断提高。,四主要研究内容描述生存过程比较生存时间分析影响生存时间的因素,五基本方法1非参数法：不考虑资料的分布形式，只根据样本提供的顺序统计量对生存率进行估计，常用有K-M法、Lifetable法2参数法：假定生存时间服从特定的参数分布，根据已知分布特点对影响生存时间的因素进行分析，常用有：指数分布法、对数正态回归分析法3半参数法：兼有参数法和非参数法的特点，分析影响生存时间和生存率的因素，常用有Cox比例风险回归模型,第二节生存时间资料的特点及描述,一、起始事件与失效事件起始事件：研究对象生存过程的起始特征事件失效事件(终点事件、死亡事件)：研究者关心的特定结局相对而言的，都由研究目的决定，失效事件并非一定死亡，死亡也并非一定发生了失效事件。分析的基石,二、生存时间（survivaltime）任何两个有联系事件之间的时间间隔，用t表示狭义：发病至死亡的时间跨度广义：起始事件到终点事件的时间间隔起始事件终点事件服药痊愈手术切除死亡出院复发,三、完全数据与删失数据完全数据：可获得准确的生存时间，所提供信息完全删失数据(截尾值)：随访中未能观察到病人的终点事件，右上角加“+”原因：失访死于其它疾病观察结束时病人尚存活,四、数据结构生存信息起始时间终点时间结局及原因失访存活死亡,预后信息一般信息个人史及既往史临床资料检验结果组织病理分子生物学影像学,与亚饭店,五、资料特点1含截尾数据2分布类型复杂六、资料的基本要求1.终点事件例数不宜太少2.截尾原因无偏性3.生存时间尽可能记录精确,七、统计描述指标1.死亡率、死亡概率、生存概率2.生存率及其标准误3.生存曲线4.中位生存时间及四分位数间距,1.死亡率、死亡概率、生存概率(1)死亡率(mortalityrate，deathrate)表示某单位时间内的死亡强度。年平均人口数=(年初人口数+年末人口数)2,(2)死亡概率(mortalityprobability)指死于某时段内的可能性大小。d/n0年内有删失，分母用校正人口数：校正人口数=年初人口数删失例数/2n0c/2,(3)生存概率(survivalprobability)指某单位时段开始时存活的个体到该时段结束时仍存活的可能性的大小。(n0-d)/n0分子为年末尚存人数，若年内有删失，分母用校正人口数。,60例肺癌病人，术后每年死亡10例，无删失。N=60,术后年内年内年初年死亡死亡生存(t+1)年年死亡截尾观察平均率概率概率生存率数例数例数例数例数p=1-qtdcn0n=no-d/2m=d/nq=d/n0(n0-d)/n0(n0-d)/N(1)(2)(3)(4)(5)(6)(7)(8)(9)01006055.185.167.833.83311005045.222.200.800.66721004035.286.250.750.50031003025.400.333.667.33341002015.667.500.500.167注意：死亡率与死亡概率的分子相同，但分母不同；生存概率与生存率的分子相同，但分母不同,2.生存率生存率(survivalrate)(累积生存概率cumulativeprobabilityofsurvival)指病人经历t个单位时间后仍存活的概率。生存概率指单个时段的概率，生存率指从0t多个时段的积累概率。反映了病人的疾病严重程度、进展快慢或凶险程度等。,公式1用于完全数据公式2用于删失数据以及完全数据S(t)=P(Ttk)=p1p2pk,3生存曲线（survivalcurve）以观察（随访）时间为横轴，以生存率为纵轴，将各个时间点所对应的生存率连接在一起的曲线图。描述生存过程比较生存过程中位生存期,(月）,4中位生存时间（mediansurvivaltime)即生存时间中位数，表示50%的个体可存活的时间(集中趋势)计算：图解法线性内插法：找出两个生存率S(ti-1)和S(ti)，使S(ti-1)0.5,S(ti)8)=0.565,t=10,P(X10)=0.478(8-10):(8-t)=(0.565-0.478):(0.565-0.5),t=9.494(月)生存率的计算方法,第三节小样本生存率的Kaplan-Meier估计,Kaplan-Meier法(K-M法、乘积极限法)Kaplan-Meier法由Kaplan和Meier于1958年提出，直接用概率乘法定理估计生存率，故称乘积极限法（product-limitmethod），是一种非参数法，适用于小样本和大样本。,基本思想：将生存时间由小到大依次排列，在每个死亡点上，计算其期初人数、死亡人数、死亡概率、生存概率和生存率。本法只估计死亡时点的生存率,(a)研究终止在475天(b)研究终止在474天图12-3乘积极限法生存曲线(阶梯形)及其半数生存期（Md=158天）,第四节大样本生存率的寿命表法估计,寿命表法(LifeTables)基本原理:是先求出病人在治疗后活过n年后再活过下一年的可能性，统计学上称为生存概率，根据概率论的乘法定律将逐年生存概率相乘即可得出一定年限的生存率。适用于大样本资料,K-M法和Lifetable法比较,第五节生存曲线比较的假设检验,log-rank检验（对数秩检验、时序检验）基本原理：对每组生存数据根据其在各个时刻中尚存活的对象数计算期望死亡数，然后将期望死亡总数与实际死亡数进行比较将生存曲线作为整体对曲线与曲线间的差异进行比较，偏重于远期疗效评价，属非参数检验,比较两组或多组生存曲线，实为一单因素分析要求各组生存曲线不能交叉，如交叉提示存在混杂因素，应采用分层分析方法或多因素方法来校正当假设检验有统计意义时，可从以下几方面来评价各组效应大小：生存曲线图目测判断、中位生存期比较,第六节Cox比例风险回归模型,一概述简称Cox回归模型，是一种允许资料有“删失（或截尾）”数据存在的，可以同时分析众多因素对生存时间影响的多变量生存分析方法，是一种半参数方法。,该模型由英国统计学家D.R.Cox于1972年提出，主要用于肿瘤和其它慢性病的预后分析，也可用于队列研究的病因探索。优点：多因素分析方法不考虑生存时间分布利用截尾数据,二基本形式h(t,X)t时刻风险函数、风险率或瞬时死亡率h0(t)基准风险函数，即所有协变量都取0时t时刻风险函数。X1、X2、Xp协变量、影响因素、预后因素。1、2、p回归系数。,p：在其它协变量不变的情况下，协变量Xp每改变一个测定单位所引起相当危险度的自然对数的改变量RR:相当危险度，两个分别具有协变量Xi与Xj的病人其风险函数（危险度）之比RR=expj(xi-xj),如协变量treat的=-1.617，RRtreat=0.199,表示treat变量水平2与1比较，treat=2的危险度是treat=1的0.199倍，提示治疗方案2优于治疗方案1。而age的=0.119，RRage=1.127，表明年龄每增加一岁，死亡的可能性增加1.127倍。,0，RR1，说明变量X增加时，危险率增加，即X是危险因素。0，RR1，说明变量X增加时，危险率下降，即X是保护因素。=0，RR=1，说明变量X增加时，危险率不变，即X是危险无关因素。,三参数估计与假设检验参数估计最大似然法假设检验(模型中变量的剔除和引入)似然比检验得分检验Wald检验,四因素筛选与最优模型的建立因素初步筛选2检验、log-rank检验等单变量Cox模型分析（P=0.10）直接进行逐步Cox模型分析,最佳模型建立1、Enter2、Forward:Conditional3、Forward:LR4、Forward:Wald5、Backward:Conditional6、Backward:LR7、Backward:Wald,五检验水准检验水准包括引入的检验水准和剔除的检验水准。一般地，剔除引入初步的、探索性的研究，或变量数较少时，可取0.10或0.15。设计严谨的、证实性的研究，或变量数较多时，可取0.05。,六统计描述1回归系数和标准回归系数标准回归系数：观察值经过标准化变换后所求得。相对值，比较自变量对模型贡献，反映因素对生存时间影响强度2相对危险度及可信区间无序分类变量，转换亚变量分析,3个体预后指数（PI）PI=预后指数越小，预后越好；预后指数越大，预后越差。据PI的大小可将病人分为不同的危险组，绘制分组的生存曲线，直观比较各危险组的生存情况。,例如：传统治疗方式，淋巴结转移，PI2.6929传统治疗方式，淋巴结无转移，PI1.7616新治疗方式，淋巴结转移，PI0.9313新治疗方式，淋巴结无转移，PI0,七分析结果（结果解释）与生存相关的因素因素作用大小及方向：保护因素还是危险因素、相对危险度的大小。因素作用大小排序：标准化回归系数的绝对值。,八结果报告1起始时间与终止时间2结局事件3删失数据的种类与原因4生存率估计值、置信区间及其统计学方法,5生存率比较的统计学方法(如logrank法)及其p值6给出Cox回归模型，呈现解释变量与风险之间的联系风险比（hazardratio）及其置信区间,九应用范围因素分析具有logistic回归模型的所有优点同时考虑生存结局和生存时间，可处理截尾数据不考虑基准风险的分布,生存预测强调设计的重要性专业知识角度：选择疾病种类、终点事件、影响因素及结果的专业解释等。统计学角度：样本例数、因素的赋值、生存时间准确、因素筛选方法、结果的统计学解释等。,十注意事项应用条件除生存资料的基本要求外，还要求因素对生存时间的作用不随时间变化（比例风险假定）。设计阶段资料的代表性和可靠性将尽可能的因素都包括在调查分析中,筛选变量容易收集费用明确规定生存时间样本量为观察协变量的520倍模型拟合多元共线性：自变量间不独立，相互间有一定线性依存关系,选入模型的变量是统计学上的有关变量，不一定都与生存时间有因果关系。其中某些变量可能只是伴随关系而已，未选入模型的变量不一定全是无关变量，要考虑是否模型内的某些变量代替了它的作用，或因例数不够，或试验中对该变量进行了控制而引起。,模型拟合优度考察：据预后指数PI分组，比较各组基于Cox模型的生存曲线与基于kaplan-Meier法估计的生存曲线，如两组曲线吻合较好，表明Cox模型拟合较好。生存率分析：生存曲线不能随意延长，也不能轻易用于预测预报，经过大量研究所得的生存曲线才有可能推广应用。,十一常见问题1生存时间的计算和选择2生存率计算方法使用不当3忽略结局，采用t检验、线性回归分析生存时间4忽略生存时间，采用Logistic回归分析死亡率,5采用平均生存时间表示生存时间的平均水平6采用常规t检验或方差分析进行组间生存时间比较7丢弃截尾数据，只考虑确切数据8将截尾数据当作确切数据处理,例某研究者追踪100名癌症患者，经治疗后的生存情况，随访第1年有30人死亡，随访第2年时，有20人死亡，有40人无法联系失访。方法一：去掉截尾数据方法二：不去掉截尾数据总例数N=60N=100生存数生存率生存数生存率1年3030/60=50%7070/100=70%,第七节实例分析,1研究设计目的专业：描述缺血性脑卒中患者病后生存状况，并分析影响死亡的危险因素统计学：估计：如根据缺血性脑卒中患者治疗后的生存时间资料，估计不同时间的生存率、生存曲线以及中位生存期等。,比较：对不同因素对生存率的影响进行比较影响因素分析：目的是为了探索和了解影响生存时间的因素，如病人的年龄、性别、病程、严重程度等。预测：具有不同因素水平的个体生存预测,方法：回顾性研究确定起始事件、终点事件、生存时间起始事件：患者发病终点事件：因缺血性脑卒中死亡生存时间：患者发病至死亡、失访