我正在估计具有时变协变量的加速失效时间生存模型，并希望计算我估计模型的偏差残差。
对于固定时间协变量，这很容易通过survival::survreg和实现stats::residuals。这是从这个网页中获取的一个例子（数据是包ovarian中包含的数据集survival）：
survregLogLogistic <- survreg(Surv(futime, fustat) ~ ecog.ps + rx,
data=ovarian, dist = "loglogistic")
devRes <- residuals(survregLogLogistic, type="deviance")
但是，survival::survreg无法处理随时间变化的协变量（有关如何格式化此类数据以在 R 中进行生存分析的示例，请参见此处）。因此，我需要使用flexsurv::flexsurvregor eha::aftreg，如下所示（没有随时间变化的协变量，但这足以重现问题）：
flexsurvregLogLogistic <- flexsurvreg(Surv(futime, fustat) ~ ecog.ps + rx,
data=ovarian, dist = "llogis")
估计工作正常，但residuals命令不起作用（survival::residuals使用 one 仔细查看函数可以看到函数或methods(residuals)函数都没有子函数）。这是残差函数的输出：flexsurv::flexsurvregeha::aftreg
> 残差（flexsurvregLogLogistic，type="deviance"）
无效的
现在为eha::aftreg：
aftregLogLogistic <- aftreg(Surv(futime, fustat) ~ ecog.ps + rx,
data=ovarian, dist = "loglogistic")
> 残差(aftregLogLogistic, type="deviance")
无效的
如何计算任一对象的偏差（或潜在的鞅）flexsurv::flexsurvreg残差eha::aftreg？