4生存分析.ppt

生存分析 主要内容回顾 生存分析中的基本概念生存率的估计和生存曲线 统计描述 生存曲线的log rank检验 单因素分析 Cox比例风险回归模型 多因素分析 疾病确诊死亡治疗开始痊愈手术复发接触毒物出现毒性反映吸毒HIV感染HIV感染自杀未遂出生发病 起始事件和终点事件 生存时间 survivaltime 又称失效时间 failuretime 定义为观察对象从起始事件到终点事件所经历的时间 用t表示 生存时间的分布不是正态分布 多为正偏峰分布 图102名黑色素瘤患者的生存时间分布 time Month 0 5 10 15 20 25 30 35 40 45 0 10 20 30 40 生存时间的分布为非正态分布对数正态分布指数分布 Weibull分布Gamma分布 logistic分布 含有删失数据 censoreddata 删失数据提供的信息是不完全的 uncompleteddata 数据特点 数据结构 7名患者手术到死亡的时间39 152020 2630 事件结局变量 T 生存时间 Y 删失指示量 Y 1表示未删失 Y 0表示删失 主要内容 生存分析中的基本概念生存率的估计和生存曲线生存曲线的log rank检验Cox比例风险回归模型 生存率的估计 小样本时 不需要对生存时间进行分组 常用乘积极限法 product limited PL 也称Kaplan Meier K M 法估计生存率 大样本时 需要对生存时间进行分组 常用寿命表 lifetimetable 法估计生存率 属于非参数方法 图29 5儿童横纹肌肉瘤加化疗组 treat 与对照组 control 缓解曲线比较 生存曲线的比较 log rank检验 Dataaa inputtygroup cards 2011212301160290119021001191210012012121120121501241215112412160130121811311224113412300142123611441240115312451159129026212 Proclifetestmethod plplots s timet y 1 stratagroup run 绘制生存曲线 PL法估计生存曲线 单因素分析方法Log rank检验 多因素分析方法Cox回归分析 Cox回归分析 比例危险率回归 英国统计学家D R Cox于1972年提出的比例危险率回归模型 proportionalhazardsregressionmodel 简称Cox回归 主要用于肿瘤与其他慢性病的预后分析 也可用于队列研究的病因探索 模型结构h t x 具有协变量X的个体在t时刻的风险函数 表示这些个体在t时刻的瞬时危险率或死亡率 h0 t 称为基准风险函数 是与时间有关的任意函数 表示所有X取值为0时的个体在t时刻的瞬时危险率或瞬时死亡率 X 与生存时间有关的协变量或交互项 可以定性或定量在整个观察时间不随时间的变化而变化 半参数模型 1 建立偏似然函数用最大似然估计方法估计系数2 回归系数的检验方法似然比检验Score检验Wald检验 回归系数的估计与检验 其他协变量不变时 变量每变化一个单位所引起的相对危险度的自然对数改变量 j ln RRj 危险无关因素 危险因素 保护因素 回归系数 j的含义 Cox回归分析的应用实例 调用Cox模型分析模块 分析数据集为cox 建立生存时间为t 截尾指示变量为d的Cox模型 d取值为1时表示截尾 选入和剔出水平均为0 05 计算筛选因素的相对危险度及其95 的可信区间 治疗方式 x4 相对危险度为5 820 说明传统的治疗方式和新的治疗方式相比 病人死亡的风险为5 820倍 相对危险度的95 的可信区间为1 989 17 031 淋巴结是否转移 x5 的RR值的含义与治疗方式相同 例3某医药公司与中山医科大学协作 研究天花粉注射治疗绒癌的疗效问题 将16只体表接种绒癌成功的裸鼠随机分为4组 其中一组作空白对照 其余三组分别注射天花粉 甲药和乙药 经整理 数据如表所示 试作Cox回归分析 dataaa inputtdv0tr1tr2tr3vitCdday cards 1925000118171600011919370001181619000118142510010181318100111716251001114910100111592201011151025010111114250101113123701011121737001119142900111121313001011217310011110 procphreg modelday d 0 tdv0tr1tr2tr3 selection forward run h t x h0 t exp 0 42td 2 94tr1 1 其他因素保持不变的情况下 带瘤天数具有统计学意义 x2 6 65 P 0 05 带瘤天数每增加一个单位 死亡风险增加0 522倍 2 其他因素保持不变的情况下 注射天花粉具有统计学意义 x2 7 53 P 0 05 注射天花粉者死亡风险是不注射者的0 053倍 例 某医院泌尿外科医师选择1996 2000年间经手术治疗的膀胱肿瘤患者 对