第十四章生存分析

第十四章生存分析第一节LifeTables过程14.1.1主要功能14.1.2实例操作第二节Kaplan-Meier过程14.2.1主要功能14.2.2实例操作第三节CoxRegression过程14.3.1主要功能14.3.2实例操作人治愈率、病死率等指标来考核，因为其无法在短时间内明确判断预后情况，为此，只能对患者进行长期随访，统计一定时期后的生存或死亡情况以判断诊疗效果。这就是生存分析。第一节LifeTables过程14.1.1主要功能调用此过程时，系统将采用即寿命表分析法，完成对病例随访资料在任意指定时点的生存状况评价。返回目录返回全书目录14.1.2实例操作［例14-1］用中药+化疗（中药组，16例）和单纯化疗（对照组，10例）两种疗法治疗白血病患者后，随访记录存活情况如下所示，试比较两组的生存率。中药组对照组随访月数是否死亡随访月数是否死亡1021213186192698643943124否是是否否是是否是是是是否否否否21371161113177是否是是否否否否否否14.1.2.1数据准备激活数据管理窗口，定义变量名：随访月数的变量名为TIME，是否死亡的变量名为DEATH，分组（即中药组与对照组）的变量名为GROUP。输入原始数据：随访月数按原数值；是否死亡的，是为1，否为0；分组的，中药组为1，对照组为2。14.1.2.2统计分析激活Statistics菜单选Survival中的LifeTables...项，弹出LifeTables对话框（图14.1）。从对话框左侧的变量列表中选time，点击钮使之进入time框；在DisplayTimeIntervals栏中定义需要显示生存率的时点，本例要求从0个月显示至48个月，间隔为2个月，故在0through框中输入48，在by框中输入2。选death，点击钮使之进入Status框，点击DefineEvent...钮弹出LifeTables:DefineEventforStatusVariable对话框，在Singlevalue栏中输入1，表明death=1为发生死亡事件者；点击Continue钮返回LifeTables对话框。选group，点击钮使之进入Factor框，点击DefineRange...钮，弹出LifeTables:DefineRangeforFactorVariable对话框，定义分组的范围，在Mininum框中输入1，在Maxinum框中输入2，点击Continue钮返回LifeTables对话框。图14.1生存资料的寿命表分析对话框点击Options...钮弹出LifeTables:Options对话框，在Plot栏中选Survival项，要求绘制生存率曲线图；在CompareLevelsofFirstFactor栏中选Overall项，要求作组间生存状况的比较。之后点击Continue钮返回LifeTables对话框，再点击OK钮即完成分析。14.1.2.3结果解释在结果输出窗口中将看到如下统计数据：共有26个观察对象进入分析。系统先显示中药组（group=1）的生存状况寿命表，按用户指定，从0月起，隔2个月直至42个月（原指定从0—48个月，但因42个月后，生存概率已为0，故42个月后至48个月的生存状况不再显示），分别显示进入该时点例数（NumberEntrngthisIntrvl）、从该时点失去的例数（NumberWdrawnDurongIntrvl）、该时点暴露于死亡危险的例数（NumberExposdtoRisk）、该时点死亡的例数（NumberofTermnlEvents）、该时点死亡概率（PropnTerminating）、该时点生存概率（PropnSurviving）、该时点末生存率（PropnSurvatEnd）、单位时点的累积概率（CumulProbabilityDensty）、该时点风险比例（HazardRate）、生存率的标准误（SEofCumulSurviving）、单位时点累积概率的标准误（SEofProbabilityDensty）、风险比例的标准误（SEofHazardRate）。如本例，用中药+化疗的方式治疗白血病患者，至8个月时，死亡率为17.39%，生存概率为82.61%，生存率为66.38%，风险比例为9.52%。至42个月时，生存概率和生存率均为0，此时风险比例为100%。中药组的50%生存率在19.44个月。对照组同类结果的显示，因在16个月时生存概率已为0，故仅从0月起，隔2个月至16个月止。分析显示，单纯用化疗，白血病患者的半数生存率约在16个月多一点，比中药组少三个月。Thissubfilecontains:26observationsLifeTableSurvivalVariableTIMEforGROUP=1NumberNumberNumberNumberCumulSEofSEofIntrvlEntrngWdrawnExposdofPropnPropnPropnProba-CumulProba-SEofStartthisDuringtoTermnlTermi-Sur-SurvbilityHazardSur-bilityHazardTimeIntrvlIntrvlRiskEventsnatingvivingatEndDenstyRatevivingDenstyRate------------------------------------------------------------------------------.016.0.016.0.0.00001.00001.0000.0000.0000.0000.0000.00002.016.0.016.01.0.0625.9375.9375.0313.0323.0605.0303.03224.015.01.014.5.0.00001.0000.9375.0000.0000.0605.0000.00006.014.0.014.02.0.1429.8571.8036.0670.0769.1019.0441.05428.012.01.011.52.0.1739.8261.6638.0699.0952.1231.0458.067010.09.01.08.5.0.00001.0000.6638.0000.0000.1231.0000.000012.08.01.07.51.0.1333.8667.5753.0443.0714.1348.0420.071214.06.0.06.0.0.00001.0000.5753.0000.0000.1348.0000.000016.06.0.06.0.0.00001.0000.5753.0000.0000.1348.0000.000018.06.01.05.51.0.1818.8182.4707.0523.1000.1453.0489.099520.04.0.04.0.0.00001.0000.4707.0000.0000.1453.0000.000022.04.0.04.0.0.00001.0000.4707.0000.0000.1453.0000.000024.04.01.03.5.0.00001.0000.4707.0000.0000.1453.0000.000026.03.01.02.5.0.00001.0000.4707.0000.0000.1453.0000.000028.02.0.02.0.0.00001.0000.4707.0000.0000.1453.0000.000030.02.01.01.5.0.00001.0000.4707.0000.0000.1453.0000.000032.01.0.01.0.0.00001.0000.4707.0000.0000.1453.0000.000034.01.0.01.0.0.00001.0000.4707.0000.0000.1453.0000.000036.01.0.01.0.0.00001.0000.4707.0000.0000.1453.0000.000038.01.0.01.0.0.00001.0000.4707.0000.0000.1453.0000.000040.01.0.01.0.0.00001.0000.4707.0000.0000.1453.0000.000042.01.0.01.01.01.0000.0000.0000.23541.0000.0000.0727.0000Themediansurvivaltimeforthesedatais19.44LifeTableSurvivalVariableTIMEforGROUP=2NumberNumberNumberNumberCumulSEofSEofIntrvlEntrngWdrawnExposdofPropnPropnPropnProba-CumulProba-SEofStartthisDuringtoTermnlTermi-Sur-SurvbilityHazardSur-bilityHazardTimeIntrvlIntrvlRiskEventsnatingvivingatEndDenstyRatevivingDenstyRate------------------------------------------------------------------------------.010.01.09.5.0.00001.00001.0000.0000.0000.0000.0000.00002.09.01.08.51.0.1176.8824.8824.0588.0625.1105.0553.06244.07.0.07.0.0.00001.0000.8824.0000.0000.1105.0000.00006.07.02.06.01.0.1667.8333.7353.0735.0909.1628.0678.09058.04.0.04.0.0.00001.0000.7353.0000.0000.1628.0000.000010.04.01.03.51.0.2857.7143.5252.1050.1667.2122.0918.164312.02.01.01.5.0.00001.0000.5252.0000.0000.2122.0000.000014.01.0.01.0.0.00001.0000.5252.0000.0000.2122.0000.000016.01.01.0.5.0.00001.0000.5252.0000.0000.2122.0000.0000Themediansurvivaltimeforthesedatais16.00+接着显示两组比较的结果。系统采用Gehan比分检验法，得u=0.012，P=0.9113，即中药组与对照组的生存率无差别。ComparisonofsurvivalexperienceusingtheWilcoxon(Gehan)statisticSurvivalVariableTIMEgroupedbyGROUPOverallcomparisonstatistic.012D.F.1Prob..9113GrouplabelTotalNUncenCenPctCenMeanScore1168850.00.18752103770.00-.3000最后，系统输出生存率曲线图（图14.2）。从图中可见，对照组（group=2）在8个月前一段时点的生存率均较中药组（group=1）略低，而8-12个月这一段其生存率又较中药组略高，12个月后再又下降。但在治疗中加用中药，对个别患者而言，20个月后依然有一定的生存率。图14.2中药组与对照组生存率曲线的比较返回目录返回全书目录第二节Kaplan-Meier过程14.2.1主要功能调用此过程，系统将采用Kaplan-Meier方法，对病例随访资料进行生存分析，在对应于每一实际观察事件时点上，作生存率的评