使用数据集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 = surv.