多元统计学课件.pptVIP

Cox’s proportional hazard regression 比例风险模型——Cox回归 是一种允许资料有“截尾或终检”数据存在的，可以同时分析众多因素对生存时间影响的多变量生存分析方法。 Cox回归 第一节 生存时间 一．生存时间概念： 狭义的角度来讲，生存时间是患某种疾病的病人发病到死亡所经历的时间跨度。 广义的角度来讲，可以把生存时间定义为从某种起始事件到达某种终点事件所经历的跨度。 例如： 我们可以把下列事件作为起始事件和终点事件： 起始事件 终点事件 疾病确诊 死亡 治疗开始 治愈 症状缓解 疾病恶化 接触毒物 出现毒性反应 接触危险因素 发病 二．两类生存时间数据：（一）完全数据： 在随访工作中，当观察到了某患者的明确结局时，该病人所提供的关于生存时间的信息是完整的。 我们把到达了明确结局的病人的生存时间数据称为完全数据。用符号“t”表示。 （二）截尾数据  1．定义：在随访工作中，由于某种原因未能观察到病人的明确结局，不知道病人的确切生存时间，随访资料不完整。这样的资料称为截尾数据。 2．价值：提供了观察期内的信息，生存时间不会短于观察时间。用符号“t+”表示。 3．产生截尾现象的原因： （1）病人失访 （2）病人的生存期超过了研究的终止期。 （3）在动物实验中，有时预先规定观察期限。 （三）图示： ? ? ? ? ? ? ? ? 三．生存时间资料的整理： 设总共观察了n例病人的生存时间，记第i名病人的生存时间为ti，则全部n例病人的生存时间可以记为 t1,t2,…,tn 如图1.1中的生存时间（年）按病人号的顺序排列为 2,4+,3,1+ 如果将这些生存时间按由小到大的顺序排列，则得到一个有序的生存时间列为t(1)≤t(2)≤…≤t(n) 图1.1中的生存时间有序系列为：1+,2,3,4+。 几种不同类型的生存时间的例子： 例1.1 20名行输卵管结扎术的妇女经峡部-峡部吻合术后的受孕时间（月）为： 1,1,2,3,3,4,4,4,6,6,8,9,9,10,11,12,13,15,17,18。 此例中的生存时间均为完全数据，并已按由小到大的顺序排列整理。 例1.2 23名行输卵管结扎术的妇女经壶腹部-壶腹部吻合术后的受孕时间（月）为： 1,3,5,5,5,6,6,6,7,8,10,10,14+,17,19+,20+,22+,26+,31+,34,34+,44,59。 此例为有序生存时间资料，并包括有截尾数据。 可在随访期内任何时点上发生截尾类型称为随机截尾（random censoring）。本例中的截尾就是随机截尾。 例1.3 15只雌性大白鼠接触毒物DHG后观察12周，其生存时间为：4,6,8,9,9,10,10,10,11,12,12,12+,12+,12+,12+。 此例是一个定时截尾数据的资料，所有到第12周未死亡的动物的生存时间都属于截尾数据。 例1.4 当观察单位较多时，可以按一定的时间区间分段整理。设在第i个时间内（ti,ti+1）开始时的病人数为ni，在此区间内的死亡数为di，则在第i+1个区间开始时的病人数为ni+1=ni-di。分组资料情况见下表。 四．生存分析的方法： 一般可以分为参数、非参数、半参数三类。 1、生存时间的分布符合某一特定类型，如对数正态分布、weibull分布、指数分布等，则可以用特定的分布函数分析，这称之为参数法。 2、若不知道生存时间的分布类型，而对分布或其某些特征作推断，就只能用半参数或非参数法。 如：用寿命表估计期望寿命；用Kaplan-meier法求生存率，作生存曲线；用logrank检验等作不同组别生存过程差异的显著检验，均属非参数法。 3、由于Cox 比例风险模型不是直接利用时间t的全部信息，而只是利用了时间t所提供的顺序统计量的信息，故有人称之为半参数方法。 第二节 Cox比例风险回归 一．模型结构： 设有n名病人，第i名病人的生存时间为ti，同时该病人具有一组伴随变量xi1,xi2,xi3,…,xip。该病人生存到时间ti的风险函数hi(t)是其基础风险函数h0(t)与相应伴随变量的函数的乘积，写成数学表达式为 ? hi(t)=h0(t)·exp(β1xi1+…+βpxip)? ? Cox证实伴随变量的函数具有指数形式，故Cox的比例风险回归可写为 hi(t)=h0