### for循环批量绘图
for (i in vars){
splots <- list()
km_fit <- survfit(Surv(Survival_months,Vital_Status)~mydata[[i]], data=mydata)
splots[[1]]<-ggsurvplot(km_fit,
xlab = "Time,mo",
ylab="Proportion Alive",
pval = T,
conf.int = F,##置信带
risk.table = T,
legend.title = i,
legend.labs = levels(mydata[[i]]),##
#surv.median.line = "hv",# 中位生存
palette="lancet")
## width=6.95,height=6.5
res<-arrange_ggsurvplots(splots, print = F,
ncol = 1, nrow = 1, risk.table.height = 0.25)
ggsave(paste(i,"All_surv.pdf",sep = "_"), res,width=7,height = 6)
# Arrange multiple ggsurvplots and print the output
}
网友评论