11
15

Vertical-likelihood Monte Carlo

Abstract

In this review, we address the use of Monte Carlo methods for approximating definite integrals of the form Z=L(x)dP(x)Z = \int L(x) d P(x), where LL is a target function (often a likelihood) and PP a finite measure. We present vertical-likelihood Monte Carlo, which is an approach for designing the importance function g(x)g(x) used in importance sampling. Our approach exploits a duality between two random variables: the random draw XgX \sim g, and the corresponding random likelihood ordinate YL(X)Y\equiv L(X) of the draw. It is natural to specify g(x)g(x) and ask: what is the the implied distribution of YY? In this paper, we take up the opposite question: what should the distribution of YY be so that the implied importance function g(x)g(x) is good for approximating ZZ? Our answer turns out to unite seven seemingly disparate classes of algorithms under the vertical-likelihood perspective: importance sampling, slice sampling, simulated annealing/tempering, the harmonic-mean estimator, the vertical-density sampler, nested sampling, and energy-level sampling (a suite of related methods from statistical physics). In particular, we give an alterate presentation of nested sampling, paying special attention to the connection between this method and the vertical-likelihood perspective articulated here. As an alternative to nested sampling, we describe an MCMC method based on re-weighted slice sampling. This method's convergence properties are studied, and two examples demonstrate the promise of the overall approach.

View on arXiv
Comments on this paper