Nested sampling cross-checks using order statistics
- 11 August 2020
- journal article
- research article
- Published by Oxford University Press (OUP) in Monthly Notices of the Royal Astronomical Society
- Vol. 497 (4) , 5256-5263
Abstract
Nested sampling (NS) is an invaluable tool in data analysis in modern astrophysics, cosmology, gravitational wave astronomy, and particle physics. We identify a previously unused property of NS related to order statistics: the insertion indexes of new live points into the existing live points should be uniformly distributed. This observation enabled us to create a novel cross-check of single NS runs. The tests can detect when an NS run failed to sample new live points from the constrained prior and plateaus in the likelihood function, which break an assumption of NS and thus leads to unreliable results. We applied our cross-check to NS runs on toy functions with known analytic results in 2–50 dimensions, showing that our approach can detect problematic runs on a variety of likelihoods, settings, and dimensions. As an example of a realistic application, we cross-checked NS runs performed in the context of cosmological model selection. Since the cross-check is simple, we recommend that it become a mandatory test for every applicable NS run.Keywords
All Related Versions
Funding Information
- National Natural Science Foundation of China (11950410509)
- University of Adelaide
- Science and Technology Facilities Council (G102229)
This publication has 43 references indexed in Scilit:
- Conservative constraints on early cosmology with MONTE PYTHONJournal of Cosmology and Astroparticle Physics, 2013
- Nested sampling for parameter inference in systems biology: application to an exemplar circadian modelBMC Systems Biology, 2013
- Bayesian analysis of inflation. II. Model selection and constraints on reheatingPhysical Review D, 2012
- CONSTRAINTS ON COSMIC-RAY PROPAGATION MODELS FROM A GLOBAL BAYESIAN ANALYSISThe Astrophysical Journal, 2011
- Efficient Sampling of Atomic Configurational SpacesThe Journal of Physical Chemistry B, 2010
- MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physicsMonthly Notices of the Royal Astronomical Society, 2009
- The impact of priors and observables on parameter inferences in the constrained MSSMJournal of High Energy Physics, 2008
- Multimodal nested sampling: an efficient and robust alternative to Markov Chain Monte Carlo methods for astronomical data analysesMonthly Notices of the Royal Astronomical Society, 2008
- Information criteria for astrophysical model selectionMonthly Notices of the Royal Astronomical Society: Letters, 2007
- Bayes FactorsJournal of the American Statistical Association, 1995