Cahiers du CEREMADE

Unité Mixte de Recherche du C.N.R.S. N°7534
Abstract : Nested sampling is a novel simulation method for approximating marginal likelihoods, proposed by Skilling (1006,2007). We establish that nested sampling leads to an error that vanishes at the standard Monte Carlo rate $N^{-1/2}$, where $N$ is a tuning parameter that is proportional to the computational effort, and that this error is asymptotically Gaussian. We show that the corresponding asymptotic variance typically grows linearly with the dimension of the parameter. We use these results to discuss the applicability and efficiency of nested sampling in realistic problems, including posterior distributions for mixtures. We propose an extension of nested sampling that makes it possible to avoid resorting to MCMC to obtain the simulated points. We study two alternative methods for computing marginal likelihood, which, in contrast with nested sampling, are based on draws from the posterior distribution and we conduct a comparison with nested sampling on several realistic examples
Contemplating Evidence: properties, extensions of, and alternatives to Nested Sampling
CHOPIN Nicolas, ROBERT Christian P.
Université de PARIS - DAUPHINE
Place du Maréchal de Lattre De Tassigny - 75775 PARIS CEDEX 16 - FRANCE
Téléphone : +33 (0)1 44-05-49-23 - fax : +33 (0)1 44-05-45-99