## Stable Limit Laws for Harmonic Mean Estimators of Marginal Likelihoods (2010)

### BibTeX

@MISC{Wolpert10stablelimit,

author = {Robert L. Wolpert and Scott C. Schmidler},

title = {Stable Limit Laws for Harmonic Mean Estimators of Marginal Likelihoods},

year = {2010}

}

### OpenURL

### Abstract

The task of calculating marginal likelihoods arises in a wide array of statistical inference problems, including the evaluation of Bayes factors for model selection and hypothesis testing. Although Markov chain Monte Carlo methods have simplified many posterior calculations needed for practical Bayesian analysis, the evaluation of marginal likelihoods remains difficult. We consider the behavior of the wellknown harmonic mean estimator (Newton and Raftery, 1994) of the marginal likelihood, which converges almost-surely but may have infinite variance and so may not obey a central limit theorem. We give examples illustrating the convergence in distribution of the harmonic mean estimator to a one-sided stable law with characteristic exponent 1 < α < 2. While the harmonic mean estimator does converge almost surely, we show that it does so at rate n −ǫ where ǫ = 1 − α −1 is often as small as 0.10 or 0.01. In such a case, the reduction of Monte Carlo sampling error by a factor of two requires increasing the Monte Carlo sample size by a factor of 2 1/ǫ, or in excess of 2.5 ·10 30 when ǫ = 0.01. We explore the possibility of estimating the parameters of the limiting stable distribution to provide accelerated convergence.