SPSS生存分析過程

1、某公司在ASPSSSurViVa I住存分析)菜單SPSS SUrViv菜單包括 Life TabIes過程、KaPIan-Meie過程、COX RegreSSiQ過程、COX w/Time-Dep COV 過程。這里只介紹Life TabIeSS程和Kaplan-Meie過程。Life Table過程Life Table過程用于：1、估計某生存時間的生存率。2、繪制各種曲線如生存函數(shù)、風險函數(shù)曲線等。3、對某一研究因素不同水平的生存時間分布進行比較，控制另一因素后對研究因素不同水平的生存時間分布進行比較，包括從總體上比較和不同水平之間進行兩兩比較。一、建立數(shù)據(jù)文件定義兩個列變量：時間變量：取

2、名“ time label 標上“ SUrViVaI time(week)。 ”生存狀態(tài)變量：取名“ StatUS，并賦值：O= “刪失”，1= “死亡”。二、操作過程從菜單選擇1、Analyze=>Suwival =>Life TabIeS2、Time框：選入 time3、 DiSPlay Time Intervals框：在by前面的框內(nèi)填入生存時間上限，本例填入 20 （此區(qū)間必須包括 生存時間的最大值）；在by后面的框內(nèi)填入生存時間的組距，本例填入 5,以保證結(jié)果列出15-”的組 段。4、 StatuS匡：選入StatUS擊define eve nt鈕,在Si ngle Va

3、Iue框右邊的空格中輸入 15、單擊OPtiOn按鈕，彈出對話框：Life TabIe（S）輸出壽命表，系統(tǒng)默認Plots:選SUrViVaI （累積生存函數(shù)曲線）擊 Continue6、單擊OK鈕二 Lire TablesXOKPtee匚acelHeIP附：界面說明IITilnft! CFrDiSPIAy Time IntD0 LZ 5Iatus:5(1)Ftctor:Factor圖1壽命表主對話框Time框選入生存時間變量【DiSPlay Time Interva】框欲輸出生存時間范圍及組距C在by前面的框內(nèi)填入生存時間上限，本例填入200 （此區(qū)間必須包括生存時間的最大值）；在by后面的

4、框內(nèi)填入生存時間的組距，本例填入20,以保證結(jié)果列出“ 100的組段?！維tatUS框選入生存狀態(tài)變量，并定義終結(jié)事件的標記值。選入變量“ StatUS后，【Define EVent鈕被激活變黑，擊該按鈕，彈出定義終結(jié)事件標記值的對話 框（圖1）。對二分類變量，一般以死亡、復發(fā)、惡化等表示終結(jié)事件。本例以死亡為終結(jié)事件，其標 記值為1，故在Single value框內(nèi)填入1。擊COntinUe鈕。若生存狀態(tài)變量取值為一連續(xù)型變量，如反應 變量為收縮壓，則在下面的 Range of VaIUeSI中輸入140 through 400此處上限輸入400是我隨便輸入 的一個上限，目的是為了定義高血壓

5、患者，實際上恐怕沒有人的血壓能達到400,這樣才能保證包括所有的高血壓病例，具體情況具體分析。頁腳內(nèi)容13圖2定義終結(jié)事件標記值的對話框【Facto】框定義第1層因素，即分組因素?！綛y Facto框定義第2層因素，即分層因素?！綩ptionS選項擊OPtiOns按鈕，彈出選項對話框。（圖3）圖3壽命表選項對話框Life Table(S)輸出壽命表，系統(tǒng)默認。PlOt :統(tǒng)計圖。SUrViVat累積生存函數(shù)曲線。HaZard累積風險函數(shù)散點圖。One minus SUrViVaI生存函數(shù)被1減后的曲線。Log SUrViVaI對數(shù)累積生存函數(shù)曲線。DenSity密度函數(shù)散點圖。1)中的fact

6、or框中所COmPare LeVeIS of FirSt FaCt對第1層因素不同水平的比較，即主對話框(圖選入的因素。None不做比較。系統(tǒng)默認。Overall:整體比較。Pairwise:兩兩比較。壽命表各個指標的意義IntrVl Start Time生存時間的組段下限。Number Entrng this IntrVI進入該組段的觀察例數(shù)。Number WdraWn DUring Intrval 該組段的刪失例數(shù)。NUmber EXPOSed to RiSk暴露于危險因素的例數(shù)，即有效觀察例數(shù)（校正觀察人數(shù)）NUmber of Term nl EVe nts終結(jié)事件的例數(shù)，即死亡例數(shù)。P

7、ropn Terminating終結(jié)事件比例，即死亡比例。Propn SUrViVing 生存比例。Cumul Propn SUrV at EnC至本組段上限的累積生存率。PrObabiIity Den Sty概率密度。所有個體在時點t后單位時間內(nèi)死亡概率的估計值。HaZard Rate風險率?；钸^時點t的個體，在時點t后單位時間內(nèi)死亡概率的估計值。SE Of Cumul SUrViVing累積生存率的標準誤。SE of PrObabiIity De nsty概率密度的標準誤。SE Of HaZard Rate風險率的標準誤。KaP Ia n-Meie 過 程采用乘積極限法（PrOdUCt-I

8、imit estimates來估計生存率，同時還可以對一個因素進行檢驗。適用 于以個體為單位來收的未分組生存資料，是最基本的一種生存分析方法。KaP Ia n-Meie 法用于：1估計某研究因素不同水平的中位生存時間。2、比較該研究因素不同水平的生存時間有無差異。3、控制一分層因素后對研究因素不同水平的生存時間比較（此時將按分層因素的不同水平對研 究因素對生存時間的影響分別進行分析）。操作過程：1. An alyze=>Survival =>Kapla n-Meier2. Time框：選入 time3. StatuS框：選入StatuS擊define eventS鈕,在Single

9、 Value框右邊的空格中輸入 14. FaCtOr框：選入 group；5. COmPare factor列 表框：TeSt Statistics:擇 Log rank Breslow Tarone-WareLinear trend for factor IeVeI選 Pooled over Strata或 PairWiSe over Strata6. OPtiOn列表框：StatiStiCS選 SurViVal table(s) Mean and median SurViValPlots:選 SUrViVaI單擊OK鈕三、界面說明圖1 KaP Ia n-Meie法主對話框Time框選入生存

10、時間變量【StatuS框選入生存狀態(tài)變量【Facto】框選入分組變量。【Strata框選入分層變量?！綥ableS CaSeS框給個體標記?！綜OmPare FaCtOr 選項擊COmPare FactO按鈕，彈出選項對話框。（圖2）圖2分組因素水平間比較對話框TeSt Statistics：驗統(tǒng)計量。Log rank檢驗生存分布是否相同，各時間點權(quán)重一樣BreSlOW檢驗生存分布是否相同，以各時間點的觀察例數(shù)為權(quán)重。TarOne-Ware檢驗生存分布是否相同，以各時間點的觀察例數(shù)的平方根為權(quán)重。Linear trend for factor levels分組因素水平間的線性趨勢檢驗。Pool

11、ed over Strata水平間的整體比較。系統(tǒng)默認。FOreaCh StratUm按分層變量，對每一層進行分組因素各水平間的整體比較。PairWiSe over Strata分組因素各水平間的兩兩比較。PairWiSe for each StratU m按分層變量，對每一層進行分組因素各水平間的兩兩比較?！維aVe選項擊SaVe按鈕，彈出SaVe NeW VariabIeS保存新變量）對話框（圖 7）。圖3 保存新變量對話框SUrViVat累積生存率估計。Standard error of SUrViVaI累積生存率估計的標準誤。HaZard累積風險函數(shù)估計。CumuIatiVe eve nts終結(jié)事件的累積頻數(shù)。在各水平內(nèi)，按生存時間和生存狀態(tài)排序。【OptionS選項擊OPtiOns按鈕,彈出選項對話框(圖8)。圖4K-M法選項對話框StatiStiCS 統(tǒng)計量。SUrViVaI table(s)生存分析表。Me