生存分析概述及实例分析学习教案.pptx

《生存分析概述及实例分析学习教案.pptx》由会员分享，可在线阅读，更多相关《生存分析概述及实例分析学习教案.pptx（37页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、会计学1生存生存(shngcn)分析概述及实例分析分析概述及实例分析第一页，共37页。目录(ml)1概述2常用术语3分析方法4案例分析第1页/共37页第二页，共37页。概述(i sh)第2页/共37页第三页，共37页。定义(dngy)生存分析是研究生存现象和响应时间数据及其统计规律的一门学科。由于最初研究的关键事件是死亡，故称为生存分析。生存分析是统计科学的重要分支，其研究的两个重要变元为“事件”和“寿命(shumng)”。事件：生存分析中定义的事件有死亡、损坏、失败、解雇、病发等等。例如病人的死亡，产品的失效，疾病的发生，职员被解雇。寿命(shumng)：从记录开始到事件发生的时间。第3页/

2、共37页第四页，共37页。特点(tdin)生存分析的优点在于(ziy)其能够处理删失数据。生存分析的统计资料以生存时间为反应变量，此类资料的生存时间变量大多不服从正态分布，且由于删失值的存在，不适合用传统的分析方法处理。此时就应选用生存分析的方法。第4页/共37页第五页，共37页。研究(ynji)内容生存分析(fnx)研究的内容主要有以下两个方面：一 对生存过程的描述 二 分析(fnx)生存过程的影响因素并对生存的结局加以预测第5页/共37页第六页，共37页。应用领域生存分析虽然源自医学领域，但其在生物学，保险学，可靠性工程学，经济学，教育学，社会学等领域都有广泛的应用。比如：医疗科学中病人的

3、去世 保险行业中的赔偿 可靠性工程中产品的失效 金融领域中银行账户从开立到取消的时间的研究 教育行业中学生的中途(zhngt)退学 客户关系管理中的客户流失第6页/共37页第七页，共37页。常用(chn yn)术语第7页/共37页第八页，共37页。生存(shngcn)数据生存分析中所分析的数据通常称为分析数据，一般度量的是某个事件发生所经历的时间长度。生存数据可以分为完全(wnqun)数据和删失数据。完全(wnqun)数据：指提供了完整信息的数据。例如，在研究产品的失效时间时，某个样品从进入研究直到失效都在我们的观察中，可以得到该样品的具体失效时间，这就是一个完全(wnqun)数据。删失数据：

4、是指在观测期内，我们并没有看见个体的状态发生改变，无法确定个体具体的生存时间。又分为左删失数据，右删失数据，区间删失数据。第8页/共37页第九页，共37页。ABCDE删失数据(shj)示意图起始 观测时间(shjin)区间 终止死亡(swng)死亡死亡未知退出完全数据：A,观测期内死亡右删失数据：B，观测未终止时因故退出C，观测终止时尚未死亡左删失数据：D，死亡时间在某一时刻之前，具体时间未知区间删失数据：E，死亡时间位于某一区间，具体时间未知第9页/共37页第十页，共37页。生存(shngcn)函数生存函数（survival function），又称为累积生存率，我们用符号T表示个体的生存时

5、间（从开始记录到事件发生的时间），用 t 表示观测时间，将生存函数记作 S（t），是指个体生存时间大于 t 的概率。S（t）=P（T t），显然(xinrn)S（t）是非升函数，且S（0）=1,S(）=0，第10页/共37页第十一页，共37页。风险(fngxin)函数风险(fngxin)函数（hazard function），又称为瞬时死亡率，记作 h（t）。是指在t时刻存活的个体，在t+t 时刻死亡的概率。显然，h（t）非负，且无上限。第11页/共37页第十二页，共37页。分析方法第12页/共37页第十三页，共37页。按照是否使用(shyng)参数，可以将生存分析中的分析方法分为三类：参数方

6、法：若已经证明某事件的发展可以用某个参数模型很好地拟合，就可以用参数方法做该事件的生存分析。常用的参数模型有指数分布模型、对数分布模型、正态分布模型，威泊分布模型等。非参数方法：当被研究事件不能被参数模型很好地拟合时，可以采用非参数方法研究它的生存特征。常用的非参数方法包括生命表分析和K-M分析。半参数方法：它比参数模型灵活，与非参数方法相比更容易对分析结果进行解释。生存分析中使用(shyng)的半参数模型是Cox比例风险模型。第13页/共37页第十四页，共37页。非参数(cnsh)方法l生命表分析(fnx)lK-M分析(fnx)第14页/共37页第十五页，共37页。生命表分析(fnx)生命表

7、分析将观测时间分成时间段，按时间段逐个统计事件发生的情况，以此估计生存函数。假设共有k个时间段 t 0,t 1),t 1,t 2),t k-1,t k),每个区间中事件发生的次数(csh)分别为 d 1，d 2,d k,每个区间中的个体总数分别为 n 1,n 2,n k,所以在第 i 个区间个体存活的概率为（n i-d i）/n i ，而个体可以从第一个区间存活到第 i 个区间的概率(累积生存率）为：i=1,2，k，且S（t）为递减函数。第15页/共37页第十六页，共37页。K-M分析(fnx)Kaplan-Meier分析，也称为乘积(chngj)极限分析，是Kaplan和Meier在1958

8、年提出的一种估计生存函数的非参数方法。与生命表分析不同，K-M分析以事件发生的时间点将观测区间分段，用来估计生存函数。下举例说明其具体的分析过程。第16页/共37页第十七页，共37页。个体编号个体编号生存状态生存状态存活时间存活时间/小小时时1F312S653F1504S2205F300 下表记录了5个实验对象的存活(cn hu)时间，其中F代表失效，S代表存活(cn hu)，2和4为右删失数据。用S（t）表示实验对象(duxing)的累积存活概率，分时间段计算如右：t0,31):此区间5个实验对象均存活，故 S（t）=5/5=1.31,65):个体(gt)1在31小时死亡，故本区间 S（t）

9、=14/5=0.865,150):个体(gt)2在65小时退出实验，本区间无个体(gt)死亡，S（t）=0.84/4=0.8.150,220):个体(gt)3在150小时死亡，S（t）=0.82/3=0.53.220,300):个体(gt)4在220小时退出实验，本区间无个体(gt)死亡，S（t）=0.532/2=0.53.t=300时，个体(gt)5死亡，S（300）=0第17页/共37页第十八页，共37页。以SPSS对上例进行K-M分析(fnx)，结果如下：1.输入数据2.进行(jnxng)K-M分析第18页/共37页第十九页，共37页。参数设置第19页/共37页第二十页，共37页。输出(

10、shch)结果K-M分析(fnx)生存函数图第20页/共37页第二十一页，共37页。生命表分析(fnx)与K-M分析(fnx)的比较l生命表分析适用于大样本的情况，特别是没有个体数据的情形，主要优点(yudin)是对生存时间的分布没有要求。lK-M分析中时间区间的划分是以事件的发生为依据的，因此必须知道每个个体的生存时间数据，适用于小样本的情况。第21页/共37页第二十二页，共37页。半参数(cnsh)方法生存(shngcn)分析中我们常常遇到个体的生存(shngcn)状况受到多种因素影响的情况。这些对生存(shngcn)时间有影响的变量称为协变量。在分析生存(shngcn)数据时要将协变量的

11、影响考虑进去。Cox半参数模型就很好地解决了这个问题。它假定风险函数由两部分构成：基准风险函数和协变量线性组合的指数。Cox半参数模型又分为独立协变量比例风险模型和时间相依性协变量比例风险模型两种。二者的区别在于协变量的取值是否和时间有关。第22页/共37页第二十三页，共37页。Cox独立协变量比例(bl)风险模型该模型可写成如下形式：式中，Z1，Z2，Zm为协变量，这里的协变量与时间无关，1，2，m为对应协变量的未知参数。h 0（t）是基准风险函数。实际应用中常比较两个(lin)不同个体风险函数的比率，即危险率。可以证明危险率为常数，因此该模型又称为比例风险模型。当协变量与时间有关时，危险率

12、将不再是常数，此时称为时间相依性比例风险模型。第23页/共37页第二十四页，共37页。案例(n l)分析第24页/共37页第二十五页，共37页。l 两组小白鼠用来检验(jinyn)癌症的治疗状况。一组使用传统治疗方法，另一组使用试验方法，试验人员记录了小白鼠的存活时间及状态：Days为存活时间或观测时间；Status表示生存状态，取值1表示死亡，0表示存活；Group表示治疗方法，取0表示传统疗法，取1表示试验疗法，共有64组数据。第25页/共37页第二十六页，共37页。原始数据如下(rxi)：第26页/共37页第二十七页，共37页。首先用生命表分析方法对数据(shj)进行处理：1.输入(sh

13、r)数据2.选择(xunz)生命表分析第27页/共37页第二十八页，共37页。3.设置(shzh)参数第28页/共37页第二十九页，共37页。4.输出(shch)结果可以看出，大约在200天时(tinsh)两种治疗方法的生存函数相交，在200天以前传统治疗方法的存活率较高，而在200天以后试验方法的治疗效果明显优于传统治疗方法。中位数生存时间是生存率为50%时，生存时间的平均水平。从中位数生存时间来看，传统治疗方法的中位数为241天，试验方法的中位数为266天，明显高于传统治疗方法。可以(ky)判断试验方法的疗效相比传统治疗方法有所提高。第29页/共37页第三十页，共37页。用K-M方法对数据

14、进行(jnxng)处理,结果如下：生存函数(hnsh)分布和生命表分析的结果相似。K-M方法可以记录删失数据，且由于分段较多整体呈现密集的锯齿，而生命表分析的分布则较为平缓。在结果检验上，不同检验方法结果有所差异，其中Log Rank检验的p值小于0.05，表明两种治疗方法有显著性差异。第30页/共37页第三十一页，共37页。除了治疗方法对小白鼠的生存状况有影响，其他因素如性别，年龄，体重等都可能对其生存时间造成影响。加入这些数据后，用Cox独立(dl)协变量比例风险模型重新分析。1.输入(shr)数据2.设置(shzh)参数第31页/共37页第三十二页，共37页。3.输出(shch)结果 分

15、类变量是指不连续的变量，此例中治疗方法的值只取0和1，性别只取F和M。不同于体重、年龄这些连续变量，分类变量在计算风险比例时以参考类别作为参照(cnzho)。如在本案例中治疗方法这一因子以试验方法作为参照(cnzho)。计算结果为传统方法的风险率相对于参考的倍数。第32页/共37页第三十三页，共37页。第33页/共37页第三十四页，共37页。上表为模型系数的综合检验结果。可以看到p=0，小于0.05，说明这些因素(yn s)中有些变量对白鼠的生存时间有显著影响。第34页/共37页第三十五页，共37页。此表给出了各个变量的单个模型系数检验结果，可以看到体重变量的p=0，说明体重对风险函数有极显著影响。体重每增加1（盎司），风险大约为原来的1/3；治疗方法的p=0.068，影响几乎显著，传统方法的风险为试验方法的1.75倍；而性别和年龄(ninlng)对风险的影响微弱。结论：综合几种不同的分析方法，我认为综合以上三种模型的分析结果，我认为在本实验中白鼠体重是影响其生存(shngcn)时间的重要因素，治疗方法的不同也有比较重要的影响，而年龄和性别几乎不对其生存(shngcn)时间产生影响。第35页/共37页第三十六页，共37页。THANKSn n谢谢(xi xie)观看第36页/共37页第三十七页，共37页。