COX模型操作流程與解析

1、Cox Regression過程例：下表資料摘自衛(wèi)生統(tǒng)計學（倪宗瓚主編）218頁的表17.6附表 裸鼠絨癌療效觀察數(shù)據(jù)編號 帶菌天數(shù) 瘤體大小 天花粉 甲藥 乙藥 結局（死=1） 生存天數(shù)No X1 X2 X3 X4 X5 d day 1 19 25 0 0 0 1 8 2 17 16 0 0 0 1 9 3 19 37 0 0 0 1 8 4 16 19 0 0 0 1 8 5 14 25 1 0 0 0 18 6 13 18 1 0 0 1 17 7 16 25 1 0 0 1 14 8 9 10 1 0 0 1 15 9 9 22 0 1 0 1 1510 10 25 0 1 0 1 1

2、111 14 25 0 1 0 1 1312 12 37 0 1 0 1 1213 17 37 0 0 1 1 914 14 29 0 0 1 1 1215 13 13 0 0 1 1 1216 17 31 0 0 1 1 10一、界面說明 Cox回歸主對話框【Time】框、【Status】框前文已經介紹過了，這里我就不再廢話嘮叨的了。Block 1 of 1右邊的Next鈕被激活。這個按鈕用于確定不同自變量進入回歸方程的方法，詳見Method框的內容。用同一種方法進入回歸方程的自變量在同一個Covariates框內?！綜ovariates】框選入自/協(xié)變量，即選入你認為可能對生存時間有影響的

3、變量?！綧ethod】框選擇自變量進入Cox回歸方程的方法，SPSS提供下面幾種方法：Enter: Covariates框內的全部變量均進入回歸模型。Forward: Conditional: 基于條件參數(shù)估計的向前法。Forward: LR: 基于偏最大似然估計的向前法。Forward: Wald: 基于Wald統(tǒng)計量的向前法。Backward: Conditional: 基于條件參數(shù)估計的后退法。Backward: LR: 基于偏最大似然估計的后退法。Backward: Wald: 基于Wald統(tǒng)計量的后退法?！維trata】框定義分層因素,將生存時間按分層因素分別進行Cox回歸?！綜at

4、egorical】選項用于告訴系統(tǒng)，Covariates框內的變量中哪些是分類變量或字符型變量。系統(tǒng)默認字符型變量為分類變量，數(shù)字型變量為連續(xù)型變量。選入自變量后，categorical鈕被激活。按categorical鈕，進入確定分類變量的對話框。見圖10。圖10 確定分類變量對話框左邊的Covariates框中列出了剛剛被選取的自變量，將分類變量選入Categorical Covariates框中。此時Change Contrast框被激活，請你選擇比較方法，即計算參數(shù)OR/i的方法。當選入分類變量后，Change Contrast框被激活，此時可選擇比較方法。SPSS提供下面幾種比較方法

5、。Indicator:指示對比。用于指定某一分類變量的基線，即參照水平。這樣計算出來的參數(shù)OR/i是以該變量的第一個或最后一個水平為基準水平（取決于下面的reference category中你選擇的是last還是first）。在這里SPSS自動創(chuàng)建亞變量，對照水平在對比分類矩陣中用0行代表。Simple:差別對比?？捎嬎阍摲诸愖兞康母魉脚c參照水平相比的OR值。參照水平自己當然就不用跟自己相比了。對于本例來說，Simple與Indicator選項是一樣的，前提是下面的Reference Category中你所選擇的同是last（或first）。Difference:差別對比。分類變量欲比較水

6、平與其前面的各水平平均值進行比較，當然也不包括第一水平。與Helmert法相反，因此也叫反Helmert法。如3水平與1、2水平的平均值相比，下同。Helmert:赫爾默特對比。分類變量欲比較水平與其后面各水平平均值進行比較，當然不包括最后一個水平。Repeated:重復對比。分類變量的各水平與其前面相鄰的水平相比較（第一水平除外）。Polynomial:多項式對比。僅用于數(shù)字型的分類變量。無效假設是假設各水平是等距離的（可以是線性的關系，也可以是立方、四次方的關系）。例如年齡每增加10歲，死亡風險的增加值是一樣的，但實際情況常常與此相反，如在20歲與60歲年齡段，年齡都增加10歲，所增加的死

7、亡風險肯定是不一樣的，具體情況需根據(jù)各人的研究課題，專業(yè)而定。Deviation:離差對比。除了所規(guī)定的參照水平外，其余每個水平均與總體水平相比。Reference category：如果你選擇了Deviation, Simple, 或Indicator三個選項，就必須選擇First或Last作為參照水平。完成上述選擇后，擊change鈕，確認選擇。你若對上面寫的一段不感興趣的話，可跳過去，直接用系統(tǒng)默認的選項。【Plots】選項圖11 Cox回歸統(tǒng)計圖對話框 Survival：累積生存函數(shù)曲線。 Hazard：累積風險函數(shù)曲線。 Log minus log：對數(shù)累積生存函數(shù)乘以-1后再取對數(shù)

8、。 One minus survival：生存函數(shù)被1減后的曲線。Change Value：系統(tǒng)默認用各變量的均數(shù)進行作圖。若用分類變量的其它水平進行作圖，則選定該變量，此時Change Value鈕被激活，按Value鈕，在其右邊的框內輸入你所想要用于作圖的值。擊Change。Separate Line for：輸入分類變量的名稱，此時可以用分類變量的不同水平進行作圖，此分類變量必須包括在前面的自變量框中?！維ave】存為新變量圖12 Cox回歸存為新變量對話框l Survival：生存函數(shù)。Function：累積生存函數(shù)估計值。Standard error：累積生存函數(shù)估計值的標準誤。Lo

9、g minus log：對數(shù)累積生存函數(shù)乘以-1后再取對數(shù)。l Diagnostics：回歸診斷。Hazard function Cox-Snell：殘差。Partial residual：偏殘差。Dfbeta(s)：剔除某一觀察單位后的回歸系數(shù)變化量。X*Beta：線性預測得分?！綩ptions】選項擊Options按鈕，彈出選項對話框。圖13 Cox回歸選項對話框l Model Statistics：模型統(tǒng)計量。CI for exp() 95%：相對危險度的可信區(qū)間。系統(tǒng)默認95%可信區(qū)間。Correlation of estimates：回歸系數(shù)的相關陣。l Display model：

10、輸出模型方式。l At each step：輸出每一步的模型。系統(tǒng)默認。l At last step：輸出最后一步的模型。l Probability for Stepwise：模型保留變量的顯著性水平。Entry：系統(tǒng)默認選入變量為P0.05。Removal：系統(tǒng)默認剔除變量為P0.10。Maximum Iterations：最大迭代次數(shù)，系統(tǒng)默認20次。Display baseline function：輸出風險基準函數(shù)以及基于各協(xié)變量均值的生存函數(shù)與風險函數(shù)。COX分析操作流程與結果解析一、 建庫定義7個列變量：時間變量： 取名“day”。生存狀態(tài)變量：取名“status”，并賦值：0=“

11、截尾”，1=“死亡”。X1(取名和賦值略)X2(取名和賦值略)X3(取名和賦值略)X4(取名和賦值略)X5(取名和賦值略)二、分析操作步驟：1. Analyze=Survival =Cox regression2. Time框：選入day3. Status框：選入結局；擊define events鈕,在single value框右邊的空格中輸入1；4. Covariate框：選入x1,x2,x3,x4,x5；5 Method: 下拉列表：Enter6. Plots 列表框：Plot Type:選survival;7Option列表框： Model Statistics: 選CI for exp():輸出exp( )的95%可信區(qū)間。8 單擊OK鈕三、結果解釋：Cox Regression上表輸出總例數(shù)、刪失例數(shù)、失訪例數(shù)。Block 0: Beginning Block模型擬合迭代過程，可