如何使用R语言做生存分析?
首先,我们这里是介绍使用survival这个包来做。
软件的安装:install.packages(survival)
加载包:library(survival)
我们这里利用R自带的数据进行分析。
attach(melanom)
names(melanom)
Kaplan-Meier分析
survfit(Surv(days,status==1))
surv.all<-survfit(Surv(days,status==1))
summary(surv.all)
画生存曲线
plot(surv.all)
接着根据性别进行分组画生存曲线
surv_sex<-survfit(Surv(days,status==1))~sex)
plot(surv_sex,conf.int=T,col=c("red","green"))
接着进行对数秩检验
survfit(Surv(days,status==1)~sex)
Cox分析
summary(coxph(Surv(days,status==1)~sex))
本次分享就到这里,希望本文对大家有所帮助。
赞 (0)