在包含时间相关的协变量时使用SAS绘制危险图

Question

我使用proc Phreg和编码处理方法在SAS中建立了Cox比例风险模型，该模型具有时间依赖的协变量。

我对绘制估计的危险率感兴趣，但我可以找到的图形选项似乎不支持时间相关的协变量。我的结果变量编码为（开始，停止）*死亡（0）。我之前使用了Proc lifetest for hazard graphs，但是当时间变量用这种方式表示时，它无法运行。

我非常感谢任何建议，谢谢

Answer 1

出于某种原因，PHREG只会绘制累积危险。这种方法将为您提供基线危险函数的图表。

首先使用PHREG中的BASELINE语句来获取幸存者函数估计值，如下所示：

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

然后将输出数据集提供给一个名为smooth的宏，由Paul Allison编写，他是使用SAS的生存分析：实用指南的作者。该宏需要PROC IML，它与基础SAS分开许可，因此如果您无法访问PROC IML，那么您可以查看他的代码并使用一系列数据步骤或您拥有的内容进行调整。

可以找到宏here。您可能还想考虑重写绘图部分以使用像SGPLOT这样的ODS图形程序，但这只是个人偏好。

查看this related and informative UCLA page以获取有关使用宏的更多信息。同样有用的是PHREG中BASELINE statement的SAS文档。