是否可以在情节中绘制审查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)
提前致谢