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、估計(jì)某生存時(shí)間的生存率。2、繪制各種曲線如生存函數(shù)、風(fēng)險(xiǎn)函數(shù)曲線等。3、對(duì)某一研究因素不同水平的生存時(shí)間分布進(jìn)行比較，控制另一因素后對(duì)研究因素不同水平的生存時(shí)間分布進(jìn)行比較，包括從總體上比較和不同水平之間進(jìn)行兩兩比較。一、建立數(shù)據(jù)文件 定義兩個(gè)列變量：時(shí)間變量

2、： 取名“time”，label標(biāo)上“survival time(week)”。生存狀態(tài)變量：取名“status”，并賦值：0=“刪失”，1=“死亡”。二、操作過程 從菜單選擇1、Analyze=Survival =Life Tables2、Time框：選入time3、Display Time Intervals框：在by前面的框內(nèi)填入生存時(shí)間上限，本例填入20（此區(qū)間必須包括生存時(shí)間的最大值）；在by后面的框內(nèi)填入生存時(shí)間的組距，本例填入5，以保證結(jié)果列出“15-”的組段。4、Status框：選入status；擊define events鈕,在single value框右邊的空格中輸入15、

3、單擊Option按鈕，彈出對(duì)話框：lLife Table(s) 輸出壽命表，系統(tǒng)默認(rèn)lPlots: 選Survival（累積生存函數(shù)曲線）擊Continue6、單擊OK鈕附：界面說明圖1 壽命表主對(duì)話框【Time】框選入生存時(shí)間變量?！綝isplay Time Intervals】框欲輸出生存時(shí)間范圍及組距。在by前面的框內(nèi)填入生存時(shí)間上限，本例填入200（此區(qū)間必須包括生存時(shí)間的最大值）；在by后面的框內(nèi)填入生存時(shí)間的組距，本例填入20，以保證結(jié)果列出“100-”的組段?！維tatus】框選入生存狀態(tài)變量，并定義終結(jié)事件的標(biāo)記值。選入變量“Status”后，【Define Event】鈕被激

4、活變黑，擊該按鈕，彈出定義終結(jié)事件標(biāo)記值的對(duì)話框（圖1）。對(duì)二分類變量，一般以死亡、復(fù)發(fā)、惡化等表示終結(jié)事件。本例以死亡為終結(jié)事件，其標(biāo)記值為1，故在Single value 框內(nèi)填入1。擊Continue鈕。若生存狀態(tài)變量取值為一連續(xù)型變量，如反應(yīng)變量為收縮壓，則在下面的Range of values 框中輸入140 through 400，此處上限輸入400是我隨便輸入的一個(gè)上限，目的是為了定義高血壓患者，實(shí)際上恐怕沒有人的血壓能達(dá)到400，這樣才能保證包括所有的高血壓病例，具體情況具體分析。圖2 定義終結(jié)事件標(biāo)記值的對(duì)話框【Factor】框定義第1層因素，即分組因素。【By Factor

5、】框定義第2層因素，即分層因素?！綩ptions】選項(xiàng)擊Options按鈕，彈出選項(xiàng)對(duì)話框。（圖3）圖3 壽命表選項(xiàng)對(duì)話框 Life Table(s)：輸出壽命表，系統(tǒng)默認(rèn)。Plot：統(tǒng)計(jì)圖。 Survival：累積生存函數(shù)曲線。 Hazard：累積風(fēng)險(xiǎn)函數(shù)散點(diǎn)圖。 One minus survival：生存函數(shù)被1減后的曲線。 Log survival：對(duì)數(shù)累積生存函數(shù)曲線。 Density：密度函數(shù)散點(diǎn)圖。Compare Levels of First Factor：對(duì)第1層因素不同水平的比較，即主對(duì)話框（圖1）中的factor框中所選入的因素。 None：不做比較。系統(tǒng)默認(rèn)。 Overa

6、ll：整體比較。 Pairwise：兩兩比較。 壽命表各個(gè)指標(biāo)的意義Intrvl Start Time：生存時(shí)間的組段下限。Number Entrng this Intrvl：進(jìn)入該組段的觀察例數(shù)。Number Wdrawn During Intrval：該組段的刪失例數(shù)。Number Exposed to Risk：暴露于危險(xiǎn)因素的例數(shù)，即有效觀察例數(shù)（校正觀察人數(shù)）。Number of Termnl Events：終結(jié)事件的例數(shù)，即死亡例數(shù)。Propn Terminating：終結(jié)事件比例，即死亡比例。Propn Surviving：生存比例。Cumul Propn Surv at End

7、：至本組段上限的累積生存率。Probability Densty：概率密度。所有個(gè)體在時(shí)點(diǎn)t后單位時(shí)間內(nèi)死亡概率的估計(jì)值。Hazard Rate：風(fēng)險(xiǎn)率?；钸^時(shí)點(diǎn)t的個(gè)體，在時(shí)點(diǎn)t后單位時(shí)間內(nèi)死亡概率的估計(jì)值。SE of Cumul Surviving：累積生存率的標(biāo)準(zhǔn)誤。SE of Probability Densty：概率密度的標(biāo)準(zhǔn)誤。SE of Hazard Rate：風(fēng)險(xiǎn)率的標(biāo)準(zhǔn)誤。 Kaplan-Meier過程采用乘積極限法（Product-limit estimates）來估計(jì)生存率，同時(shí)還可以對(duì)一個(gè)因素進(jìn)行檢驗(yàn)。適用于以個(gè)體為單位來收的未分組生存資料，是最基本的一種生存分析方法。

8、Kaplan-Meier法用于：1、估計(jì)某研究因素不同水平的中位生存時(shí)間。2、比較該研究因素不同水平的生存時(shí)間有無差異。3、控制一分層因素后對(duì)研究因素不同水平的生存時(shí)間比較（此時(shí)將按分層因素的不同水平對(duì)研究因素對(duì)生存時(shí)間的影響分別進(jìn)行分析）。操作過程：1.Analyze=Survival =Kaplan-Meier2. Time框：選入 time3. Status框：選入status；擊define events鈕,在single value框右邊的空格中輸入14. Factor框：選入group；5. Compare factors列表框：l Test Statistics:選擇Log ra

9、nk、Breslow、Tarone-Warel Linear trend for factor levels:選Pooled over strata或Pairwise over strata6. Option列表框：l Statistics: 選Survival table(s)、Mean and median、Survivall Plots: 選Survival單擊OK鈕三、界面說明圖1 Kaplan-Meier法主對(duì)話框【Time】框選入生存時(shí)間變量。【Status】框選入生存狀態(tài)變量?！綟actor】框選入分組變量?！維trata】框選入分層變量?！綥ables Cases】框給個(gè)體標(biāo)記

10、。【Compare Factor】選項(xiàng)擊Compare Factor按鈕，彈出選項(xiàng)對(duì)話框。（圖2）圖2 分組因素水平間比較對(duì)話框l Test Statistics:檢驗(yàn)統(tǒng)計(jì)量。Log rank：檢驗(yàn)生存分布是否相同，各時(shí)間點(diǎn)權(quán)重一樣。Breslow：檢驗(yàn)生存分布是否相同，以各時(shí)間點(diǎn)的觀察例數(shù)為權(quán)重。Tarone-Ware：檢驗(yàn)生存分布是否相同，以各時(shí)間點(diǎn)的觀察例數(shù)的平方根為權(quán)重。Linear trend for factor levels：分組因素水平間的線性趨勢檢驗(yàn)。 Pooled over strata：水平間的整體比較。系統(tǒng)默認(rèn)。 For each stratum：按分層變量，對(duì)每一層進(jìn)

11、行分組因素各水平間的整體比較。 Pairwise over strata：分組因素各水平間的兩兩比較。 Pairwise for each stratum：按分層變量，對(duì)每一層進(jìn)行分組因素各水平間的兩兩比較?！維ave】 選項(xiàng)擊Save按鈕，彈出Save New Variables（保存新變量）對(duì)話框（圖7）。圖3 保存新變量對(duì)話框Survival：累積生存率估計(jì)。Standard error of survival：累積生存率估計(jì)的標(biāo)準(zhǔn)誤。Hazard：累積風(fēng)險(xiǎn)函數(shù)估計(jì)。Cumulative events：終結(jié)事件的累積頻數(shù)。在各水平內(nèi)，按生存時(shí)間和生存狀態(tài)排序?！綩ptions】選項(xiàng)擊Options按鈕,彈出選項(xiàng)對(duì)話框（圖8）。圖4 K-M法選項(xiàng)對(duì)話框lStatistics：統(tǒng)計(jì)量。Survival table(s)：生存分析表。Mean and median survival：平均生存時(shí)間和中位生存時(shí)間及其標(biāo)準(zhǔn)誤和可信區(qū)間。Q