R Flexsurv- оценка функции выживания

Question

У меня есть кривые Капплана-Мейера на выживание, и я хочу экстраполировать различные модельные кривые (например, Weibull, Gompertz,...), используя пакет flexsurv. Мне удалось выполнить экстраполяцию, но я не нашел решения для создания матрицы графика экстраполяции.

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)

Я хотел бы создать вывод кривой выживания КМ и экстраполяции, как на графиках.

Например, с кривой КМ (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)

Я хотел бы сделать то же самое с экстраполяцией Гомперца, но я не могу получить доступ к результатам, применяя эту кривую для каждой временной точки. Любая помощь??? Спасибо!

Answer 1

Вы можете получить доступ к прогнозируемым значениям Гомперца в est колонка его резюме:

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

Если вам нужно рассчитать функцию в разные моменты времени, чем исходные данные, вы можете использовать

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