17 生存分析-SPSS

2020/2/281生存分析的SPSS实现预防系卫生统计学教研室张彦琦Tel:023-6875234313996478364E-mail:[e-mail]2020/2/282SPSS生存分析过程名寿命表法Analyze→Survival→LifeTablesKaplan-Meier法Analyze→Survival→Kaplan-MeierCox回归分析Analyze→Survival→CoxRegression2020/2/283生存分析寿命表法Kaplan-Meier法Cox回归2020/2/284寿命表法2020/2/285实例分析生存率估计——寿命表法例17-1有人研究了甲种手术方法治疗某病的生存情况，定义从手术后到死亡为生存时间，得到的生存时间（月）如下，其中有“+”者表示截尾数据，表示仍生存或失访，括号内为重复死亡数，试作生存分析。135（3）6（3）7810（3）14+1719+20+22+26+31+3434+44592020/2/286数据格式23行2列：生存时间、结局（1：完全数据、2：删失数据）2020/2/287操作步骤Analyze→Survival→LifeTables2020/2/288LifeTables对话框生存时间结局变量定义生存时间范围及组距定义结局变量2020/2/289LifeTables:DefineEvent对话框2020/2/2810LifeTables:Options对话框寿命表生存曲线风险函数曲线生存曲线被1减后的曲线2020/2/2811结果输出2020/2/2812寿命表LifeTablea23023.0005.22.78.78.09.036.014.04.0218018.0007.39.61.48.10.051.016.08.0311110.5001.10.90.43.10.008.007.02.02937.5000.001.00.43.10.000.000.00.00615.5000.001.00.43.10.000.000.00.00524.0001.25.75.32.12.018.016.05.05202.0000.001.00.32.12.000.000.00.00202.0001.50.50.16.13.027.022.11.10101.0000.001.00.16.13.000.000.00.00101.00011.00.00.00.00.027.022.33.00IntervalStartTime061218243036424854NumberEnteringIntervalNumberWithdrawingduringIntervalNumberExposedtoRiskNumberofTerminalEventsProportionTerminatingProportionSurvivingCumulativeProportionSurvivingatEndofIntervalStd.ErrorofCumulativeProportionSurvivingatEndofIntervalProbabilityDensityStd.ErrorofProbabilityDensityHazardRateStd.ErrorofHazardRateThemediansurvivaltimeis11.5714a.2020/2/2813LifeTablea23023.0005.22.78.78.09.036.014.04.0218018.0007.39.61.48.10.051.016.08.0311110.5001.10.90.43.10.008.007.02.02937.5000.001.00.43.10.000.000.00.00615.5000.001.00.43.10.000.000.00.00524.0001.25.75.32.12.018.016.05.05202.0000.001.00.32.12.000.000.00.00202.0001.50.50.16.13.027.022.11.10101.0000.001.00.16.13.000.000.00.00101.00011.00.00.00.00.027.022.33.00IntervalStartTime061218243036424854NumberEnteringIntervalNumberWithdrawingduringIntervalNumberExposedtoRiskNumberofTerminalEventsProportionTerminatingProportionSurvivingCumulativeProportionSurvivingatEndofIntervalStd.ErrorofCumulativeProportionSurvivingatEndofIntervalProbabilityDensityStd.ErrorofProbabilityDensityHazardRateStd.ErrorofHazardRateThemediansurvivaltimeis11.5714a.寿命表生存时间组段下限初期观察例数期内删失例数有效观察例数期内死亡例数条件死亡率中位生存时间为11.57月。2020/2/2814LifeTablea23023.0005.22.78.78.09.036.014.04.0218018.0007.39.61.48.10.051.016.08.0311110.5001.10.90.43.10.008.007.02.02937.5000.001.00.43.10.000.000.00.00615.5000.001.00.43.10.000.000.00.00524.0001.25.75.32.12.018.016.05.05202.0000.001.00.32.12.000.000.00.00202.0001.50.50.16.13.027.022.11.10101.0000.001.00.16.13.000.000.00.00101.00011.00.00.00.00.027.022.33.00IntervalStartTime061218243036424854NumberEnteringIntervalNumberWithdrawingduringIntervalNumberExposedtoRiskNumberofTerminalEventsProportionTerminatingProportionSurvivingCumulativeProportionSurvivingatEndofIntervalStd.ErrorofCumulativeProportionSurvivingatEndofIntervalProbabilityDensityStd.ErrorofProbabilityDensityHazardRateStd.ErrorofHazardRateThemediansurvivaltimeis11.5714a.寿命表条件生存率累计生存率累计生存率的标准误概率密度函数概率密度函数的标准误风险率风险率标准误2020/2/2815生存曲线中位生存时间2020/2/2816Kaplan-Meier法2020/2/2817实例分析生存率比较——Kaplan-Meier法例17-2假定我们用乙种手术方式治疗了与例17-1病情一致的同种疾病20例。试比较两种手术方式的生存时间。135（3）6（3）7810（3）14+1719+20+22+26+31+3434+44591（2）23（2）4（3）6（2）89（2）1011121314151718甲种手术乙种手术2020/2/2818数据格式43行3列：生存时间、结局（1：完全数据、2：删失数据）、分组（1：甲组、2：乙组）2020/2/2819操作步骤Analyze→Survival→Kaplan-Meier2020/2/2820Kaplan-Meier对话框生存时间结局变量定义结局变量分组变量组间比较2020/2/2821Kaplan-Meier:DefineEvent对话框2020/2/2822Kaplan-Meier:CompareFactorLevels对话框对数秩检验2020/2/2823Kaplan-Meier：Options对话框寿命表平均和中位生存时间生存曲线2020/2/2824结果输出2020/2/2825基本统计描述CaseProcessingSummary2316730.4%20200.0%4336716.3%分组甲组乙组OverallTotalNNofEventsNPercentCensored例数死亡例数删失例数删失率2020/2/2826寿命表SurvivalTable1.0001.00.957.0431223.0001.00.913.0592215.0001.00..3205.0001.00..4195.0001.00.783.0865186.0001.00..6176.0001.00..7166.0001.00.652.0998157.0001.00.609.1029148.0001.00.565.103101310.0001.00..111210.0001.00.478.104121114.000.00..121017.0001.00.430.10413919.000.00..13820.000.00..13722.000.00..13626.000.00..13531.000.00..13434.0001.00.323.12214334.000.00..14244.0001.00.161.12915159.0001.00.000.0001601.0001.00..1191.0001.00.900.0672182.0001.00.850.0803173.0001.00..4163.0001.00.750.0975154.0001.00..6144.0001.00..7134.0001.00.600.1108126.0001.00..9116.0001.00.500.11210108.0001.00.450.1111199.0001.00..1289.0001.00.350.10713710.0001.00.300.10214611.0001.00.250.09715512.0001.00.200.08916413.0001.00.150.08017315.0001.00.100.06718217.0001.00.050.04919118.0001.00.000.00020012345678910111213141516171819202122231234567891011121314151617181920分组甲组乙组TimeStatusEstimateStd.ErrorCumulativeProportionSurvivingattheTimeNofCumulativeEventsNofRemainingCases生存时间结局累积生存率及其标准误累计死亡例数存活例数2020/2/2827寿命表SurvivalTable1.0001.00.957.0431223.0001.00.913.0592215.0001.00..3205.0001.00..4195.0001.00.783.0865186.0001.00..6176.0001.00..7166.0001.00.652.0998157.0001.00.609.1029148.0001.00.565.103101310.0001.00..111210.0001.00.478.104121114.000.00..121017.0001.00.430.10413919.000.00..13820.000.00..13722.000.00..13