时间戳

首页 > 解决方案 > 使用 r 中的 survminer 包将风险表和 95% 置信区间添加到调整后的生存曲线

r - 使用 r 中的 survminer 包将风险表和 95% 置信区间添加到调整后的生存曲线

问题描述

我想在调整后的生存曲线中添加风险表和 95% 置信区间。类似于下面的生存曲线（链接）。我在这里看到了一些相关的代码，但没有提到我需要什么。

另一个问题，如果我在调整中使用获得的权重变量（来自iptw 包）是否正确（因为我在实际数据集中有一个名为 weight 的列。

这是我的代码：

library(survival);library(survminer)
data(lung);names(lung)
#fit <- coxph( Surv(time, status==2) ~ ph.karno + strata(sex), data = lung )
lung$sex <- ifelse(lung$sex == 1,  "Male", "Female")
fit <- coxph(Surv(time, status) ~  ph.ecog + age +strata(sex),  data = lung)

ggadjustedcurves(fit,
                 variable = "sex", 
                 data = lung,
                 method = "average",
                 palette = c("#E69F00", "#56B4E9"),
                 size = 1.3,
                 legend = "right",
                 legend.title = expression(bold("Legend title")),
                 xlab = "Time",font.legend = 12) + theme(legend.text.align = 0.5)

在此处输入图像描述

标签： rsurvival-analysisadjustment

解决方案

简短的回答：目前没有 R-Package 允许您直接绘制带有置信区间的混杂因素调整的生存曲线。然而ate，在R-Package 中调用了一个函数riskRegression，可用于在某些时间点以 95% 的置信区间计算调整后的生存概率估计值。如果你做了一些编码，你肯定可以使用那个来得到你想要的。

或者，您可以等待更长时间。我目前正在研究一个 R-Package，它实现了各种方法的混杂因素调整生存曲线（以及存在竞争风险时的累积发生率函数），并带有置信区间、假设检验、风险表等。

推荐阅读