我想使用参数生存模型,假设使用Gompertz分布并使用R中的flexsurv
包来量化生物医学研究中的关联。但是,据我所知,没有内置功能(例如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