第十九章 比例风险模型——Cox回归.ppt

第十八章生存分析和COX回归;第一节 基本概念; 生存期不同于一般指标的二个特点： 1.有截尾数据（censored data) 随访中未能知道病人的确切生存时间，只知道病人的生存时间大于某时间。 （1）病人失访或因其他原因而死亡---失访 （2）到了研究的终止期病人尚未死亡---终访 截尾数据可记为t+,如： 4+ = 生存时间大于4年。 虽然截尾数据提供的信息是不完全的，但不能删去，因为这不仅损失了资料，而且会造成偏性。;2. 生存期的资料一般不服从正态分布。 由于上述原因，常用的统计方法不适用，而要用特殊的统计方法。 生存分析是指对于生存期这一指标进行分析的一系列特殊的统计方法。 ; 生存时间不一定专用于死与活的情况,生存时间(存活时间)可定义为从某种起始事件到达某终点事件所经历的时间跨度。例如急性白血病病人从治疗开始到复发为止之间的缓解期;冠心病病人在两次发作之间的时间间隔;已作输卵管结扎的妇女从施行输卵管吻合手术后至受孕的时间间隔;在流行病学研究中,从开始接触危险因素到发病所经历的时间等都可作为生存时间用作生存分析。 有时还收集一些有关因素(称为自变量或协变量), 以分析这些协变量是否对生存时间有影响,影响的大小，是缩短或延长生存时间。这可以通过Cox回归进行分析,因此,Cox 回归可看成带有协变量的生存分析。; 包括： (1)开始观察日期，终止观察日期---生存时间 (2)结局（最终的观察到的是死亡还是存活） 死于该病---完全数据 存活或死于其他原因---截尾数据 每个生存期数据要用2个变量表示：观察到的生存时间和是否截尾（如：用1表示截尾，用0表示死亡；4+ 用4，1表示；4用4，0表示）。 (3)协变量---各种影响生存期长短的因素。;第二节 描述生存时间分布规律的函数;;§1.2 概率密度函数 (Probability Density Function) 简称为密度函数,记为f(t),其定义为: f(t)=lim (一个病人在区间(t,t+△t)内死亡概率/△t) 它表示死亡速率的大小。如以t为横坐,f(t) 为纵坐标作出的曲线称为密度曲线,由曲线上可看出不同时间的死亡速率及死亡高峰时间。纵坐标越大,其死亡速率越高,如曲线呈现单调下降,则死亡速率越来越小,如呈现峰值,则为死亡高峰。;§1.3 风险函数(Hazard Function) 用h(t)表示,其定义为: h(t)=lim(在时间t生存的病人死于区间(t,△t)的概率/△t) 由于计算h(t)时,用到了生存到时间t,这一条件,故上式极限式中分子部分是一个条件概率。可将h(t)称为生存到时间t的病人在时间t的瞬时死亡率或条件死亡速率或年龄别死亡速率。当用t作横坐标,h(t)为纵坐标所绘的曲线,如递增,则表示条件死亡速率随时间而增加,如平行于横轴,则表示没有随时间而加速(或减少)死亡的情况。; 风险函数的不同情况： 常数， 如：死于飞机失事。 下降， 如：急性损伤。 上升， 如：持续接触危险因素。 澡盆样，如：人的一生。 ; 生存分析目的： （1）估计生存函数。 （2）比较各组的生存函数。 （3）研究影响生存期长短的因素。;第三节 生存率的估计方法;;一. 乘积极限法(Product-Limit Method); 例22.1 用某中药加化疗(中药组)和化疗(对照组)两种疗法治疗白血病后, 随访记录各患者的生存时间,不带+号者表示已死亡,即完全数据,带+ 号者表示尚存活,即截尾数据,试作生存分析。时间单位为月。 中药组 10,2+,12+,13,18,6+,19+,26,9+,8+,6+,43+,9,4,31,24 对照组 2+,13,7+,11+,6,1,11,3,17,7 ;资料中药组积限法计算生存率 ───────────────────────────────────── 时间 状态 期初人数 死亡人数 条件生存率 累积生 di ∑di/ni(ni-di)累积生存 ti si ni di (ni-di)/ni 存率^S(ti)ni(ni-di) 率标准误 ① ② ③ ④ ⑤ ⑥ ⑦ ⑧ ⑨=⑥√⑧ ───────────────────────────────────── 2 活 4 死 15 1 0.9333 0.93