生存分析.docx

1、生存分析SPSS Survival(生存分析)SPSS Survival菜单包括Life Tables过程、Kaplan-Meier(卡普兰-梅尔)过程、Cox Regression过程、Cox w/Time-Dep Cov(含时间依存变量的Cox模型)过程。这里只介绍Life Tables过程和Kaplan-Meier过程。一、Kaplan-Meier过程采用乘积极限法（Product-limit estimates）来估计生存率，同时还可以对一个因素进行检验。适用于以个体为单位来收的小样本或大样本且有精确生存时间的生存资料，是最基本的一种生存分析方法。 Kaplan-Meier法用于：1、

2、估计某研究因素不同水平的中位生存时间。2、比较该研究因素不同水平的生存时间有无差异。3、控制一个分层因素后对研究因素不同水平的生存时间比较（此时将按分层因素的不同水平对研究因素对生存时间的影响分别进行分析）。二、Life Tables过程采用寿命表法（Life Tables）来估计生存率。适用于将整个观测时间划分为很多小的时间段，对于每个时间段，计算所有活到某时间段起点的病例在该时间段内死亡（出现结局）的概率。因此，当资料是按照固定的时间间隔收集（比如一个月随访一次）时，随访结果只有该年或该月期间的若干观察人数、发生失效事件人数（即出现预期观察结果的人数）和截尾人数，每位患者的确切生存时间无法

3、知道，此时采用寿命表法来分析。它用于大样本或粗略的存在时间资料。Life Tables过程用于：1、估计某生存时间的生存率。如估计术后1年生存率、2年生存率、3年生存率、4年生存率、5年生存率，等等。2、绘制各种曲线如生存函数、风险函数曲线等。3、对某一研究因素不同水平的生存时间分布进行比较，控制另一因素后对研究因素不同水平的生存时间分布进行比较，包括从总体上比较和不同水平之间进行两两比较。练习，P460第1题（1）求各年生存概率：各年生存概率=某年活满一年人数/某年年初人口数（2）求逐年生存率：某年生存率=某年生存概率与该年前各年生存概率的乘积结果为：术后年份死亡数年初人口数生存数生存概率生

4、存率1（0）10100900.9000 0.9000 2（1）1090800.8890 0.8001 3（2）1080700.8750 0.7009 作业：P460第2题题依分析：小样本且已知精确生存时间，比较两种疗法的生存时间有无差异。故采用Kaplan-Meier法。H0：两总体生存率相同（或两总体生存曲线相同）H1：两总体生存率不相同（或两总体生存曲线不相同）a=0.05SPSS过程:1、对变量进行定义2、输入数据timestatusgroup5117011311131123113011300138114211421145011123123127121012151215122312301

5、21.Analyze=Survival =Kaplan-Meier2.Time框：选入 “time”3.Status框：选入“status”； 击define events钮,在single value框右边的空格中输入“1”（0=“截尾或生存”，1=“死亡”等阳性结果）4.Factor框：选入“group”5. Compare factors列表框（分组因素水平间比较）： Test Statistics:选择Log- rank、Breslow、Tarone-Ware Linear trend for factor levels:选Pooled over strata（水平间的整体比较）6.

6、Option列表框 Statistics: 选Survival table(s)、Mean and median Survival Plots: 选Survival单击OK钮结果解释:1、生存分析表2、生存时间估计甲疗法生存时间均数为29.520周，中位生存时间为38.000天；乙疗法生存时间均数为11.889周，中位生存时间为10.000天。3、水平间的整体比较（即对整条生存曲线进行比较）3种检验方法的检验统计量 P 均小于0.05，2种疗法的生存率的差异有统计学意义。习惯取log-rank统计量。所以，本例2=7.628，P=0.006。因PSurvival =Life Tables3、T

7、ime框：选入time4、Display Time Intervals框（输出生存时间范围及组距）：在by前面的框内填入生存时间上限“20”（此区间必须包括生存时间的最大值）；在by后面的框内填入生存时间的组距（即时间间隔）“5”。5、Status框：选入status；击define events钮,在single value框右边的空格中输入“1”（死亡）6、单击Option按钮，弹出对话框： Life Table(s) 输出寿命表，系统默认 Plots: 选Survival（累积生存函数曲线）击Continue7、单击OK钮结果解释:1、生存分析表Life Table(a)Interval

8、 Start TimeNumber Entering IntervalNumber Withdrawing during IntervalNumber Exposed to RiskNumber of Terminal EventsProportion TerminatingProportion SurvivingCumulative Proportion Surviving at End of IntervalStd. Error of Cumulative Proportion Surviving at End of IntervalProbability DensityStd. Erro

9、r of Probability DensityHazard RateStd. Error of Hazard Rate0949189407310.780.220.220.010.1560.0030.250.01520016192520.270.730.160.010.0120.0020.060.01101327594.5140.150.850.140.010.0050.0010.030.0115433326.5100.380.620.090.010.010.0030.090.03The median survival time is 3.21寿命表生存时间的组段下限进入该组段的观察例数该组段

10、的截尾例数。暴露于危险因素的例数，即有效观察例数终结事件的例数，即死亡例数死亡比例生存比例至本组段上限的累积生存率累积生存率标准误概率密度概率密度标准误风险率风险率标准误0949189407310.780.220.220.010.1560.0030.250.01520016192520.270.730.160.010.0120.0020.060.01101327594.5140.150.850.140.010.0050.0010.030.0115433326.5100.380.620.090.010.010.0030.090.03中位生存时间为 3.21年。所以，5年生存率为22%；10年生存率为16%；15年生存率为14%；20年生存率为9%。2、生存曲线