有時候需要將得到的臨床數(shù)據(jù)畫出生存曲線進行生存分析。則需要用到以下代碼

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')