SPSS生存分析过程

1、SPSS Survival（生存分析）菜单SPSS Survival 菜单包括 Life Tables 过程、Kaplan-Meier 过程、Cox Regression过程、Cox w/Time-Dep Cov 过程。这里只介绍 Life Tables 过 程和 Kaplan-Meier 过程。Life Tables 过程Life Tables过程用于：1、估计某生存时间的生存率。2、绘制各种曲线如生存函数、风险函数曲线等。3、对某一研究因素不同水平的生存时间分布进行比较，控制另一因素后对研究因素不同水平的生存时间分布进行比较，包括从总体上比较和不同水平之间进行两两比较。一、建立数据文件定义

2、两个列变虽：时间变虽：取名 “time ”，labe l 标上 a survival time（week）:生存状态变虽：取名“ status”，并赋值：0= “删失”，1= “死亡”。二、操作过程从菜单选择1、Analyze=>Survival =>Life Tables2、Time 框：选入 time3、Display Time Intervals框：在by前面的框内填入生存时间上限，本例填入20 （此区间必须包括生存时间的最大值）；在by后面的框内填入生 存时间的组距，本例填入 5,以保证结果列出 15-”的组段4、Status 框：选入 status;击 define ev

3、ents钮,在 single value框右边的空格中输入15、单击Option按钮，弹出对话框：Life Table(s)输出寿命表，系统默认Plots:选Survival (累积生存函数曲线)击 Continue6、单击OK钮附：界面说明j Life TabledOKFasteResetCancelHelpI ID i splay Timw Intervals0河 如SFactor:SABy Factor:Jpti tmM-图1寿命表主对话框Time框选入生存时间变虽【Display Time Intervals 】框欲输出生存时间范围及组距。在by前面的框内填入生存时间上限，本例填入20

4、0 （此区间必须包括生存时间的最大值）；在by后面的框内填入生存时间的组距，本例填入20,以保证结果列出 “ 100”的组段。【Status】框选入生存状态变虽，并定义终结事件的标记值。选入变H “ Status”后，【Define Event钮被激活变黑，击该按钮，弹出 定义终结事件标记值的对话框（图1）。对二分类变虽，一般以死亡、复发、恶化等表示终结事件。 本例以死亡为终结事件，其标记值为1,故在Singlevalue框内填入1。击Continue钮。若生存状态变H取值为一连续型变H, 如反应变虽为收缩压，则在下面的Range of values框中输入140 through400,此处上

5、限输入 400是我随便输入的一个上限，目的是为了定义高血 压患者，实际上恐怕没有人的血压能达到400,这样才能保证包括所有的高血压病例，具体情况具体分析。图2定义终结事件标记值的对话框【Factor框定义第1层因素，即分组因素【By Factor框定义第2层因素，即分层因素。【Options选项击Options按钮，弹出选项对话框。（图3）图3寿命表选项对话框Life Table（s）:输出寿命表，系统默认。Plot:统计图。Survival :累积生存函数曲线。Hazard :累积风险函数散点图。One minus survival :生存函数被 1减后的曲线。Log survival:对数

6、累积生存函数曲线。Density:密度函数散点图。Compare Levels of First Factor :对第1层因素不同水平的比较，即主对话框（图1）中1的factor框中1所选入的因素。None:不做比较。系统默认。Overall :整体比较。Pairwise:两两比较。寿命表各个指标的意义Intrvl Start Time :生存时间的组段下限。Number Entrng this Intrvl :进入该组段的观察例数。Number Wdrawn During Intrval :该组段的删失例数。Number Exposed to Risk :暴露于危险因素的例数，即有效观察例数

7、（校 正观察人数）。Number of Termnl Events :终结事件的例数，即死亡例数。Propn Terminating :终结事件比例，即死亡比例。Propn Surviving :生存比例。Cumul Propn Surv at End :至本组段上限的累积生存率。Probability Densty :概率密度。所有个体在时点t后单位时间内死亡概率的估计值。Hazard Rate :风险率。活过时点t的个体，在时点t后单位时间内死亡概 率的估计值。SE of Cumul Surviving :累积生存率的标准误。SE of Probability Densty :概率密度的标

8、准误。SE of Hazard Rate :风险率的标准误。Kaplan-Meier 过程采用乘积极限法(Product-limit estimates)来估计生存率，同时还可以对一个因素进行检验。适用于以个体为单位来收的未分组生存资料，是最 基本的一种生存分析方法。Kaplan-Meier 法用于：1、估计某研究因素不同水平的中位生存时间。2、比较该研究因素不同水平的生存时间有无差异。3、控制一分层因素后对研究因素不同水平的生存时间比较(此时将按分层因素的不同水平对研究因素对生存时间的影响分别进行分析)。操作过程：1. Analyze=>Survival =>Kaplan-Mei

9、er2. Time 框：选入 time3. Status 框： 选入 status; 击 define events钮，在 single value框右边的 空格中输入14. Factor 框：选入 group;5. Compare factors 列表框：Test Statistics选择 Log rank、Breslow、Tarone-WareLinear trend for factor levels:选 Pooled over strata或 Pairwise overstrata6. Option列表框：Statistics:选 Survival table(s)、Mean and

10、median、SurvivalPlots:选 Survival单击OK钮三、界面说明图1 Kaplan-Meier 法主对话框Time框选入生存时间变虽【Status】框选入生存状态变虽【Factor框选入分组变H。【Strata 框选入分层变虽。【Lables Cases框给个体标记。【Compare Factor 】 选项击Compare Factor按钮，弹出选项对话框。（图2）图2分组因素水平间比较对话框Test Statistics 检验统计 H。Log rank :检验生存分布是否相同，各时间点权重一样。Breslow:检验生存分布是否相同，以各时间点的观察例数为权重。Tarone

11、-Ware :检验生存分布是否相同，以各时间点的观察例数的平方 根为权重。Linear trend for factor levels :分组因素水平间的线性趋势检验。Pooled over strata:水平间的整体比较。系统默认。For each stratum :按分层变H,对每一层进行分组因素各水平间的整 体比较。Pairwise over strata :分组因素各水平间的两两比较。Pairwise for each stratum :按分层变对每一层进行分组因素各水平间的两两比较【Save】选项击Save按钮，弹出Save New Variables （保存新变虽）对话框（图7）图3保存新变H对话框Survival :累积生存率估计。Standard error of survival :累积生存率估计的标准误。Hazard :累积风险函数估计。Cumulative events :终结事件的累积频数。在各水平内，按生存时间和 生存状态排序。【Options选项击Options按钮，弹出选项对话框（图 8）。图4 K-M法选项对话框Statistics*:统计H。Survival table(s):生存分析表。Mean a