如何提取 ggcuminc 图中特定百分位数的时间？

Question

我正在绘制累积发生函数 (CIF)，我希望能够看到特定百分位数的时间。

这是我的代码。我正在使用 tidycmprsk 包，然后添加 ggcuming 来绘制 CIF：

cif_RFS_overall <- tidycmprsk::cuminc(Surv(time, status) ~ 1, data = df) %>%  
  ggcuminc(outcome = "relapse") +
  labs(
    x = "Months"
  ) + 
  add_risktable(risktable_stats = "n.risk") +  
  scale_ggsurvfit() 

例如，我如何获得百分位 50% 或 10% 的时间？

Answer 1

您可以使用

add_quantile()

函数提供绘图上注释的分位数信息。您可以指示将在其中绘制线段的数值。在这种情况下，对于 10% 的累积发生率，您可以提供

y_value = .1

作为参数。

这是使用试验数据的完整示例。该数据集包括 200 名接受药物 A 或药物 B 的患者以及肿瘤对治疗的反应结果。

library(tidycmprsk)
library(ggsurvfit)

cuminc(Surv(ttdeath, death_cr) ~ trt, trial) %>%
  ggcuminc(outcome = "death from cancer") +
  add_risktable(risktable_stats = "n.risk") +
  add_quantile(y_value = .1) +
  scale_ggsurvfit()

剧情