library(survival)
Lung <- lung

Kaplan-Meier曲线

#基本流程：Surv()提取生存时间；survfit()生存模型;survdiff():多组时分析差异
#绘图：plot

Lusurv <- Surv(time=Lung$time,event=Lung$status)
#输出结局时间（+：还没有死亡）

Lufit <- survfit(Lusurv~Lung$sex)
#观察sex对生存时间的关系

绘图及美化

plot(Lufit,cof.int="none",col=c("red","blue"),lwd=2,mark.time=T) 
#cof.int 是否显示置信区间
#mark.time=T 显示末端时间点
abline(h=0.5) 
#添加中位线
abline(v=270,lty=3)
legend("bottomleft",c("Male","Female"),col=c("red","blue"),lwd=2)
text(900,0.9,"p=0.001")

survdiff(Lusurv~Lung$sex)
#两组之间是否有统计学差异