Gaussbock: Fast parallel-iterative cosmological parameter estimation with Bayesian nonparametrics
arXiv:1905.09800 · doi:10.3847/1538-4357/ab93cb
Abstract
We present and apply Gaussbock, a new embarrassingly parallel iterative algorithm for cosmological parameter estimation designed for an era of cheap parallel computing resources. Gaussbock uses Bayesian nonparametrics and truncated importance sampling to accurately draw samples from posterior distributions with an orders-of-magnitude speed-up in wall time over alternative methods. Contemporary problems in this area often suffer from both increased computational costs due to high-dimensional parameter spaces and consequent excessive time requirements, as well as the need for fine tuning of proposal distributions or sampling parameters. Gaussbock is designed specifically with these issues in mind. We explore and validate the performance and convergence of the algorithm on a fast approximation to the Dark Energy Survey Year 1 (DES Y1) posterior, finding reasonable scaling behavior with the number of parameters. We then test on the full DES Y1 posterior using large-scale supercomputing facilities, and recover reasonable agreement with previous chains, although the algorithm can underestimate the tails of poorly-constrained parameters. Additionally, we discuss and demonstrate how Gaussbock recovers complex posterior shapes very well at lower dimensions, but faces challenges to perform well on such distributions in higher dimensions. In addition, we provide the community with a user-friendly software tool for accelerated cosmological parameter estimation based on the methodology described in this paper.
arXiv comment19 pages, 10 figures, accepted for publication in ApJ
References in corpus (41)
- emcee: The MCMC Hammer
- Variational Inference: A Review for Statisticians
- Cosmological parameters from CMB and other data: a Monte-Carlo approach
- MCMC using Hamiltonian dynamics
- MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
- Stable clustering, the halo model and nonlinear cosmological power spectra
- Improved cosmological constraints from a joint analysis of the SDSS-II and SNLS supernova samples
- Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
- Bayes in the sky: Bayesian inference and model selection in cosmology
- Dark Energy Survey Year 1 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing
- Revising the Halofit Model for the Nonlinear Matter Power Spectrum
- The Dark Energy Survey: more than dark energy - an overview
- Cosmology and Fundamental Physics with the Euclid Satellite
- CMB power spectrum parameter degeneracies in the era of precision cosmology
- PolyChord: next-generation nested sampling
- A Nested Sampling Algorithm for Cosmological Model Selection
- Constraints on cosmic-ray propagation models from a global Bayesian analysis
- CosmoSIS: modular cosmological parameter estimation
- Inference of the cosmological parameters from gravitational waves: application to second generation interferometers
- A Short History of Markov Chain Monte Carlo: Subjective Recollections from Incomplete Data
- CosmoLike - Cosmological Likelihood Analyses for Photometric Galaxy Surveys
- Cosmology from Cosmic Shear with DES Science Verification Data
- II: Bayesian Methods for Cosmological Parameter Estimation from Cosmic Microwave Background Measurements
- Massive optimal data compression and density estimation for scalable, likelihood-free inference in cosmology
- Cosmology constraints from shear peak statistics in Dark Energy Survey Science Verification data
- Properties of Nested Sampling
- The Age of the Universe and the Cosmological Constant Determined from Cosmic Microwave Background Anisotropy Measurements
- Approximate Bayesian Computation for Forward Modeling in Cosmology
- CosmoHammer: Cosmological parameter estimation with the MCMC Hammer
- Estimation of cosmological parameters using adaptive importance sampling
- cosmoabc: Likelihood-free inference via Population Monte Carlo Approximate Bayesian Computation
- Cosmological inference using only gravitational wave observations of binary neutron stars
- Comparison of sampling techniques for Bayesian parameter estimation
- Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
- Sampling Errors in Nested Sampling Parameter Estimation
- Sampling using a `bank' of clues
- Constraints on the sum of neutrino masses using cosmological data including the latest extended Baryon Oscillation Spectroscopic Survey DR14 quasar sample
- (Lack of) Cosmological evidence for dark radiation after Planck
- On statistical uncertainty in nested sampling
- Bayesian model comparison in cosmology with Population Monte Carlo
- Dark energy constraints and correlations with systematics from CFHTLS weak lensing, SNLS supernovae Ia and WMAP5