0

我有一个 Cox 比例风险模型，使用 R 中的以下代码来预测死亡率。添加协变量 A、B 和 C 只是为了避免混淆（即年龄、性别、种族），但我们对预测变量 X 非常感兴趣。其中 X 是一种降低死亡风险的药物（我们的假设）。

所以我想制作一个cox生存图来显示两组（X = 0和X = 1）之间的比较，并调整其他协变量（A，B，C）。

我已经开始了第一步

cox.model <- coxph(Surv(time, dead) ~ A + B + C + X, data = df)

然后我用这段代码制作了考克斯图

ggsurvplot(survfit(cox.model), data = df, palette = "#2E9FDF",
       ggtheme = theme_minimal(),
       risk.table = TRUE)

它给了我一个只有一条（累积）线的图，而不是（X = 1和X = 0）之间的比较，甚至风险表也只是一组。

我的问题是如何在该图中显示比较？还如何在风险表中显示两组？

注意：当我做 KM 时，我很容易在不调整协变量 A+ B+ C 的情况下进行比较。

4

0 回答 0