R语言 生存分析学习笔记

1.加载相应的包

library("survival")

library("survminer")

#使用生存包中提供的肺癌数据

data("lung")

head(lung)

列名：机构代码，生存时间（天），状态（1审查，2死亡），年龄，性别，ph.ecog(0-5:好-死)，ph.karno(0-100:坏-好)，pat.karno(同前，前一个是由医生评，后一个是由病人评)，膳食卡路里，减轻的体重

fit<-survfit(Surv(time,status)~sex,data=lung)#计算生存曲线
print(fit)

 

默认情况下，函数print（）显示生存曲线的简短摘要。它显示观察次数，事件数量，中位数生存和中位数的置信限。

plot(fit)#画图

可参考文章：http://blog.csdn.net/shmilyringpull/article/details/17529637