醫(yī)學(xué)統(tǒng)計(jì)學(xué)：生存分析SPSS

生存分析SPSS過程

（SPSSofSurvivalAnalysis）

1.何為生存分析？

生存分析（survivalanalysis）是將事件的結(jié)果（終點(diǎn)事件）和出現(xiàn)結(jié)果經(jīng)歷的時(shí)間結(jié)合起來分析的一種統(tǒng)計(jì)分析方法。2.生存分析的目的：描述生存過程：估計(jì)不同時(shí)間的總體生存率，計(jì)算中位生存期，繪制生存函數(shù)曲線。統(tǒng)計(jì)方法包括Kaplan-Meier（K-M）法、壽命表法。比較：比較不同處理組的生存率，如比較不同療法治療腦瘤的生存率，以了解哪種治療方案較優(yōu)。統(tǒng)計(jì)方法log-rank檢驗(yàn)等。影響因素分析：研究某個(gè)或某些因素對(duì)生存率或生存時(shí)間的影響作用。如為改善腦瘤病人的預(yù)后，應(yīng)了解影響病人預(yù)后的主要因素，包括病人的年齡、性別、病程、腫瘤分期、治療方案等。統(tǒng)計(jì)方法cox比例風(fēng)險(xiǎn)回歸模型等。預(yù)測(cè)：建立cox回歸預(yù)測(cè)模型。生存分析的理論復(fù)習(xí)生存分析(SurvivalAnalysis)菜單壽命表（LifeTables）過程Lifetables過程用于（小樣本和大樣本資料）：估計(jì)某生存時(shí)間的生存率，以及中位生存時(shí)間。繪制各種曲線：如生存函數(shù)、風(fēng)險(xiǎn)函數(shù)曲線等。對(duì)某一研究因素不同水平的生存時(shí)間分布的比較?？刂屏硪粋€(gè)因素后對(duì)研究因素不同水平的生存時(shí)間分布的比較。對(duì)多組生存時(shí)間分布進(jìn)行兩兩比較。（比較總體生存時(shí)間分布采用wilcoxon檢驗(yàn)）CompanyLogo實(shí)例分析例1：為了比較不同手術(shù)方法治療腎上腺腫瘤的療效，某研究者隨機(jī)將43例病人分成兩組，甲組23例、乙組20例的生存時(shí)間（月）如下所示：其中有“+”者是刪失數(shù)據(jù)，表示病人仍生存或失訪，括號(hào)內(nèi)為死亡人數(shù)。（1）計(jì)算甲、乙兩法術(shù)后10月的生存率和標(biāo)準(zhǔn)誤。（2）估計(jì)兩組的中位生存期。（3）繪制各組生存函數(shù)曲線。（4）比較兩組的總體生存時(shí)間分布有無差別。CompanyLogo一、建立數(shù)據(jù)文件（data-01.sav）定義5個(gè)變量：生存時(shí)間變量：t，值標(biāo)簽“生存時(shí)間（月）”生存狀態(tài)變量：status，取值“1=死亡，0=刪失或存活”頻數(shù)變量：freq，值標(biāo)簽“人數(shù)”分組變量：group，取值“1=甲組，2=乙組”生存時(shí)間序號(hào)變量(可無）：i頻數(shù)格式的數(shù)據(jù)文件，生存分析之前需要采用“加權(quán)個(gè)案”過程指定頻數(shù)變量：二、生存分析操作過程主菜單：分析Analyze

生存Survival壽命表Lifetables對(duì)話框參數(shù)設(shè)置：時(shí)間time框：選入“t”。顯示時(shí)間間隔Displaytimeintervals框：步長(zhǎng)by前面填入最大生存時(shí)間的上限（必須包括生存時(shí)間最大值），步長(zhǎng)by后面填入生存時(shí)間的組距。本例上限填“60”，組距填“1”。狀態(tài)status框：選入“status”，擊defineevents鈕，在singlevalue框填入“1”因子factor框：選入“group”，定義最小值“1”，最大值“2”。單擊選項(xiàng)option按鈕，彈出對(duì)話框：1）√壽命表，系統(tǒng)默認(rèn)。2）圖：√生存函數(shù)3）比較第一個(gè)因子的水平：√整體比較三、主要輸出結(jié)果壽命表Intervalstarttime:生存時(shí)間區(qū)間的下限NumberEnteringInterval：該組段的觀察例數(shù)NumberWithdrawingduringInterval：該組段的刪失例數(shù)NumberExposedtoRisk：有效（校正）觀察例數(shù)NumberofTerminalEvents：陽(yáng)性事件（死亡）例數(shù)ProportionTerminating：死亡概率

ProportionSurviving：生存概率=1-死亡概率CumulativeProportionSurvivingatEndofInterval：至本組段上限的累積生存率Std.ErrorofCumulativeProportionSurvivingatEndofInterval：累積生存率的標(biāo)準(zhǔn)誤

ProbabilityDensity：概率密度Std.ErrorofProbabilityDensity：概率密度標(biāo)準(zhǔn)誤HazardRate：風(fēng)險(xiǎn)率Std.ErrorofHazardRate ：風(fēng)險(xiǎn)率的標(biāo)準(zhǔn)誤三、主要輸出結(jié)果10月生存率的估計(jì)：

甲法48%，標(biāo)準(zhǔn)誤0.1

乙法30%，標(biāo)準(zhǔn)誤0.1

兩組的中位生存期估計(jì)：3.繪制生存曲線：4.兩組生存時(shí)間分布的比較：結(jié)論：（1）從中位生存時(shí)間來看，甲組（M=10.75）大于乙組（M=8）；（2）從生存曲線分布看，甲組療效優(yōu)于乙組；（3）從兩組生存時(shí)間分布的整體比較結(jié)果看，兩組生存時(shí)間分布差異有統(tǒng)計(jì)學(xué)意義，甲組優(yōu)于乙組。Kaplan-Meier過程Kaplan-Meier過程用于(尤其小樣本資料）：估計(jì)各生存時(shí)間的生存率以及中位生存時(shí)間。繪制各種曲線：如生存函數(shù)、風(fēng)險(xiǎn)函數(shù)曲線等。比較某研究因素不同水平的生存時(shí)間有無差異?？刂颇硞€(gè)分層因素后對(duì)研究因素不同水平的生存時(shí)間分布進(jìn)行比較。對(duì)多組生存時(shí)間分布進(jìn)行兩兩比較。（各總體分布比較采用Log-rank等非參數(shù)方法）CompanyLogo實(shí)例分析例2：（數(shù)據(jù)同例1）為了比較不同手術(shù)方法治療腎上腺腫瘤的療效，某研究者隨機(jī)將43例病人分成兩組，甲組23例、乙組20例的生存時(shí)間（月）如下所示：其中有“+”者是刪失數(shù)據(jù)，表示病人仍生存或失訪，括號(hào)內(nèi)為死亡人數(shù)。（1）計(jì)算甲、乙兩法各生存時(shí)間的生存率和標(biāo)準(zhǔn)誤。（2）估計(jì)兩組的中位生存期。（3）繪制各組生存函數(shù)曲線。（4）比較兩組的總體生存時(shí)間分布有無差別。一、建立數(shù)據(jù)文件（同前）二、操作過程主菜單：分析Analyze

生存SurvivalKaplan-Meier對(duì)話框參數(shù)設(shè)置：時(shí)間time框：選入“t”。狀態(tài)status框：選入“status”，擊defineevents鈕，在singlevalue框填入“1”。因子factor框：選入“group”。單擊選項(xiàng)option按鈕，彈出對(duì)話框：1）統(tǒng)計(jì)量：√生存分析表（系統(tǒng)默認(rèn)）；√均值和中位生存時(shí)間

（系統(tǒng)默認(rèn)）2）圖：√生存函數(shù)2）水平間的兩兩比較?！?.單擊比較因子CompareFactor按鈕，彈出對(duì)話框：1）檢驗(yàn)統(tǒng)計(jì)量TestStatistics：

用于檢驗(yàn)時(shí)間分布是否相同

√對(duì)數(shù)秩Log-rank：各時(shí)間點(diǎn)的權(quán)重一樣。Breslow：按各時(shí)間點(diǎn)的觀察例數(shù)賦權(quán)。Tarone-Ware：按各時(shí)間點(diǎn)觀察例數(shù)的平方根賦權(quán)。分層后，各層作組間兩兩比較不分層，組間兩兩比較不分層，作各組間總的差異比較，默認(rèn)項(xiàng)6.單擊Save按鈕，彈出保存新變量Savenewvariables對(duì)話框：√√三、主要輸出結(jié)果生存表：略

兩組的中位生存期估計(jì)：3.繪制生存曲線：4.兩組生存時(shí)間分布的比較：結(jié)論：（1）從中位生存時(shí)間來看，甲組（M=24.23）大于乙組（M=7.8）；（2）從生存曲線分布看，甲組療效優(yōu)于乙組；（3）從兩組生存時(shí)間分布的整體比較結(jié)果看，兩組生存時(shí)間分布差異有統(tǒng)計(jì)學(xué)意義，甲組優(yōu)于乙組。Cox回歸過程COX回歸過程用于研究各種因素（稱為協(xié)變量或伴隨變量等）對(duì)于生存期長(zhǎng)短的關(guān)系，進(jìn)行多因素分析。

h(t,x)=h0(t)exp(1x1+2x2+??????+mxm

） X1，X2，????，Xm是協(xié)變量

1，2，??????，m是偏回歸系數(shù)，由樣本估計(jì)而得。

i>0表示該協(xié)變量是危險(xiǎn)因素，越大使生存時(shí)間越短

i<0表示該協(xié)變量是保護(hù)因素，越大使生存時(shí)間越長(zhǎng)實(shí)例分析例3：為探討某惡性腫瘤的預(yù)后，某研究者收集了63例患者的生存時(shí)間、生存結(jié)局及影響因素。影響因素包括病人年齡、性別、組織學(xué)類型、治療方式、淋巴結(jié)轉(zhuǎn)移、腫瘤浸潤(rùn)程度，生存時(shí)間以月計(jì)算。變量的賦值和所收集的資料分別見表17-8和表17-9。試用Cox回歸模型進(jìn)行分析。NoX1X2X3X4X5X6tYNoX1X2X3X4X5X6tY1540011052033620001012002570100051034401110140135800011351355010010261443111101030363311000120054801000713757111001200640010006003848100101200744010005803928000103183600011291405410110120193911101700413501011711042010016704247000101811142010006604349101101200124210110870444301000120013511110085045481100015114550100182046440001041154911101760476011100120016521110174048400001016117481110063049320100124118541011110105044000111911938010001000514810010120020401110166152720101024121380001093053420001021221900010241546310110120023671011093055550110012124370011090056390001051254310010151574400010120026490001031584211100120027501111187059740001171285311100120060610101040129321110012006145101101080304601001120062380100024131431011012006362000101613244101101200表17-963名某惡性腫瘤患者的生存時(shí)間（月）及影響因素一、建立數(shù)據(jù)文件（data-03.sav）二、操作過程主菜單：分析Analyze

生存SurvivalCoxRegressionCox回歸主對(duì)話框二、操作過程

主對(duì)話框參數(shù)設(shè)置：時(shí)間time框：選入“t”。狀態(tài)status框：選入“y”，擊defineevents鈕，在singlevalue框填入“1”協(xié)變量Covariates框：選入“x1~x6”。方法Method框：選擇自變量進(jìn)入Cox模型的方法，SPSS提供以下7種方法?！獭?.分析例數(shù)描述三、主要輸出結(jié)果2-1.模型檢驗(yàn)（全變量模型）結(jié)果提示：（1）對(duì)模型總體檢驗(yàn)有顯著意義（P=0.003），即至少有1個(gè)自變量的總體回歸系數(shù)不為0。2-2.模型檢驗(yàn)（逐步回歸模型，Method=向前法LR，自變量進(jìn)入P≤0.05，剔除P>0.10）（2）采用逐步回歸法進(jìn)行Cox模型分析的結(jié)果提示：模型擬合自變量進(jìn)入和剔除的檢驗(yàn)水準(zhǔn)分別為0.05和0.1時(shí)，篩選后的最佳模型包含兩個(gè)協(xié)變量，即X4（治療方式）和X5（淋巴結(jié)是否轉(zhuǎn)移），該擬合模型總體檢驗(yàn)提示具有統(tǒng)計(jì)學(xué)意義（整體卡方=17.594，P<0.001）。3.參數(shù)估計(jì)（逐步回歸模型，Method=向前法LR，自變量進(jìn)入P≤0.05，剔除P>0.10）B：偏回歸系數(shù)，SE：偏回歸系數(shù)的標(biāo)準(zhǔn)誤Wald：用于檢驗(yàn)總體偏回歸系數(shù)與0有無顯著差異，v=1時(shí)，W=（B/SE）2。Exp(B)：相對(duì)危險(xiǎn)度估計(jì)（RR值）（3）X4（治療方式）對(duì)生存時(shí)間有影響，采用新療法病人的死亡風(fēng)險(xiǎn)降至傳統(tǒng)療法的17.2%（RR的95%CI為0.059~0.503）。（4）X5（淋巴結(jié)是否轉(zhuǎn)移）對(duì)生存時(shí)間也有影響，有淋巴結(jié)轉(zhuǎn)移病人的死亡風(fēng)險(xiǎn)為無淋巴結(jié)轉(zhuǎn)移病人的2.538倍（RR的95%CI為1.062~6.066）。4.自/協(xié)變量的均值（選入2個(gè)變量的模型）5.生存曲線（選入2個(gè)變量的模型）練習(xí)1：某臨床試驗(yàn)對(duì)20名第Ⅲ或第Ⅳ期黑色素瘤的患者進(jìn)行隨訪研究，截至研究期結(jié)束，記錄的生存資料見表1(ex-01.sav)。應(yīng)用SPSS軟件建立數(shù)據(jù)文件，并計(jì)算100周生存率。表120名第Ⅲ/Ⅳ期黑色素瘤患者的生存時(shí)間（周）12.815.624.0+26.429.230