r - 使用 ggcompetingrisks 将审查形状添加到累积发生率图中

Question

是否可以在情节中绘制审查survminer ggcompetingrisks()？

我正在从cmprsk包装中进行竞争风险分析。该图显示了 4 种不同事件（局部/区域/远处进展或死亡）的累积发生率。如果这些事件均未发生，则使用审查。我希望在情节中显示审查，希望作为一个勾号“|” 或所有 4 行上的“+”。

我知道使用时survminer ggsurvplot()，这很容易。累积发病率图也可以吗？

# Packages
library(cmprsk)
library(survminer)

# Dummy data
set.seed(2)
failure_time <- rexp(100)
status <- factor(sample(0:4, 100, replace=TRUE), 0:4,
                 c('censor', 'A', 'B','C', 'D'))
df <- data.frame(failure_time,status)

# How many censored?
table(df$status)

# Cmprsk fit 
fit <- cuminc(ftime = df$failure_time, fstatus = df$status, cencode = "censor")

# Survminer plot
ggcompetingrisks(fit, multiple_panels = FALSE)

# I was hoping this might work:
ggcompetingrisks(fit, multiple_panels = FALSE, censor = TRUE, censor.shape = "+", censor.size = 4.5)

提前致谢

标签： rggplot2survival-analysissurvminer