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

Question

我想将风险表和95%的置信区间添加到调整后的生存曲线中。类似于下面的生存曲线(Link.我看到了一些相关的代码here，但它没有提到我需要什么。

另一个问题是，如果我在调整中使用获得的权重变量(从iptw package获得)，是否正确(因为我在实际数据集中有一个名为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)

​

​

Answer 1

简而言之:目前还没有R-Package可以让你直接用置信区间绘制混杂调整后的生存曲线。但是，在riskRegression R-Package中有一个名为ate的函数，可用于以95%的置信区间计算某些时间点的调整后的生存概率估计。如果你做了一些编码，你可以明确地使用它来获得你想要的东西。

或者，您可以等待更长时间。我目前正在开发一个R-Package，它实现了混淆器调整的生存曲线(以及存在竞争风险时的累积发病率函数)的各种方法，包括置信区间、假设测试、风险表等。