如何使用 survreg 模型预测特定时间点的存活率？

问题描述

数据

library(survival)
kidney

型号

model = survreg(Surv(time,censored) ~ sex + age,data = kidney)

Call:
survreg(formula = Surv(time,data = kidney)

Coefficients:
(Intercept)   sexfemale         age 
 8.44411429 -0.89481679 -0.02170266 

Scale= 1.653512 

Loglik(model)= -122.1   Loglik(intercept only)= -122.7
    Chisq= 1.21 on 2 degrees of freedom,p= 0.547 
n= 76 

如何预测多个时间点（例如 30、90、182 天）两性的存活率（加上 95% CI）？

在不同的尺度（例如原始时间尺度、概率）上做这件事有什么技巧吗？

示例代码或示例将不胜感激。

解决方法

您可以使用 survminer 包。示例：

library(survival)
library(survminer)

f1 <- survfit(Surv(time,status) ~ sex,data = kidney)

res.sum <- surv_summary(f1,data = kidney)

# define desired time points
times <- c(30,90,182)

summary(f1,times)