Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
read the original abstract
In performing a Bayesian analysis of astronomical data, two difficult problems often emerge. First, in estimating the parameters of some model for the data, the resulting posterior distribution may be multimodal or exhibit pronounced (curving) degeneracies, which can cause problems for traditional MCMC sampling methods. Second, in selecting between a set of competing models, calculation of the Bayesian evidence for each model is computationally expensive. The nested sampling method introduced by Skilling (2004), has greatly reduced the computational expense of calculating evidences and also produces posterior inferences as a by-product. This method has been applied successfully in cosmological applications by Mukherjee et al. (2006), but their implementation was efficient only for unimodal distributions without pronounced degeneracies. Shaw et al. (2007), recently introduced a clustered nested sampling method which is significantly more efficient in sampling from multimodal posteriors and also determines the expectation and variance of the final evidence from a single run of the algorithm, hence providing a further increase in efficiency. In this paper, we build on the work of Shaw et al. and present three new methods for sampling and evidence evaluation from distributions that may contain multiple modes and significant degeneracies; we also present an even more efficient technique for estimating the uncertainty on the evaluated evidence. These methods lead to a further substantial improvement in sampling efficiency and robustness, and are applied to toy problems to demonstrate the accuracy and economy of the evidence calculation and parameter estimation. Finally, we discuss the use of these methods in performing Bayesian object detection in astronomical datasets.
This paper has not been read by Pith yet.
Forward citations
Cited by 10 Pith papers
-
Symbolic Classification-Enabled LHC Limits Online BSM Global Fits
Symbolic regression produces an approximate classifier for LHC exclusion limits that enables their direct inclusion during pMSSM global fits.
-
Numerical polology: towards next-generation model-building for cosmology
Numerical polology framework samples coupling space to discover ghost-free tensor field theories up to rank three for cosmology, then applies resulting priors to black hole superradiance, dynamical dark energy, and GW data.
-
Non-holomorphic $S^{\prime}_{4}$ modular symmetry for leptons and leptogenesis
36 viable non-holomorphic S'4 modular models for leptons are identified via numerical scans, with two yielding successful unflavored thermal leptogenesis from the real part of τ while fitting neutrino data.
-
Analytical and Machine Learning Methods for Model Discernment at CE$\nu$NS Experiments
Shape correlations in CEνNS allow likelihood and CNN analyses to discriminate sterile neutrinos from NSI and approximately localize sterile parameters in favorable regions.
-
The Majoron Cosmological Window: Dark Matter and Thermal Leptogenesis
The minimal majoron framework permits simultaneous majoron dark matter and thermal leptogenesis in a constrained cosmological window set by freeze-in production, warm dark matter bounds, and indirect detection limits.
-
Constraints on Primordial Black Holes from Galactic Diffuse Synchrotron Emissions
Galactic synchrotron emissions above 20 MHz can set tighter upper limits on the abundance of primordial black holes with masses above 10^16 grams than previous cosmic-ray electron data.
-
GetDist: a Python package for analysing Monte Carlo samples
GetDist implements boundary-corrected KDE with automatic smoothing for analyzing weighted and correlated Monte Carlo samples, plus plotting and diagnostic tools.
-
CMB constraints on ultra-light primordial black holes with extended mass distributions
For ultra-light PBHs with extended mass distributions, new CMB-derived 95% exclusion bounds on f_PBH are obtained by jointly varying ΛCDM parameters, yielding f_PBH < 1.6 × 10^{-5} for a uniform distribution between 1...
-
The impact of FRB dispersion measure probability distribution functions on cosmographic estimates
Cosmographic constraints from 106 FRBs on H0, q0, j0 depend sensitively on whether a Gaussian or quasi-Gaussian PDF is assumed for DM_IGM and on whether the baryon mass fraction is fixed or free.
-
BSMArt 2: simpler and faster parameter space scans
BSMArt version 2 adds new scanning algorithms including Affine MC, MLScanner, and CMA-ES variants to simplify and accelerate parameter space exploration in new physics models, demonstrated on soft lepton excess search...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.