Clustered nested sampling: efficient Bayesian inference for cosmology
Ruth G. Shaw, M. Bridges, Michael P. Hobson · arXiv (Cornell University) · 2007
Bayesian model selection provides the cosmologist with an exacting tool to distinguish between competing models based purely on the data, via the Bayesian evidence. Previous methods to calculate this quantity either lacked general applicability or were computationally demanding. However, nested sampling (Skilling 2004), which was recently applied successfully to cosmology by Muhkerjee et al. 2006, overcomes both of these impediments. Their implementation restricts the parameter space sampled, and thus improves the efficiency, using a decreasing ellipsoidal bound in the n-dimensional parameter space centred on the maximum likelihood point. However, if the likelihood function contains any multi-modality, separated over a significant portion of the parameter space then the ellipse is prevented from constraining the sampling region by less than the distance between the likelihood peaks. In this paper we introduce a method of clustered nested sampling whereby ellipsoidal clusters can form on any peaks identified –thus improving the efficiency by a factor which is equal to the ratio of the volumes enclosed by the set of small clustered ellipsoids and the large single ellipse that would necessarily be required without clustering. In addition we have implemented a method for determining