醫(yī)學(xué)研究的數(shù)據(jù)管理與分析——第十一章(1)

1、第十一章 SPSS軟件應(yīng)用（四）第一節(jié)第一節(jié) 生存分析生存分析隨訪(fǎng)研究隨訪(fǎng)研究隨訪(fǎng)研究隨訪(fǎng)研究(follow-up study)(follow-up study)是醫(yī)學(xué)中常用前瞻性研是醫(yī)學(xué)中常用前瞻性研究究. .例例: :兩種方法腎移植病人術(shù)后腎的生存時(shí)間和結(jié)局兩種方法腎移植病人術(shù)后腎的生存時(shí)間和結(jié)局( (生生存率存率) )比較比較. .例例: :不同方法對(duì)某病人不同方法對(duì)某病人( (癌癥、反復(fù)發(fā)作疾?。┥鏁r(shí)癌癥、反復(fù)發(fā)作疾病）生存時(shí)間間( (緩解時(shí)間緩解時(shí)間) )與結(jié)局與結(jié)局( (生存率生存率) )比較比較. .該類(lèi)數(shù)據(jù)通過(guò)隨訪(fǎng)得到，稱(chēng)為隨訪(fǎng)資料。該類(lèi)數(shù)據(jù)通過(guò)隨訪(fǎng)得到，稱(chēng)為隨訪(fǎng)資料。隨訪(fǎng)研

2、究隨訪(fǎng)研究(follow-up study)(follow-up study)示意圖示意圖合格的合格的研究對(duì)象研究對(duì)象 出現(xiàn)結(jié)果出現(xiàn)結(jié)果 尚未出現(xiàn)結(jié)果尚未出現(xiàn)結(jié)果 失訪(fǎng)、脫落失訪(fǎng)、脫落試驗(yàn)組試驗(yàn)組對(duì)照組對(duì)照組伴隨因素干擾因素生存分析生存分析在臨床醫(yī)學(xué)中, 對(duì)病人療效考查：1.治療結(jié)局？2.生存時(shí)間？ 生存分析生存分析（survival analysis）是將是將事件的結(jié)果事件的結(jié)果（終（終點(diǎn)事件）和出現(xiàn)結(jié)果經(jīng)歷的點(diǎn)事件）和出現(xiàn)結(jié)果經(jīng)歷的時(shí)間時(shí)間結(jié)合起來(lái)分析的一種統(tǒng)結(jié)合起來(lái)分析的一種統(tǒng)計(jì)分析方法。計(jì)分析方法。一、隨訪(fǎng)數(shù)據(jù)概念一、隨訪(fǎng)數(shù)據(jù)概念 1.分析的變量（分析的變量（y） 1） 結(jié)局事件：指結(jié)

3、局出現(xiàn)的特征，如疾病的死結(jié)局事件：指結(jié)局出現(xiàn)的特征，如疾病的死亡、復(fù)發(fā)、發(fā)生亡、復(fù)發(fā)、發(fā)生( y=1或或0) 。 2）時(shí)間間隔變量）時(shí)間間隔變量 記為（記為（t） 生存時(shí)間生存時(shí)間t=結(jié)局事件出現(xiàn)日期結(jié)局事件出現(xiàn)日期 事件的起事件的起始日期始日期 （起始日期可規(guī)定（起始日期可規(guī)定:如診斷、用藥、手術(shù)日期如診斷、用藥、手術(shù)日期等）等）, t的單位：可用年、月、周表示的單位：可用年、月、周表示生存分析的基本概念生存分析的基本概念 起始事件 隨訪(fǎng)時(shí)間 終點(diǎn)事件疾病確診 死亡疾病確診 痊愈治療開(kāi)始 死亡治療開(kāi)始 痊愈癥狀緩解 疾病惡化接觸毒物 出現(xiàn)毒性反映接觸危險(xiǎn)因素 發(fā)病2.2.截尾數(shù)據(jù)截尾數(shù)據(jù) 觀(guān)

4、察過(guò)程中個(gè)體因其他原因未觀(guān)察到明確觀(guān)察過(guò)程中個(gè)體因其他原因未觀(guān)察到明確的結(jié)局，的結(jié)局， 稱(chēng)為截尾或刪失數(shù)據(jù)（稱(chēng)為截尾或刪失數(shù)據(jù)（ censored censored datadata）。）。截尾原因有：截尾原因有：失訪(fǎng)，退出研究，如其他原因死亡。研究時(shí)間結(jié)失訪(fǎng)，退出研究，如其他原因死亡。研究時(shí)間結(jié)束，未出現(xiàn)結(jié)局事件。束，未出現(xiàn)結(jié)局事件。截尾值（截尾值（censored valuecensored value）：）：時(shí)間（時(shí)間（t t）= =截尾事件日期截尾事件日期 起始事件日期起始事件日期記為記為t t+ +。（例：。（例：1010+ +月）月）3.3.生存數(shù)據(jù)的特點(diǎn)生存數(shù)據(jù)的特點(diǎn)1 1）完全

5、數(shù)據(jù)：研究對(duì)象在規(guī)定研究期間提供確切的）完全數(shù)據(jù)：研究對(duì)象在規(guī)定研究期間提供確切的“時(shí)間時(shí)間和結(jié)局和結(jié)局”。2 2）截尾數(shù)據(jù)）截尾數(shù)據(jù)(t(t+ +) ) ：截尾數(shù)據(jù)雖然提供的信息不完全，但：截尾數(shù)據(jù)雖然提供的信息不完全，但提供了部分信息，如提供了部分信息，如 t=10t=10+ +年年9 9年。年。3 3）生存數(shù)據(jù)的結(jié)果變量（）生存數(shù)據(jù)的結(jié)果變量（Y Y ）有兩個(gè)：）有兩個(gè)： 時(shí)間（時(shí)間（t t）值）值 ，t t0 0 結(jié)局狀態(tài)（結(jié)局狀態(tài)（y y ）=“ =“ 如死亡或截尾值如死亡或截尾值” 臨床隨訪(fǎng)研究中的完全數(shù)據(jù)與截尾數(shù)據(jù)臨床隨訪(fǎng)研究中的完全數(shù)據(jù)與截尾數(shù)據(jù)觀(guān)察起點(diǎn)觀(guān)察起點(diǎn)觀(guān)察終點(diǎn)觀(guān)察終點(diǎn)

6、5 6123 4 起點(diǎn)事件 終點(diǎn)事件 失訪(fǎng)生存時(shí)間資料的特點(diǎn)生存時(shí)間資料的特點(diǎn) 2個(gè)效應(yīng)變量（1）生存時(shí)間(天數(shù))，（2）結(jié)局(死亡與否、是否陽(yáng)性等)截尾數(shù)據(jù)(censored data)：截尾數(shù)據(jù)提供的信息是不完全的 (uncompleted data)分布類(lèi)型復(fù)雜：生存時(shí)間分布常呈非正態(tài)分布對(duì)數(shù)正態(tài)分布指數(shù)分布， Weibull分布Gamma分布，logistic分布， 2個(gè)效應(yīng)變量（1）生存時(shí)間(天數(shù))，（2）結(jié)局(死亡與否、是否陽(yáng)性等)兩種錯(cuò)誤的做法：兩種錯(cuò)誤的做法：錯(cuò)誤1：忽略生存時(shí)間，采用Logistic回歸分析死亡率錯(cuò)誤2：忽略結(jié)局，采用t檢驗(yàn)、線(xiàn)性回歸分析生存時(shí)間兩種錯(cuò)誤的做法

7、：兩種錯(cuò)誤的做法：錯(cuò)誤1：丟棄丟棄截尾截尾數(shù)據(jù)，只考慮確切數(shù)據(jù)。（損失數(shù)據(jù)，只考慮確切數(shù)據(jù)。（損失了信息）了信息）錯(cuò)誤2：將將截尾截尾數(shù)據(jù)當(dāng)作確切數(shù)據(jù)處理。（低估了數(shù)據(jù)當(dāng)作確切數(shù)據(jù)處理。（低估了生存時(shí)間的平均水平）。生存時(shí)間的平均水平）。截尾( (刪失刪失) )數(shù)據(jù)兩種錯(cuò)誤的做法：兩種錯(cuò)誤的做法：錯(cuò)誤1：采用平均生存時(shí)間而不是采用采用平均生存時(shí)間而不是采用中位生存時(shí)間來(lái)表示生存時(shí)間的平均水平。來(lái)表示生存時(shí)間的平均水平。錯(cuò)誤2：采用常規(guī)采用常規(guī)t t檢驗(yàn)或方差分析進(jìn)行組間比較。檢驗(yàn)或方差分析進(jìn)行組間比較。（應(yīng)采（應(yīng)采用用log-ranklog-rank檢驗(yàn)檢驗(yàn)比較幾組生存時(shí)間比較幾組生存時(shí)間 ）

8、非正態(tài)數(shù)據(jù)生存時(shí)間資料整理和記錄生存時(shí)間資料整理和記錄例：某人研究手術(shù)方法治療例：某人研究手術(shù)方法治療2323位腎上腺腫瘤病人的生存時(shí)位腎上腺腫瘤病人的生存時(shí)間（月）如下：間（月）如下： 1 1，3 3，5 5（3 3），），6 6（3 3），），7 7，8 8，1010（2 2），），1414+ +，1717，1919+ +，2020+ +，2222+ +，2626+ +，3131+ +，3434，3434+ +，4444，5959注注：( )( )括號(hào)內(nèi)的數(shù)為相同時(shí)間點(diǎn)的人數(shù)括號(hào)內(nèi)的數(shù)為相同時(shí)間點(diǎn)的人數(shù)數(shù)據(jù)另一種記錄數(shù)據(jù)另一種記錄: :對(duì)象編號(hào)對(duì)象編號(hào) 生存時(shí)間生存時(shí)間(t) (t) 結(jié)局

9、狀態(tài)結(jié)局狀態(tài)(0(0為截尾值為截尾值) ) 1 1 1 1 1 1 2 14 0 2 14 0生存分析主要研究的內(nèi)容生存分析主要研究的內(nèi)容1.1.統(tǒng)計(jì)描述統(tǒng)計(jì)描述：計(jì)算不同時(shí)間點(diǎn)（：計(jì)算不同時(shí)間點(diǎn)（t t）的生存率，描述生）的生存率，描述生存過(guò)程。存過(guò)程。2.2.統(tǒng)計(jì)推斷統(tǒng)計(jì)推斷：檢驗(yàn)不同處理方式的生存過(guò)程有無(wú)統(tǒng)計(jì)：檢驗(yàn)不同處理方式的生存過(guò)程有無(wú)統(tǒng)計(jì)差別差別. .3.3.自變量（自變量（x x）對(duì)生存時(shí)間（）對(duì)生存時(shí)間（t t）的關(guān)系）的關(guān)系：影響生存時(shí)：影響生存時(shí)間的危險(xiǎn)因素分析間的危險(xiǎn)因素分析. .1 平均生存時(shí)間 ti的平均數(shù) (當(dāng)有截尾數(shù)據(jù)時(shí)，求不出)2 中位生存時(shí)間 ti的中位數(shù) :

10、 50%的病人生存, 50%的病人死亡的時(shí)間 (當(dāng)截尾較早時(shí)，求不出)1 生存函數(shù)(survival function) 生存率2 死亡率1生存率觀(guān)察病人總數(shù)的病人數(shù)生存時(shí)間tTtTPtS)()(常見(jiàn)生存函數(shù)的類(lèi)型生存曲線(xiàn)生存曲線(xiàn) 以生存時(shí)間為橫軸，生存率為縱軸繪制的曲以生存時(shí)間為橫軸，生存率為縱軸繪制的曲線(xiàn)線(xiàn) (a)研究終止在475天 (b) 研究終止在474天圖3 乘積極限法生存曲線(xiàn)(階梯形階梯形)及其半數(shù)生存期（Md=158天）0100200300400500t (day)0.00.20.40.60.81.0S(t)S(132)=0.50100200300400500t (day)0.0

11、0.20.40.60.81.0S(t)Md=158012345t (year)0.00.20.40.60.81.0S(t)Md=1.7圖 壽命表法生存曲線(xiàn)(折線(xiàn)折線(xiàn))及其半數(shù)生存期（Md=1.7年）3 風(fēng)險(xiǎn)函數(shù)(hazard function)tttttth時(shí)刻尚存的病人數(shù)在的病人數(shù)死于區(qū)間),()(常見(jiàn)風(fēng)險(xiǎn)函數(shù)的類(lèi)型4 風(fēng)險(xiǎn)比(hazard ratio) 相對(duì)危險(xiǎn)度(Relative risk)()(21thth第二組的第一組的風(fēng)險(xiǎn)比 非參數(shù)方法非參數(shù)方法 參數(shù)方法參數(shù)方法乘積乘積-極限法極限法 指數(shù)分布指數(shù)分布(又稱(chēng)又稱(chēng)Kaplan-Meier法法) Weibull分布分布?jí)勖矸椒▔勖?/p>

12、方法 Gamma分布分布Turnbull估計(jì)估計(jì) logistic分布分布 對(duì)數(shù)正態(tài)分布對(duì)數(shù)正態(tài)分布非參數(shù)方法非參數(shù)方法 參數(shù)方法參數(shù)方法時(shí)序時(shí)序(log-rank)檢驗(yàn)檢驗(yàn) 分布參數(shù)檢驗(yàn)分布參數(shù)檢驗(yàn)分層分層時(shí)序時(shí)序檢驗(yàn)檢驗(yàn)Gehan檢驗(yàn)檢驗(yàn)廣義廣義Wilcoxon檢驗(yàn)檢驗(yàn)Mantel-Haenszel檢驗(yàn)檢驗(yàn)生存分析的目的：生存分析的目的：(1)(1)描述生存過(guò)程描述生存過(guò)程：估計(jì)不同時(shí)間的總體生存率，計(jì)算中位生存期，估計(jì)不同時(shí)間的總體生存率，計(jì)算中位生存期，繪制生存函數(shù)曲線(xiàn)。統(tǒng)計(jì)方法包括繪制生存函數(shù)曲線(xiàn)。統(tǒng)計(jì)方法包括Kaplan-MeierKaplan-Meier（K-MK-M）法、壽命）

13、法、壽命表法。表法。(2)(2)比較：比較：比較不同處理組的生存率，如比較不同療法治療腦瘤的比較不同處理組的生存率，如比較不同療法治療腦瘤的生存率，以了解哪種治療方案較優(yōu)。統(tǒng)計(jì)方法生存率，以了解哪種治療方案較優(yōu)。統(tǒng)計(jì)方法log-ranklog-rank檢驗(yàn)等。檢驗(yàn)等。(3)(3)影響因素分析：影響因素分析：研究某個(gè)或某些因素對(duì)生存率或生存時(shí)間的影研究某個(gè)或某些因素對(duì)生存率或生存時(shí)間的影響作用。如為改善腦瘤病人的預(yù)后，應(yīng)了解影響病人預(yù)后的主要響作用。如為改善腦瘤病人的預(yù)后，應(yīng)了解影響病人預(yù)后的主要因素，包括病人的年齡、性別、病程、腫瘤分期、治療方案等。因素，包括病人的年齡、性別、病程、腫瘤分期、

14、治療方案等。統(tǒng)計(jì)方法統(tǒng)計(jì)方法coxcox比例風(fēng)險(xiǎn)回歸模型等。比例風(fēng)險(xiǎn)回歸模型等。(4)(4)預(yù)測(cè)：預(yù)測(cè)：建立建立coxcox回歸預(yù)測(cè)模型?；貧w預(yù)測(cè)模型。生存分析生存分析生存分析生存分析( (Survival Analysis) )菜單菜單壽命表壽命表（Life Tables）過(guò)程過(guò)程Life tables 過(guò)程用于（小樣本和大樣本資料）：1. 估計(jì)某生存時(shí)間的生存率，以及中位生存時(shí)間。2. 繪制各種曲線(xiàn)：如生存函數(shù)、風(fēng)險(xiǎn)函數(shù)曲線(xiàn)等。3. 對(duì)某一研究因素不同水平的生存時(shí)間分布的比較。4. 控制另一個(gè)因素后對(duì)研究因素不同水平的生存時(shí)間分布的比較。5. 對(duì)多組生存時(shí)間分布進(jìn)行兩兩比較。（比較總體生存

15、時(shí)間分布采用wilcoxon檢驗(yàn)）Company LogoCompany Logo實(shí)例分析實(shí)例分析例11-1-1Company LogoCompany Logo實(shí)例分析實(shí)例分析例1：為了比較不同手術(shù)方法治療腎上腺腫瘤的療效，某研究者隨機(jī)將43例病人分成兩組，甲組23例、乙組20例的生存時(shí)間（月）如下所示：其中有“+”者是刪失數(shù)據(jù)，表示病人仍生存或失訪(fǎng)，括號(hào)內(nèi)為死亡人數(shù)。（1）計(jì)算甲、乙兩法術(shù)后10月的生存率和標(biāo)準(zhǔn)誤。（2）估計(jì)兩組的中位生存期。（3）繪制各組生存函數(shù)曲線(xiàn)。（4）比較兩組的總體生存時(shí)間分布有無(wú)差別。Company LogoCompany Logo一、建立數(shù)據(jù)文件一、建立數(shù)據(jù)文件定

16、義定義5個(gè)變量：個(gè)變量：生存時(shí)間變量：t，值標(biāo)簽“生存時(shí)間（月）”生存狀態(tài)變量 ：status，取值“1=死亡，0=刪失或存活”頻數(shù)變量：freq，值標(biāo)簽“人數(shù)”分組變量：group，取值“1=甲組，2=乙組”生存時(shí)間序號(hào)變量(可無(wú)）：i二、操作過(guò)程二、操作過(guò)程主菜單：分析主菜單：分析Analyze生存生存Survival壽命表壽命表Life tables對(duì)話(huà)框參數(shù)設(shè)置：對(duì)話(huà)框參數(shù)設(shè)置：1. 時(shí)間time框：選入 “t”。2. 顯示時(shí)間間隔Display time intervals框：步長(zhǎng)by前面填入最大生存時(shí)間的上限（必須包括生存時(shí)間最大值），步長(zhǎng)by后面填入生存時(shí)間的組距。本例上限填“6

17、0”，組距填“1”。3. 狀態(tài)status框：選入“status”，擊define events 鈕，在single value框填入“1”4. 因子factor框：選入“group”，定義最小值“1”，最大值“2”。5. 單擊選項(xiàng)option按鈕，彈出對(duì)話(huà)框： 1）壽命表，系統(tǒng)默認(rèn)。 2）圖： 生存函數(shù) 3）比較第一個(gè)因子的水平： 整體比較三、主要輸出結(jié)果三、主要輸出結(jié)果1.10月生存率的估計(jì)：月生存率的估計(jì)： 甲法甲法 48%，標(biāo)準(zhǔn)誤，標(biāo)準(zhǔn)誤 0.1 乙法乙法 30%，標(biāo)準(zhǔn)誤，標(biāo)準(zhǔn)誤 0.1 2.兩組的中位生存期估計(jì)：兩組的中位生存期估計(jì)：Company LogoCompany Logo3

18、. 繪制生存曲線(xiàn)：繪制生存曲線(xiàn)：Company LogoCompany Logo4. 兩組生存時(shí)間分布的比較：兩組生存時(shí)間分布的比較：Kaplan-Meier 過(guò)程過(guò)程Kaplan-Meier過(guò)程用于(尤其小樣本資料）：1. 估計(jì)各生存時(shí)間的生存率以及中位生存時(shí)間。2. 繪制各種曲線(xiàn)：如生存函數(shù)、風(fēng)險(xiǎn)函數(shù)曲線(xiàn)等。3. 比較某研究因素不同水平的生存時(shí)間有無(wú)差異。4. 控制某個(gè)分層因素后對(duì)研究因素不同水平的生存時(shí)間分布進(jìn)行比較。5. 對(duì)多組生存時(shí)間分布進(jìn)行兩兩比較。（各總體分布比較采用Log-rank等非參數(shù)方法）Company LogoCompany Logo實(shí)例分析實(shí)例分析例11-1-2Com

19、pany LogoCompany Logo實(shí)例分析實(shí)例分析例2：（數(shù)據(jù)同例1）為了比較不同手術(shù)方法治療腎上腺腫瘤的療效，某研究者隨機(jī)將43例病人分成兩組，甲組23例、乙組20例的生存時(shí)間（月）如下所示：其中有“+”者是刪失數(shù)據(jù)，表示病人仍生存或失訪(fǎng)，括號(hào)內(nèi)為死亡人數(shù)。（1）計(jì)算甲、乙兩法各生存時(shí)間的生存率和標(biāo)準(zhǔn)誤。（2）估計(jì)兩組的中位生存期。（3）繪制各組生存函數(shù)曲線(xiàn)。（4）比較兩組的總體生存時(shí)間分布有無(wú)差別。一、建立數(shù)據(jù)文件（同前）一、建立數(shù)據(jù)文件（同前）二、操作過(guò)程二、操作過(guò)程主菜單：分析主菜單：分析Analyze生存生存SurvivalKaplan-Meier對(duì)話(huà)框參數(shù)設(shè)置：對(duì)話(huà)框參數(shù)設(shè)

20、置：1. 時(shí)間time框：選入 “t”。2. 狀態(tài)status框：選入“status”，擊define events 鈕，在single value框填入“1”。3. 因子factor框：選入“group”。4. 單擊選項(xiàng)option按鈕，彈出對(duì)話(huà)框： 1）統(tǒng)計(jì)量： 生存分析表，系統(tǒng)默認(rèn)。 均值和中位生存時(shí)間，系統(tǒng)默認(rèn)。 2）圖： 生存函數(shù)5. 單擊比較因子Compare Factor按鈕，彈出對(duì)話(huà)框： 1）檢驗(yàn)統(tǒng)計(jì)量Test Statistics： 都用于檢驗(yàn)時(shí)間分布是否相同。 對(duì)數(shù)秩Log-rank：各時(shí)間點(diǎn)的權(quán)重一樣。 Breslow：按各時(shí)間點(diǎn)的觀(guān)察例數(shù)賦權(quán)。 Tarone-Ware：

21、按各時(shí)間點(diǎn)觀(guān)察例數(shù)的平方根賦權(quán)。二、操作過(guò)程二、操作過(guò)程2）水平間的兩兩比較。水平間的兩兩比較。6. 單擊Save按鈕，彈出保存新變量Save new variables 對(duì)話(huà)框：三、主要輸出結(jié)果三、主要輸出結(jié)果1.生存表：生存表： 略略 2.兩組的中位生存期估計(jì)：兩組的中位生存期估計(jì)：3. 繪制生存曲線(xiàn)：繪制生存曲線(xiàn)：Company LogoCompany Logo4. 兩組生存時(shí)間分布的比較：兩組生存時(shí)間分布的比較：Cox回歸回歸過(guò)程過(guò)程Cox回歸過(guò)程用于： 1. 多個(gè)因素對(duì)生存時(shí)間的影響作用分析和比較 2. 生存（或死亡）風(fēng)險(xiǎn)預(yù)測(cè)Cox回歸模型：h(t/X)=h0(t) exp (1 X

22、1 + 2 X2 + + p Xp )h0(t)： 基準(zhǔn)風(fēng)險(xiǎn)函數(shù) 即所有變量取零時(shí)的t時(shí)刻的風(fēng)險(xiǎn)函數(shù)X1、X2 Xp：影響因素 變量1、2 p：回歸系數(shù)實(shí)例分析實(shí)例分析例11-1-3實(shí)例分析實(shí)例分析例3：為探討某惡性腫瘤的預(yù)后，某研究者收集了63例患者的生存時(shí)間、生存結(jié)局及影響因素。影響因素包括病人年齡、性別、組織學(xué)類(lèi)型、治療方式、淋巴結(jié)轉(zhuǎn)移、腫瘤浸潤(rùn)程度，生存時(shí)間以月計(jì)算。變量的賦值和所收集的資料分別見(jiàn)下表。試用Cox回歸模型進(jìn)行分析。NoX1X2X3X4X5X6tYNoX1X2X3X4X5X6tY15400110520336200010120025701000510344011101 40

23、135800011351355010010 26144311110103036331100012005480100071375711100120064001000600384810010120074401000580392800010 3183600011291405410110120193911101700413501011 71104201001670424700010 18111420100066043491011012001242101108704443010001200135111100850454811000 151145501001820464400010 41154911101

24、7604760111001200165211101740484000010 161174811100630493201001 2411854101111010504400011 19119380100010005148100101200204011101661527201010 241213800010930534200010 212219000102415463101101200236710110930555501100 121243700110900563900010 5125431001015157440001012002649000103158421110012002750111118

25、70597400011 712853111001200606101010 401293211100120061451011010803046010011200623801000 2413143101101200636200010 1613244101101200表表17-9 63名某惡性腫瘤患者的生存時(shí)間（月）及影響因素名某惡性腫瘤患者的生存時(shí)間（月）及影響因素一、建立數(shù)據(jù)文件一、建立數(shù)據(jù)文件二、操作過(guò)程二、操作過(guò)程主菜單：分析主菜單：分析Analyze生存生存SurvivalCox RegressionCox回歸主對(duì)話(huà)框回歸主對(duì)話(huà)框二、操作過(guò)程主對(duì)話(huà)框參數(shù)設(shè)置：主對(duì)話(huà)框參數(shù)設(shè)置：1.時(shí)間ti

26、me框：選入 “ t ”。2.狀態(tài)status框：選入“ y ”，擊define events 鈕，在single value框填入“1”3.協(xié)變量Covariates框：選入“x1x6”。4.方法Method框：選擇自變量進(jìn)入Cox模型的方法，SPSS提供以下7種方法。1.分析例數(shù)描述案例處理摘要案例處理摘要2641.3%3758.7%63100.0%0.0%0.0%0.0%0.0%63100.0%事件a刪失合計(jì)分析中可用的案例帶有缺失值的案例帶有負(fù)時(shí)間的案例層中的最早事件之前刪失的案例合計(jì)刪除的案例合計(jì)N百分比因變量: ta. 三、主要輸出結(jié)果三、主要輸出結(jié)果2-1.模型檢驗(yàn)（全變量模型）

27、模模 型型 系系 數(shù)數(shù) 的的 綜綜 合合 測(cè)測(cè) 試試a a, ,b b180.05220.1746.00321.9426.00121.9426.001-2 倍對(duì)數(shù)似然值卡方dfSig.整體 (得分)卡方dfSig.從上一步驟開(kāi)始更改卡方dfSig.從上一塊開(kāi)始更改起始?jí)K編號(hào) 0，最初的對(duì)數(shù)似然函數(shù):-2 倍對(duì)數(shù)似然值: 201.994a. 起始?jí)K編號(hào) 1. 方法 = 輸入b. 結(jié)果提示：結(jié)果提示：（1）對(duì)模型總體檢驗(yàn)有顯著意義（）對(duì)模型總體檢驗(yàn)有顯著意義（P=0.003），即至少），即至少有有1個(gè)自變量的總體回歸系數(shù)不為個(gè)自變量的總體回歸系數(shù)不為0。2-2.模型檢驗(yàn)（逐步回歸模型，Method=向前法LR，自變量進(jìn)入P0.05，剔除P0.10 ）模