當包含時間相關的協變量時使用SAS繪製危害

Question

我已經使用proc Phreg和編碼過程方法在SAS中構建了一個具有時間相關協變量的Cox比例風險模型。

我有興趣繪製估計的危險率，但與時間有關的協變量似乎不支持我可以找到的圖形選項。我的結果變量被編碼爲（開始，停止）*死亡（0）。我以前使用Proc lifetest來測試危險圖，但是當以這種方式表示時間變量時無法運行它。

我將非常感謝任何建議，謝謝

Answer 1

出於某種原因，PHREG只會繪製累積風險。這種方法將爲您提供基線危險函數的圖表。

首先使用基線聲明PHREG讓倖存者函數估計，像這樣：

proc phreg data = yourdata; 
    model (start, end)*death(0) = ...; 
    baseline out = bl survival = s .../method = pl; 
run; 

再喂輸出數據集中到一個名爲smooth宏，由保羅·阿利森，使用SAS 生存分析的作者寫的：實用指南。這個宏需要PROC IML，它是從基本SAS單獨授權的，所以如果你沒有PROC IML的訪問權限，那麼你可以看看他的代碼，並使用一系列數據步驟或者你有什麼改編。

該宏可以找到here。您可能還需要考慮重寫繪圖部分以使用SGPLOT等ODS圖形程序，但這只是個人偏好。

查看this related and informative UCLA page瞭解更多關於使用宏的信息。 PHREG中BASELINE statement的SAS文檔也很有用。