生存分析图

有时候需要将得到的临床数据画出生存曲线进行生存分析。则需要用到以下代码

rm(list = ls())
options(stringsAsFactors = F)##将字符串的列辨认成factor
a=read.table('eg.csv',sep=',',header = T,fill=T)##表中有空白
dat=a
library(ggplot2)
library(survival)
library(survminer)
table(dat$Status)
dat$Status=ifelse(dat$Status=='Dead',1,0)##进行分组
sfit=surv_fit(Surv(Days,Status)~Group,data=dat)
sfit
summary(sfit)
ggsurvplot(sfit,conf.int = F,pval = T)
ggsave('survial1')