A Comparison of two Sequential Metropolis-Hastings Algorithms with Standard Simulation Techniques in Bayesian Inference in Reliability Models Involving the Generalized Gamma Distribution

Abstract

In this paper we consider the generalized gamma distribution as introduced in Gasemyr and Natvig (1998). This distribution enters naturally in Bayesian inference in exponential survival models with left censoring. In the paper mentioned above it is shown that the weighted sum of products of generalized gamma distributions is a conjugate prior for the parameters of component lifetimes, having autopsy data in a Marshall-Olkin shock model. A corresponding result is shown in Gasemyr and Natvig (1999) for independent, exponentially distributed component lifetimes in a model with partial monitoring of components with applications to preventive system maintenance. A discussion in the present paper strongly indicates that expressing the posterior distribution in terms of the generalized gamma distribution is computationally efficient compared to using the ordinary gamma distribution in such models. Furthermore, we present two types of sequential Metropolis-Hastings algorithms that may be used in Bayesian inference in situations where exact methods are intractable. Finally these types of algorithms are compared with standard simulation techniques and analytical results in arriving at the posterior distribution of the parameters of component lifetimes in special cases of the mentioned models. It seems that one of these types of algorithms may be very favorable when prior assessments are updated by several data sets and when there are significant discrepancies between the prior assessments and the data.