R语言 不同治疗方式生存率的比较 ggplot2绘图

 使用数据集1stsurg进行生存分析，其中变量X1surg包含了5种不同的治疗方式

> data<-read.csv("1stsurg.csv",header=T)
> firstsurgtreat<-data
> str(firstsurgtreat)
'data.frame':	365 obs. of  3 variables:
 $ time  : num  24.3 24.3 24.8 25.3 25.9 ...
 $ status: int  0 0 1 0 0 0 0 0 0 0 ...
 $ X1surg: num  1 2 4 4 4 1 2 2 2 4 ...

把变量X1surg转换成因子，对各个水平加上标签

> firstsurgtreat$X1surg<-factor(firstsurgtreat$X1surg,levels=c(1,2,3,4,5),labels=c("Lobectomy in the same lobe","Sublobar resections in the same lobe","Lobectomy combined with sublobar resections of different lobes ipisilaterally","Multiple sublobar resections of different lobes ipisilaterally","Bilateral surgery"))

完成数据准备后，生存分析第一步是用Surv()函数创建生存对象

> surv.obj <- Surv(time = firstsurgtreat$time, event = firstsurgtreat$status)
> surv.obj
  [1]  24.26+  24.33+  24.79   25.28+  25.87+  25.97+  26.33+  26.63+  26.76+
 [10]  27.02+  27.06+  27.48+  27.98+  28.04+  28.24+  28.27+  28.64+  29.10+
 [19]  29.19+  29.59+  30.05+  30.12+  30.25+  30.25+  30.35+  30.58+  30.74+
......

得到不同治疗方式下生存率的估计

> surv.treat<-survfit(surv.obj~X1surg,data=firstsurgtreat)
> summary(surv.treat)
Call: survfit(formula = sur