我正在遵循3.5.1中survminer包中的肺部示例。初始代码如下:
library("survminer")
require("survival")
fit <- survfit(Surv(time, status) ~ sex, data = lung)
ggsurvplot(fit, data = lung)这就是按性别生存的情况,一切都很好。然后我突然想到,我只想用图形表示sex=1或sex=2。如何从fit中仅绘制数据的子集?
发布于 2021-03-21 03:42:33
Parfait的评论给出了答案,但总的来说:
# to get all the subtypes
ggsurvplot(survfit(Surv(time, status) ~ sex, data = lung)
#to select one subtype
ggsurvplot(survfit(Surv(time, status) ~ sex, data = subset(lung, sex==1))
#to select multiple subtypes
ggsurvplot(survfit(Surv(time, status) ~ sex, data = subset(lung, sex %in% c(1,2)))https://stackoverflow.com/questions/52747497
复制相似问题