Modeling complex life systems: Bayesian inference for Weibull failure times using adaptive MCMC
摘要
This study develops a Bayesian framework for analyzing failure times under the Weibull distribution, addressing challenges in prior specification caused by the absence of conjugate priors and multi-dimensional sufficient statistics. We propose an adaptive semi-parametric MCMC algorithm, combining a hierarchical Bayesian model with a noninformative LogNormal hyper-prior and the No-U-Turn Sampler (NUTS) in STAN. The framework is evaluated across 24 prior combinations and 72 simulated datasets. Results show that the proposed approach achieves up to 40% lower weighted relative efficiency (WRE) than classical methods such as maximum likelihood estimation (MLE), method of moments, and ordinary least squares, under both increasing and decreasing hazard rate scenarios. For large samples, priors such as Exponential–HalfCauchy yield WRE values as low as 0.83, demonstrating improved regularization. In a real-data application using prostate cancer survival data (