An Application of Adaptive Independent Chain Metropolis–Hastings Algorithms in Bayesian Hazard Rate Estimation

In this paper, we consider adaptive independent chain (AIC) Metropolis–Hastings algorithms as introduced in a special context in Gåsemyr et al. (2001) and developed theoretically in Gåsemyr (2003). The algorithms aim at producing samples from a specific target distribution Π, and are adaptive, non-Markovian versions of the Metropolis–Hastings independent chain. A certain parametric class of possible proposal distributions is fixed, and the parameters of the proposal distribution are updated periodically on the basis of the recent history of the chain, thereby obtaining proposals that get ever closer to Π. In the former paper a version of these algorithms was shown to be very efficient compared to standard simulation techniques when applied to Bayesian inference in reliability models with at most three dependent parameters. The aim of the present paper is to investigate the performance of the AIC algorithm when the number of dependent parameters and the complexity of the model increases. As a test case we consider a model treated in Arjas and Gasbarra (1996). The target distribution Π is the posterior distribution for the vector X=(X1,...,Xn) of dependent parameters, representing a piecewise constant approximation to the hazard rate X(t), where t0≤ t≤ tn. Especially, for the case n=12 it turned out that some versions of the AIC were very efficient compared to standard simulation techniques and also to the algorithm applied in Arjas and Gasbarra (1996). This includes a version of the componentwise adaptive independent chain the basic idea of which was given in Gåsemyr (2003).