生存分析PPT学习教案

1、会计学1 生存分析生存分析 疾病确诊 死亡 治疗开始 痊愈 手 术 复发 接触毒物 出现毒性反映 吸 毒 HIV感染 HIV感染 自杀未遂 出 生 发病 第1页/共27页 第2页/共27页 图 102名黑色素瘤患者的生存时间分布 time (Month) 051015202530354045 0 10 20 30 40 第3页/共27页 生存时间的分布为非正态分布 对数正态分布 指数分布， Weibull分布 Gamma分布，logistic分布， 含有删失数据(censored data) 删失数据提供的信息是不完全的 (uncompleted data) 数据特点 第4页/共27页 T39

2、1520202630 Y1011010 T：生存时间；：生存时间； Y：删失指示量，：删失指示量，Y1表示未删失；表示未删失；Y0表示删失。表示删失。 第5页/共27页 主要内容 生存分析中的基本概念生存分析中的基本概念 生存率的估计和生存曲线生存率的估计和生存曲线 生存曲线的生存曲线的log-rank检验检验 Cox比例风险回归模型比例风险回归模型 第6页/共27页 生存率的估计 小样本时，不需要对生存时间进行分组，小样本时，不需要对生存时间进行分组， 常用常用 乘积极限法乘积极限法(product-limited，PL)，也称也称 Kaplan-Meier(K-M)法法估计生存率。估计生存

3、率。 大样本时，需要对生存时间进行分组，常大样本时，需要对生存时间进行分组，常 用寿命表用寿命表(lifetime table)法估计生存率。法估计生存率。 属于非参数方法属于非参数方法 第7页/共27页 第8页/共27页 01224364860 t (month) 0.0 0.2 0.4 0.6 0.8 1.0 S(t) treat control 图29.5儿童横纹肌肉瘤加化疗组(treat)与对照组(control)缓解曲线比较 第9页/共27页 1 )( 2 2 2 v v ea j jj L 第10页/共27页 绘制生存曲 线 PL法估计生存曲 线 第11页/共27页 第12页/共2

4、7页 单因素分析方法单因素分析方法 Log-rank 检验 多因素分析方法多因素分析方法 Cox回归分析回归分析 第13页/共27页 Cox回归分析(比例危险率回归) 英国统计学家英国统计学家D.R.Cox于于1972年提出的年提出的比例比例 危险率回归模型危险率回归模型(proportional hazards regression model)，简称，简称Cox回归。回归。 主要用于肿瘤与其他慢性病的预后分析，也可用主要用于肿瘤与其他慢性病的预后分析，也可用 于队列研究的病因探索于队列研究的病因探索 第14页/共27页 01122mm h(t,x)h (t)exp(xxx ) 半参数模型？

5、 第15页/共27页 1、建立、建立偏似然函数偏似然函数用用最大似然估计最大似然估计方法估计系数方法估计系数 2、回归系数的检验方法、回归系数的检验方法 似然比检验似然比检验 Score 检验检验 Wald 检验检验 回归系数的估计与检验 第16页/共27页 少表示病人死亡的风险减取值越大，时， 无影响的取值对时， 大表示病人死亡的风险越取值越大，时， ,),(0 ),(0 ,),(0 XthX XthX XthX ii ii ii ，危险无关因素 危险因素危险因素 保护因素保护因素 exp() j RRb 回归系数 j的含义 第17页/共27页 第18页/共27页 第19页/共27页 调用调

6、用Cox模型分析模块模型分析模块 ，分析数据集为，分析数据集为cox 建立生存时间为建立生存时间为t ，截尾指示变，截尾指示变 量为量为d 的的Cox模型，模型，d 取值为取值为1时时 表示截尾表示截尾 选入和剔出水平均为选入和剔出水平均为 0.05 计算筛选因素的相对危险计算筛选因素的相对危险 度及其度及其95%的可信区间的可信区间 第20页/共27页 Summary of Stepwise Selection Variable Number Score Wald Step Entered Removed In Chi-Square Chi-Square Pr ChiSq 1 x4 1 13

7、.0399 . 0.0003 2 x5 2 4.7039 . 0.0301 筛选过程总结筛选过程总结 Analysis of Maximum Likelihood Estimat Parameter Standard Hazard 95%HR Variable DF Estimate Error Chi-Square PrChiSq Ratio Confidence Limits x4 1 1.76128 0.54785 10.3356 0.0013 5.820 1.989 17.031 x5 1 0.93133 0.44455 4.3890 0.0362 2.538 1.062 6.066

8、最大似然估计结果最大似然估计结果 治疗方式（治疗方式（x4）：相对危险度为：相对危险度为5.820，说明传统的治疗方式和新的治疗方式相，说明传统的治疗方式和新的治疗方式相 比，病人死亡的风险为比，病人死亡的风险为5.820倍，相对危险度的倍，相对危险度的95%的可信区间为的可信区间为1.98917.031； 淋巴结是否转移（淋巴结是否转移（x5）的）的RR值的含义与治疗方式相同。值的含义与治疗方式相同。 第21页/共27页 例例3 3 某医药公司与中山医科大学协作，研究天花粉注射治某医药公司与中山医科大学协作，研究天花粉注射治 疗绒癌的疗效问题。将疗绒癌的疗效问题。将1616只体表接种绒癌成功

9、的裸鼠随机分只体表接种绒癌成功的裸鼠随机分 为为4 4组，其中一组作空白对照，其余三组分别注射天花粉、甲组，其中一组作空白对照，其余三组分别注射天花粉、甲 药和乙药，经整理，数据如表所示。试作药和乙药，经整理，数据如表所示。试作CoxCox回归分析。回归分析。 第22页/共27页 proc phreg; model day*d(0)=td v0 tr1tr2 tr3/selection=forward; run; 第23页/共27页 h(t,x)=h0(t)exp(0.42td-2.94tr1) 1、其他因素保持不变的情况下，带瘤天数具有统计学意义（ x2=6.65，P0.05 ），带瘤天数每增加一个单位，死亡风险增 加0.522倍。 2、其他因素保持不变的情况下，注射天花粉具有统计学意义 （x2=7.53，P0.05 ），注射天花粉者死亡风险是不注射者的 0.053倍。 第24页/共27页 例：某医院泌尿外科医师选择例：某医院泌尿外科医师选择1996-2000年间经手术治疗的年间经手术治疗的 膀胱肿瘤患者，对可能影响膀胱肿瘤术后生存的因素进行了膀胱肿瘤患者，对可能影响膀胱肿瘤术后生存的因素进行了 调查，研究者不仅要看是否出现了感兴趣的终点事件或结局调查，研究者不仅要看是否出现了感兴趣的终点事件或结局 （死亡）