R语言-生存分析

有关生存分析的相关知识，请自行前往“算法与模型”类博客进行学习。
(http://blog.csdn.net/xiaohukun/article/details/77679134)

一、下载并加载包

进行生存分析使用的是survival 包

install.packages("survival")   #安装survival包
library(survival)      #加载survival包

二、准备数据

使用survival包自带的“pbc”数据集为例（418*20）

这里记录的time只是原始时间，还需要根据status判断生存时间的类型，转化成相应的格式。这个工作由Surv(time,event)完成,并返回生存对象

> Sur_Obj<-Surv(pbc$time,pbc$status)
> Sur_Obj
[1]  400  4500? 1012  1925  1504+ 2503  1832? 2466  2400 
 [10]   51  3762   304  3577? 1217  3584  3672?  769   131 ……

三、绘制生存曲线

函数surfit(formula)，根据填