我正在使用 carData::Rossi 数据来估计生存曲线。这是我用来生成绘图的代码:
library(carData); library(survival)
km1 <- survfit(Surv(week,arrest) ~ race,type = "kaplan-meier",data=Rossi)
library(ggplot2); library(ggfortify)
autoplot(km1) +
labs(x="Time", y="Proportion", strata= "Raza")
我得到了这个结果:
如何将分层标签更改为 Race?