Я хотел бы использовать параметрические модели выживания, предполагающие распределение Гомперца с использованием пакета flexsurv
в R, для количественного определения ассоциаций в биомедицинских исследованиях.Однако, насколько мне известно, нет встроенной функции, такой как residuals.survreg
, содержащейся в пакете survival
, для извлечения остатков, позволяющих проводить диагностику модели (т. Е. Проверять пропорциональность предположений об опасностях, влиятельных личностей или отдаленных).случаи и т. д.)
Существуют ли способы определения пропорциональности опасностей и оценки наличия влиятельных случаев, выбросов, правильной ковариатной функциональной формы и т. д. с использованием следующего примера:
require(flexsurv)
# No scientific notation
options(scipen = 999)
# Use the breast cancer data contained within the flexsurv package
d <- bc
# Fit a parametric survival model with Gompertz distrubution
fit <- flexsurvreg(Surv(rectime,censrec==1) ~ group, dist = "gompertz",
data = d)
fit