R：gtsummary 包中的 tbl_survfit() 是否不再适用于多种因素的组合？

Question

我对堆栈溢出完全陌生，所以我希望我做的一切都是正确的。如果没有，请友善:)

我在使用 gtsummary 包（顺便说一句，这是一个非常棒的包）中的 tbl_survfit() 函数时遇到了问题。当我在 survfit 对象中仅使用一个因素时，tbl_survfit() 工作得非常好，即

（我将以生存包中的数据集“Trial”为例）

library(gtsummary)
library(survival)
tbl_survfit(
  survfit(Surv(ttdeath, death) ~ trt, trial),
  times = c(12, 24))

或

tbl_survfit(
  survfit(Surv(ttdeath, death) ~ stage, trial),
  times = c(12, 24))

工作正常。然而，当我使用多种因素的组合时，例如，

tbl_survfit(
  survfit(Surv(ttdeath, death) ~ trt + stage, trial),
  times = c(12, 24))

我收到一个我不明白的错误。问题似乎出在 tbl_survfit() 上，因为 survfit 对象本身

survfit(Surv(ttdeath, death) ~ trt + stage, trial)

创建时没有任何问题。我很确定我之前已经将 tbl_survfit() 与多种因素结合使用，没有任何问题，但我不知道为什么它不再工作了。我已将 R、survival 和 gtsummary 更新到最新版本，并尝试了这里的所有内容：

非常感谢您的帮助，提前谢谢您！

Answer 1

答案在您提供的“tbl_survfit 中的常见错误来源”链接中给出。您想要为模型选择列，然后将其传递到

tbl_survfit

trial %>%
  select(ttdeath, death, trt, stage) %>%
  tbl_survfit(y = Surv(ttdeath, death), times = c(12, 24))

此处与页面上的不同之处在于在

select()

语句中添加了“stage”。您需要在模型中包含您想要的任何变量。