第十四章 生存分析
第十四章 生存分析 第一節(jié) Life Tables過程 14.1.1 主要功能 14.1.2 實例操作 第二節(jié) Kaplan-Meier過程 14.2.1 主要功能 14.2.2 實例操作 第三節(jié) Cox Regression過程 14.3.1 主要功能 14.3.2 實例操作 在臨床診療工作的評價中，慢性疾病的預(yù)后一般不適合用治愈率、病死率等指標(biāo)來考核 ，因為其無法在短時間內(nèi)明確判斷預(yù)后情況，為此，只能對患者進行長期隨訪，統(tǒng)計一 定時期后的生存或死亡情況以判斷診療效果。這就是生存分析。 第一節(jié) Life Tables過程 14.1.1 主要功能 調(diào)用此過程時，系統(tǒng)將采用即壽命表分析法，完成對病例隨訪資料在任意指定時點的生 存狀況評價。 返回目錄[pic] [pic]返回全書目錄 14.1.2 實例操作 ［例14-1］ 用中藥+化療（中藥組，16例）和單純化療（對照組，10例）兩種療法治療白血病患者后 ，隨訪記錄存活情況如下所示，試比較兩組的生存率。 |中藥組 |對照組 | |隨訪月數(shù) |是否死亡 |隨訪月數(shù) |是否死亡 | |10 |否 |2 |是 | |2 |是 |13 |否 | |12 |是 |7 |是 | |13 |否 |11 |是 | |18 |否 |6 |否 | |6 |是 |1 |否 | |19 |是 |11 |否 | |26 |否 |3 |否 | |9 |是 |17 |否 | |8 |是 |7 |否 | |6 |是 | | | |43 |是 | | | |9 |否 | | | |4 |否 | | | |31 |否 | | | |24 |否 | | | 14.1.2.1 數(shù)據(jù)準(zhǔn)備 激活數(shù)據(jù)管理窗口，定義變量名：隨訪月數(shù)的變量名為TIME，是否死亡的變量名為DEAT H，分組（即中藥組與對照組）的變量名為GROUP。輸入原始數(shù)據(jù)：隨訪月數(shù)按原數(shù)值； 是否死亡的，是為1，否為0；分組的，中藥組為1，對照組為2。 14.1.2.2 統(tǒng)計分析 激活Statistics菜單選Survival中的Life Tables...項，彈出Life Tables對話框（圖14.1）。從對話框左側(cè)的變量列表中選time，點擊(鈕使之進入time框 ；在Display Time Intervals欄中定義需要顯示生存率的時點，本例要求從0個月顯示至48個月，間隔為2個 月，故在0 through框中輸入48，在by框中輸入2。選death，點擊(鈕使之進入Status框，點擊Defi ne Event...鈕彈出Life Tables:Define Event for Status Variable對話框，在Single value欄中輸入1，表明death = 1為發(fā)生死亡事件者；點擊Continue鈕返回Life Tables對話框。選group，點擊(鈕使之進入Factor框，點擊Define Range...鈕，彈出Life Tables:Define Range for Factor Variable對話框，定義分組的范圍，在Mininum框中輸入1，在Maxinum框中輸入2，點擊 Continue鈕返回Life Tables對話框。 | | |[pic] | |圖14.1 生存資料的壽命表分析對話框 | 點擊Options...鈕彈出Life Tables: Options對話框，在Plot欄中選Survival項，要求繪制生存率曲線圖；在Compare Levels of First Factor欄中選Overall項，要求作組間生存狀況的比較。之后點擊Continue鈕返回Life Tables對話框，再點擊OK鈕即完成分析。 14.1.2.3 結(jié)果解釋 在結(jié)果輸出窗口中將看到如下統(tǒng)計數(shù)據(jù)： 共有26個觀察對象進入分析。系統(tǒng)先顯示中藥組（group = 1）的生存狀況壽命表，按用戶指定，從0月起，隔2個月直至42個月（原指定從0—48個月 ，但因42個月后，生存概率已為0，故42個月后至48個月的生存狀況不再顯示），分別顯 示進入該時點例數(shù)（Number Entrng this Intrvl）、從該時點失去的例數(shù)（Number Wdrawn Durong Intrvl）、該時點暴露于死亡危險的例數(shù)（Number Exposd to Risk）、該時點死亡的例數(shù)（Number of Termnl Events）、該時點死亡概率（Propn Terminating）、該時點生存概率（Propn Surviving）、該時點末生存率（Propn Surv at End）、單位時點的累積概率（Cumul Probability Densty）、該時點風(fēng)險比例（Hazard Rate）、生存率的標(biāo)準(zhǔn)誤（SE of Cumul Surviving）、單位時點累積概率的標(biāo)準(zhǔn)誤（SE of Probability Densty）、風(fēng)險比例的標(biāo)準(zhǔn)誤（SE of Hazard Rate）。如本例，用中藥+化療的方式治療白血病患者，至8個月時，死亡率為17.39%， 生存概率為82.61%，生存率為66.38%，風(fēng)險比例為9.52%。至42個月時，生存概率和生存 率均為0，此時風(fēng)險比例為100%。中藥組的50%生存率在19.44個月。 對照組同類結(jié)果的顯示，因在16個月時生存概率已為0，故僅從0月起，隔2個月至16個月 止。分析顯示，單純用化療，白血病患者的半數(shù)生存率約在16個月多一點，比中藥組少 三個月。 |This subfile contains: 26 observations | | | |Life Table | |Survival Variable TIME | |for GROUP = 1 | |Number Number Number Number Cumul SE of | |SE of | |Intrvl Entrng Wdrawn Exposd of Propn Propn Propn Proba- | |Cumul Proba- SE of | |Start this During to Termnl Termi- Sur- Surv bility | |Hazard Sur- bility Hazard | |Time Intrvl Intrvl Risk Events nating viving at End Densty| |Rate viving Densty Rate | |------ ------ ------ ------ ------ ------ ------ ------ | |------ ------ ------ ------ ------ | |.0 16.0 .0 16.0 .0 .0000 1.0000 1.0000 .0000 | |.0000 .0000 .0000 .0000 | |2.0 16.0 .0 16.0 1.0 .0625 .9375 .9375 .0313 | |.0323 .0605 .0303 .0322 | |4.0 15.0 1.0 14.5 .0 .0000 1.0000 .9375 .0000 | |.0000 .0605 .0000 .0000 | |6.0 14.0 .0 14.0 2.0 .1429 .8571 .8036 .0670 | |.0769 .1019 .0441 .0542 | |8.0 12.0 1.0 11.5 2.0 .1739 .8261 .6638 .0699 | |.0952 .1231 .0458 .0670 | |10.0 9.0 1.0 8.5 .0 .0000 1.0000 .6638 .0000 | |.0000 .1231 .0000 .0000 | |12.0 8.0 1.0 7.5 1.0 .1333 .8667 .5753 .0443 | |.0714 .1348 .0420 .0712 | |14.0 6.0 .0 6.0 .0 .0000 1.0000 .5753 .0000 | |.0000 .1348 .0000 .0000 | |16.0 6.0 .0 6.0 .0 .0000 1.0000 .5753 .0000 | |.0000 .1348 .0000 .0000 | |18.0 6.0 1.0 5.5 1.0 .1818 .8182 .4707 .0523 | |.1000 .1453 .0489 .0995 | |20.0 4.0 .0 4.0 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |22.0 4.0 .0 4.0 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |24.0 4.0 1.0 3.5 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |26.0 3.0 1.0 2.5 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |28.0 2.0 .0 2.0 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |30.0 2.0 1.0 1.5 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |32.0 1.0 .0 1.0 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |34.0 1.0 .0 1.0 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |36.0 1.0 .0 1.0 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |38.0 1.0 .0 1.0 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |40.0 1.0 .0 1.0 .0 .0000 1.0000 .4707 .0000 | |.0000 .1453 .0000 .0000 | |42.0 1.0 .0 1.0 1.0 1.0000 .0000 .0000 .2354 | |1.0000 .0000 .0727 .0000 | | | |The median survival time for these data is 19.44 | | | |Life Table | |Survival Variable TIME | |for GROUP = 2 | |Number Number Number Number Cumul SE of | |SE of | |Intrvl Entrng Wdrawn Exposd of Propn Propn Propn Proba- | |Cumul Proba- SE of | |Start this During to Termnl Termi- Sur- Surv bility | |Hazard Sur- bility Hazard | |Time Intrvl Intrvl Risk Events nating viving at End Densty| |Rate viving Densty Rate | |------ ------ ------ ------ ------ ------ ------ ------ | |------ ------ ------ ------ ------ | |.0 10.0 1.0 9.5 .0 .0000 1.0000 1.0000 .0000 | |.0000 .0000 .0000 .0000 | |2.0 9.0 1.0 8.5 1.0 .1176 .8824 .8824 .0588 | |.0625 .1105 .0553 .0624 | |4.0 7.0 .0 7.0 .0 .0000 1.0000 .8824 .0000 | |.0000 .1105 .0000 .0000 | |6.0 7.0 2.0 6.0 1.0 .1667 .8333 .7353 .0735 | |.0909 .1628 .0678 .0905 | |8.0 4.0 ...
第十四章 生存分析