r - R Flexsurv- 估计生存函数

Question

我有一条生存 kapplan-Meier 曲线，我想使用 flexsurv 包推断不同的模型曲线（例如 Weibull、Gompertz 等）。我已经成功地进行了外推，但我没有找到创建外推图矩阵的解决方案。

library(survival)
library(flexsurv)
kmsurvival <- survfit (Surv(time, status) ~ 1, data=lung)
summary(kmsurvival)
plot(kmsurvival, xlab="Time", ylab="Survival probability")
Gompertz<-flexsurvreg(Surv(time, status)~1, data=lung, dist="gompertz")
plot(Gompertz)

我想创建 KM 生存曲线的输出和图中的外推。

例如，KM 曲线（20 个第一个时间点）：

v1 <- rep(NA,20)
v2<-1:20
for(i in 1:20){
v1[i] <- summary(kmsurvival, i)$surv
i=i+1
}
m1KM<-data.frame(v2,v1)

我想对 Gompertz 外推法做同样的事情，但我无法访问在每个时间点应用此曲线的结果。有什么帮助？？？谢谢！

Answer 1

est您可以在其摘要列中访问 Gompertz 的预测值：

m1G <- summary(Gompertz)[[1]]
plot(est~time, data=m1G)

如果需要计算与原始数据不同时间点的函数，可以使用

t <- 0:1000
summary(Gompertz, t=t)[[1]][,"est"]