3

我正在尝试使用flexsurvreg包中的某些生存数据来估计 Weibull 分布的平均故障时间flexsurv。我需要能够估计用于模拟模型的标准误差。

以数据为例flexsurvreglung

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这样做吗?

4

1 回答 1

4

flexsurv0.2 版中,如果x是拟合模型对象,则x$cov是参数估计的协方差矩阵,在对数尺度上具有正参数。然后,您可以使用最大似然估计量的渐近正态属性。模拟大量多元法线向量,以估计值作为均值,以及这个协方差矩阵(例如使用rmvnorm包中的mvtnorm)。这为您提供了在采样不确定性下的参数估计值的复制。计算每个重复的相应平均存活率,然后取结果样本的 SD 或分位数以获得标准误差或置信区间。

于 2013-08-15T15:35:19.053 回答