竞争风险Cox比例风险模式.PDF

竸爭風險 Cox 比例風險模式 (Cox proportional hazard model with competing risks) 生統中心 李采娟 教授 存活資料中每個受試者在追蹤期間可經歷許多不同事件，當受試者只可能經 歷這許多事件中的其中之一時，即這些事件必須完全互斥 (mutually exclusive) ， 這些其他事件的機率就被稱為競爭風險 (competing risks) 。在醫學研究的實務應用， 可能存在有二種以上不同的事件，而感興趣的事件爲其中一種，如研究者感興趣 白血病病患進行骨髓移植後白血病復發 (relapse)的影響因素，然而這樣病患可能 發生的事件 (event)有二種：白血病復發 (relapse)和移植失敗造成的死亡 (death) ， 受試者不可能同時發生這兩個事件，在這種情形下，出現死亡的可能性爲竸爭風 險，移植失敗造成的死亡 (death)即爲白血病復發 (relapse)的竸爭事件 (competing event) 。什麼時候應考量競爭風險呢?一是當竸爭事件發生率高於感興趣事件發生 率；另一則是當追蹤期間高於五年。 本文章將介紹三種竸爭風險 Cox 比例風險迴歸模式，以糖尿病病人論質計 酬資料為範例，說明如何以 SAS程式進行這三種竸爭風險 Cox 比例風險迴歸模 式。 研究問題和資料範例 範例 : 糖尿病病人論質計酬資料 糖尿病論質計酬研究資料庫包含了46,289個病人，這些病人平均被追蹤了八 年，此研究想探討病患進入論質計酬時的糖化血色素值(glycated hemoglobin, HbA1c)和慢性阻塞性肺病 (chronic obstructive pulmonary disease, COPD)發生間的 相關，在探討此相關時，死亡和慢性阻塞性肺病為互斥事件，因此符合使用竸爭 風險 Cox 比例風險迴歸模式的條件，且資料追蹤超過五年，因此應考量死亡競爭 風險。 以下呈現4個糖尿病論質計酬病人原始數據，包含個人編號、性別、年齡、 糖化血色素值、糖尿病史、慢性阻塞性肺病的發生和死亡。 表1. 四個個案的原始資料 糖化血色 糖尿病史 慢性阻塞性 死亡 追蹤時間 編號 性別 年齡 素值 (年) 肺病的發生 (年) 101 0 70 7.3 4 0 0 7.6 102 1 75 10.9 8 1 0 5.6 103 1 42 8.5 1 0 0 7.5 104 1 52 7.9 7 0 1 1.9 1    方法一:特定因素風險函數(cause specific hazard function) 此方法使用 Cox 比例風險模式分別估計每種事件的風險和危險因子的風險 比，除了原本因失去追蹤或退出的設限(censor) ，其他競爭事件也視為設限。此 方法假設競爭風險為獨立，若此獨立假設不成立，目前沒有方法可在此假設不成 立下考量競爭風險。假設有 c種事件，則各種事件風險函數可用下列方程式表示 : ∑ , 其中 c=1, . . . . ,C ，為各種事件， 為自變項 ,