R语言 如何在ggcuminc plot中提取特定事件的时间?

hts6caw3  于 12个月前  发布在  其他
关注(0)|答案(1)|浏览(116)

我正在绘制累积发病率函数(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%的时间?

hjqgdpho

hjqgdpho1#

您可以使用add_quantile()函数来提供图上标注的分位数信息。您可以指示将绘制线段的数值。在这种情况下,对于10%的累积发生率,您可以提供y_value = .1作为参数。
下面是一个使用trial data的完整示例。该数据集包括200名接受药物A或药物B的患者以及肿瘤对治疗的React结果。

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()

字符串

剧情


的数据

相关问题