【发布时间】:2021-09-02 11:20:34
【问题描述】:
我找到了这个用于竞争风险分析的图示例:Customizing a competing risks plot in R with package "cmprsk"。 Allan Cameron 的版本真的很棒,但不幸的是我无法编辑分组变量的图例标签 - 即不同的线条颜色。
# Code written by Allan Cameron in the post mentioned above (slightly edited with my label names)
library(ggplot2)
library(cmprsk)
library(survminer)
# some simulated data to get started
comp.risk.data <- data.frame("tfs.days" = rweibull(n = 100, shape = 1, scale = 1)*100,
"status.tfs" = c(sample(c(0,1,1,1,1,2), size=50, replace=T)),
"Typing" = sample(c("A","B","C","D"), size=50, replace=T))
# fitting a competing risks model
CR <- cuminc(ftime = comp.risk.data$tfs.days,
fstatus = comp.risk.data$status.tfs,
cencode = 0,
group = comp.risk.data$Typing)
p <- ggcompetingrisks(fit = CR,
multiple_panels = F,
xlab = "Days",
ylab = "Cumulative incidence of event",
title = "Competing Risks Analysis") +
scale_linetype_manual(name="my genotype",values=c(1, 2),labels=c("yes","no"))
p$mapping <- aes(x = time, y = est, colour = group, linetype = event)
p + labs(linetype = "event", colour = "my group")
我能够编辑不同线型的标签(点与否,即“是”、“否”),但我无法更改值“A”、“B”、“C”、 “D”可能是“一”、“二”、“三”、“四”。
谁能帮帮我?
谢谢!
【问题讨论】: