我正在尝试使用flexsurvreg
包中的某些生存数据来估计 Weibull 分布的平均故障时间flexsurv
。我需要能够估计用于模拟模型的标准误差。
以数据为例flexsurvreg
;lung
require(flexsurv)
lungS <- Surv(lung$time,lung$status)
lungfit <- flexsurvreg(lungS~1,dist="weibull")
lungfit
Call:
flexsurvreg(formula = lungS ~ 1, dist = "weibull")
Maximum likelihood estimates:
est L95% U95%
shape 1.32 1.14 1.52
scale 418.00 372.00 469.00
N = 228, Events: 165, Censored: 63
Total time at risk: 69593
Log-likelihood = -1153.851, df = 2
AIC = 2311.702
现在,计算平均值只是将估计的参数值插入标准公式的情况,但是有没有一种简单的方法可以得出这个估计的标准误差?可以survreg
这样做吗?