DNest4: Diffusive Nested Sampling in C++ and Python
arXiv:1606.03757
Abstract
In probabilistic (Bayesian) inferences, we typically want to compute properties of the posterior distribution, describing knowledge of unknown quantities in the context of a particular dataset and the assumed prior information. The marginal likelihood, also known as the "evidence", is a key quantity in Bayesian model selection. The Diffusive Nested Sampling algorithm, a variant of Nested Sampling, is a powerful tool for generating posterior samples and estimating marginal likelihoods. It is effective at solving complex problems including many where the posterior distribution is multimodal or has strong dependencies between variables. DNest4 is an open source (MIT licensed), multi-threaded implementation of this algorithm in C++11, along with associated utilities including: i) RJObject, a class template for finite mixture models, (ii) A Python package allowing basic use without C++ coding, and iii) Experimental support for models implemented in Julia. In this paper we demonstrate DNest4 usage through examples including simple Bayesian data analysis, finite mixture models, and Approximate Bayesian Computation.
Submitted. 33 pages, 9 figures. v2 removed a duplicated figure, v3 added a comparison to other packages, v4 fixed a few minor issues
References in corpus (7)
- PolyChord: next-generation nested sampling
- Updating Probabilities
- Modeling reverberation mapping data I: improved geometric and dynamical models and comparison with cross-correlation results
- Fast Bayesian Inference for Exoplanet Discovery in Radial Velocity Data
- Dissecting magnetar variability with Bayesian hierarchical models
- Properties of the Affine Invariant Ensemble Sampler in high dimensions
- Inference for Trans-dimensional Bayesian Models with Diffusive Nested Sampling
Cited by in corpus (6)
- Dynamic nested sampling: an improved algorithm for parameter estimation and evidence calculation
- Model selection and parameter inference in phylogenetics using Nested Sampling
- Constraining the limiting brightness temperature and Doppler factors for the largest sample of radio bright blazars
- On parametrised cold dense matter equation of state inference
- Nested sampling with plateaus
- Identifying changing jets through their radio variability