DYNESTY: a dynamic nested sampling package for estimating Bayesian posteriors and evidences
Abstract
We present DYNESTY, a public, open-source, PYTHON package to estimate Bayesian posteriors and evidences (marginal likelihoods) using the dynamic nested sampling methods developed by Higson et al. By adaptively allocating samples based on posterior structure, dynamic nested sampling has the benefits of Markov chain Monte Carlo (MCMC) algorithms that focus exclusively on posterior estimation while retaining nested sampling's ability to estimate evidences and sample from complex, multimodal distributions. We provide an overview of nested sampling, its extension to dynamic nested sampling, the algorithmic challenges involved, and the various approaches taken to solve them in this and previous work. We then examine DYNESTY's performance on a variety of toy problems along with several astronomical applications. We find in particular problems DYNESTY can provide substantial improvements in sampling efficiency compared to popular MCMC approaches in the astronomical literature. More detailed statistical results related to nested sampling are also included in the appendix.
- Publication:
-
Monthly Notices of the Royal Astronomical Society
- Pub Date:
- April 2020
- DOI:
- 10.1093/mnras/staa278
- arXiv:
- arXiv:1904.02180
- Bibcode:
- 2020MNRAS.493.3132S
- Keywords:
-
- methods: data analysis;
- methods: statistical;
- Astrophysics - Instrumentation and Methods for Astrophysics;
- Statistics - Computation
- E-Print:
- 28 pages, 12 figures, submitted to MNRAS