培训课件--生存分析20160518.ppt

SAS程序 Proc phreg 选择项； Model response*censor(删失值) =变量/选择项； Run; Model语句选择项 Selection=method backward forward stepwise Details 详细输出每一步过程 Include=n 模型中必须包括前n个变量 Slentry=p1 slstay=p2 p1p2 Alpha= 置信水平 Risklimits 危险比的可信区间 proc phreg data=ex17_5; model t*y(1)=x1-x6 /rl; run; proc phreg data=ex17_5; model t*y(1)=x1-x6 /rl; run; （九）应用Cox模型注意事项 对资料的分布类型没有要求； 样本量不宜过小，一般在40例以上，变量多时，要求是变量个数的10～15倍的阳性结局事件数； 要求因素的效应不随时间变化； 统计学上有关的因子，不一定有因果关系，可能有伴随关系； 因子之间可能存在相关性，因此结合专业知识和统计分析综合权衡； Cox回归的生存率不同资料间不宜相互比较，因为基础风险函数只在同一份资料中保持相同； 生存曲线不能任意延长，也不能轻易用于预测； Cox模型与logistic回归模型的关系 两者都可以估计相应因素的相对危险度；但logistic回归只考虑了事件的结果且作为二分类变量，而没有考虑生存时间的长短。 Cox回归模型的分析步骤： 模型的选取，包括确定自变量； PH假定条件验证； 单因素分析； 多因素分析； 模型诊断（模型拟合优度、残差分析、相关性）； 统计结论。 五、小结 * 在临床随访研究中，我们着重搜集病人出现某种结果（如痊愈、复发或死亡）所经历的时间，以便比较不同的治疗措施的远期效应。 * SAS程序 Proc lifetest options; Time variable*censor (删失值); Strata variable; Freq 变量； Run; 选择项 Method=PL /LT KM LIFE缺省为PL法 在LT法中，可规定 Intervals=一组数字 规定时间段的起点 Intervals=(o to n by n2) Width=value 寿命表计算时的区间宽度 Plots=(s, ls , lls，h) 生存曲线的比较总结： log-rank检验 将两组非删失时间混合从小到大排序，得多个四格表，aj和ej分别为第j个四格表中某组死亡数的实际频数和理论频数， vj为aj的方差。 ? Group 1 Group 2 ? Event No Event a b c d Wilcoxon检验（SPSS中称为Breslow检验） aj和ej分别为第j个四格表中某组死亡数的实际频数和理论频数， vj为aj的方差。Nj为第j个四格表的例数。 log-rank检验与Wilcoxon检验的公式相似，但log-rank检验权重为1，而Wilcoxon检验的权重为Nj 。 Log-rank检验和Wilcoxon检验统计量的一般形式： 似然比检验 其中Nj是第j组的结局例数，Tj是第j组的总时间。 似然比检验与log-rank检验、Wilcoxon检验不同，通常用于服从指数分布的数据。 三种检验方法的选择： 大样本情况下，三种方法的结果是一致的； 文献中最常报道的是log-rank检验； 三种检验结果不一致时，需结合生存曲线、对数生存曲线、二次对数生存曲线等来判断。 生存曲线(S)图：以生存率估计值为纵坐标，时间为横坐标。 对数生存曲线(LS)图：以生存率估计值的对数为纵坐标，以时间为横坐标。 二次对数生存曲线(LLS)图：以生存率估计值的对数的对数为纵坐标，以时间的对数为横坐标。 三种检验方法结果不同时选哪个？ 似然比检验：通常只用于LS图近似为直线，也就是服从指数分布的情形，其它情况下效率较低。 Log-rank检验：当LLS图中曲线近似为直线时，或S图的两条曲线大致平行时，效率较高。 Wilcoxon检验：当S图中的曲线有交叉时， Wilcoxon检验的效率高于log-rank检验和似然比检验。 Logrank检验与Wilcoxon检验比较： Wilcoxon检验对生存时间较短的个体赋予较大权重，更容易检验出早期的差异。 Logrank检验对所有死亡时间赋予相等的权重，对后期曲线差别的检验更为敏感。 如果log-rank检验有意义而Wilcoxon检验无意义，表明远期的生存差异较大，近