r - 生存剧情传说

Question

嗨，我对 R 完全陌生。这是我第一次尝试。我正在制作一个按年龄分割的生存图。我不知道如何为每个年龄线指定颜色并将其放入图例中。任何人都可以帮忙吗？

require(survival)  # not loaded by default although installed by default
group <- age
kmsurvival1 <- survfit(Surv(as.numeric(time),event) ~ group)
plot(kmsurvival1, xlab="Time",ylab="Survival Probability", mark.time = F)

Answer 1

您只需要指定一个长度与绘图中的行数(即组数)相同的颜色向量。你可以这样做

N <- length(unique(group))
plot(kmsurvival1, xlab="Time",ylab="Survival Probability", mark.time = F,
col=1:N)
legend(
  "topright",
  legend=unique(group),
  col=1:N,
  horiz=FALSE,
  bty='n')

或者您可以手动指定颜色 col=c('black','blue','red') (取决于您需要多少种颜色)。

来自 ?plot.survfit 中的示例,

library(survival)
leukemia.surv <- survfit(Surv(time, status) ~ x, data = aml)
plot(leukemia.surv, lty = 2:3,col=3:4)
lLab <- gsub("x=","",names(leukemia.surv$strata))  ## legend labels
legend(
  "top",
  legend=lLab,
  col=3:4,
  lty=2:3,
  horiz=FALSE,
  bty='n')