在 r 中使用 survminer package 將風險表和 95% 置信區間添加到調整后的生存曲線

Question

我想在調整后的生存曲線中添加風險表和 95% 置信區間。 類似於下面的生存曲線（鏈接）。 我在這里看到了一些相關的代碼，但沒有提到我需要什么。

另一個問題，如果我在調整中使用獲得的權重變量（來自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 允許您直接使用 plot 混雜因素調整生存曲線和置信區間。 然而，在riskRegression R-Package 中有一個名為ate的 function 可用於在某些時間點以 95% 的置信區間計算調整后的生存概率估計值。 如果你做了一些編碼，你肯定可以使用那個來得到你想要的。

或者，您可以等待更長時間。 我目前正在研究一個 R-Package，它實現了各種方法的混雜因素調整生存曲線（以及存在競爭風險時的累積發生率函數），並帶有置信區間、假設檢驗、風險表等。