R语言-生存分析与结果的图像处理

R语言-生存分析与结果的图像处理

数据准备：

library("survival")
library("survminer")
data("lung")

调用“lung”数据集，使用head()命令调查前6行，得到以下结果：

 inst time status age sex ph.ecog ph.karno pat.karno meal.cal wt.loss
1    3  306      2  74   1       1       90       100     1175      NA
2    3  455      2  68   1       0       90        90     1225      15
3    3 1010      1  56   1       0       90        90       NA      15
4    5  210      2  57   1       1       90        60     1150      11
5    1  883      2  60   1       0      100        90       NA       0
6   12 1022      1  74   1       1       50        80      513       0

time是指生存时间，status是二值型的生存状态，其中“1”代表存活，“2”代表已经死亡。

KM非参数法：

调用survival包中的survfit来进行生存分析，得到每个时间段中的生存率。

fit <- survfit(Surv(time, status) ~ 1, data = lung)

直接查看fit的信息可以看到

> fit
Call: survfit(formula = Surv(time, status) ~ 1