r - 在 R 中使用 flexsurvreg 预测样本外

Question

我在 R 中有以下模型

library(flexsurv)

data(ovarian)

model = flexsurvreg(Surv(futime, fustat) ~ ecog.ps + rx, data = ovarian, dist='weibull')

model

predict(model,data = ovarian, type = 'response')

模型摘要如下所示 flexsurvreg model output

我正在尝试使用 R 中的预测函数预测生存时间并得到以下错误

error while trying to predict

如何使用此 flexsurvreg 模型预测预期生命周期？

我知道文档中提到了一个totlos.fs 函数，但是这个数据似乎没有一个trans 变量totlos.fs> 需要提供输出。

如果除了 totlos.fs 没有其他选择，我如何在该数据中创建一个 trans 变量并将其与现有协变量一起处理？

请指教。

Answer 1

supplementary examples doc 的第 3 节对于 flexsurv文档有一个示例，其中使用模型方程直接计算预测值。当您使用 Weibull 分布(n=2 参数)时，我相信这应该有效:

pred.model <- model.matrix(model) %*% model$res[-(1:n),"est"]

干杯