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"]