Estimation of cosmological parameters using adaptive importance sampling

Abstract

We present a Bayesian sampling algorithm called adaptive importance sampling or Population Monte Carlo (PMC), whose computational workload is easily parallelizable and thus has the potential to considerably reduce the wall-clock time required for sampling, along with providing other benefits. To assess the performance of the approach for cosmological problems, we use simulated and actual data consisting of CMB anisotropies, supernovae of type Ia, and weak cosmological lensing, and provide a comparison of results to those obtained using state-of-the-art Markov Chain Monte Carlo (MCMC). For both types of data sets, we find comparable parameter estimates for PMC and MCMC, with the advantage of a significantly lower computational time for PMC. In the case of WMAP5 data, for example, the wall-clock time reduces from several days for MCMC to a few hours using PMC on a cluster of processors. Other benefits of the PMC approach, along with potential difficulties in using the approach, are analysed and discussed.

17 pages, 11 figures

References in corpus (27)

• Efficient Computation of CMB anisotropies in closed FRW models
• Cosmological parameters from CMB and other data: a Monte-Carlo approach
• The Supernova Legacy Survey: Measurement of Omega_M, Omega_Lambda and w from the First Year Data Set
• Stable clustering, the halo model and nonlinear cosmological power spectra
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Data Processing, Sky Maps, and Basic Results
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Likelihoods and Parameters from the WMAP data
• Accurate photometric redshifts for the CFHT Legacy Survey calibrated using the VIMOS VLT Deep Survey
• Three Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Polarization Analysis
• Three-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Temperature Analysis
• Power Spectrum Tomography with Weak Lensing
• Very weak lensing in the CFHTLS Wide: Cosmology from cosmic shear in the linear regime
• Cosmic Confusion: Degeneracies among Cosmological Parameters Derived from Measurements of Microwave Background Anisotropies
• A new measure for cosmic shear
• Cosmological parameters from lensing power spectrum and bispectrum tomography
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Angular Power Spectra
• Fast and reliable MCMC for cosmological parameter estimation
• First cosmic shear results from the Canada-France-Hawaii Telescope Wide Synoptic Legacy Survey
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Galactic Foreground Emission
• Efficient Bayesian inference for multimodal problems in cosmology
• Dark Energy Constraints from the CTIO Lensing Survey
• Cosmological parameters from combined second- and third-order aperture mass statistics of cosmic shear
• Fast cosmological parameter estimation using neural networks
• Estimation of Polarized Power Spectra by Gibbs sampling
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• The nonlinear matter power spectrum
• CMB likelihood approximation by a Gaussianized Blackwell-Rao estimator

Cited by in corpus (46)

• CFHTLenS tomographic weak lensing cosmological parameter constraints: Mitigating the impact of intrinsic galaxy alignments
• Cosmology with cosmic shear observations: a review
• Evidence for the accelerated expansion of the Universe from weak lensing tomography with COSMOS
• CFHTLenS: Combined probe cosmological model comparison using 2D weak gravitational lensing
• Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy
• General Relativity and Cosmology: Unsolved Questions and Future Directions
• Cosmology with the shear-peak statistics
• Galaxy clustering in the CFHTLS-Wide: the changing relationship between galaxies and haloes since z ~ 1.2
• Mass Calibration and Cosmological Analysis of the SPT-SZ Galaxy Cluster Sample Using Velocity Dispersion ${\sigma }_v$ and X-ray $Y_{\text{X}}$ Measurements
• The galaxy-halo connection from a joint lensing, clustering and abundance analysis in the CFHTLenS/VIPERS field
• Precision calculations of the cosmic shear power spectrum projection
• CFHTLenS tomographic weak lensing: Quantifying accurate redshift distributions
• Approximate Bayesian Computation for Astronomical Model Analysis: A Case Study in Galaxy Demographics and Morphological Transformation at High Redshift
• Approximate Bayesian Computation for Forward Modeling in Cosmology
• The phenomenological approach to modeling the dark energy
• Probing the galaxy-halo connection in UltraVISTA to $z\sim 2$
• A new model to predict weak-lensing peak counts II. Parameter constraint strategies
• Cosmological constraints with weak lensing peak counts and second-order statistics in a large-field survey
• A Dust spectral energy distribution model with hierarchical Bayesian inference. I. Formalism & benchmarking
• The evolution of clustering length, large-scale bias and host halo mass at 2<z<5 in the VIMOS Ultra Deep Survey (VUDS)
• TEASING: a fast and accurate approximation for the low multipole likelihood of the Cosmic Microwave Background temperature
• Cosmological parameter estimation via iterative emulation of likelihoods
• Forecasts of non-Gaussian parameter spaces using Box-Cox transformations
• The VIMOS Ultra Deep Survey. Luminosity and stellar mass dependence of galaxy clustering at z~3
• The Galaxy--Halo Connection in High-Redshift Universe: Details and Evolution of Stellar-to-Halo Mass Ratios of Lyman Break Galaxies on CFHTLS Deep Fields
• A fast, always positive definite and normalizable approximation of non-Gaussian likelihoods
• Hierarchical probabilistic inference of cosmic shear
• Cosmological parameter estimation with Genetic Algorithms
• On the impact of intergalactic dust on cosmology with type Ia supernovae
• Bayesian model comparison in cosmology with Population Monte Carlo
• Cosmology with the largest galaxy cluster surveys: Going beyond Fisher matrix forecasts
• Properties and uncertainties of scalar field models of dark energy with barotropic equation of state
• Halo-model analysis of the clustering of photometric luminous red galaxies at $0.10\le z\le 1.05$ from the Subaru Hyper Suprime-Cam Survey
• Approximating Probability Densities by Iterated Laplace Approximations
• Initializing adaptive importance sampling with Markov chains
• A likely magnetic activity cycle for the exoplanet host M dwarf GJ 3512
• Euclid preparation. LIII. LensMC, weak lensing cosmic shear measurement with forward modelling and Markov Chain Monte Carlo sampling
• Posterior inference unchained with EL_2O
• The very wide-field $gzK$ galaxy survey -- II. The relationship between star-forming galaxies at $z\sim 2$ and their host haloes based upon HOD modelling
• Bayesian model selection for dark energy using weak lensing forecasts
• A Bayesian inference and model selection algorithm with an optimisation scheme to infer the model noise power
• Gaussbock: Fast parallel-iterative cosmological parameter estimation with Bayesian nonparametrics
• Improved Gibbs samplers for Cosmic Microwave Background power spectrum estimation
• Using parallel computation to improve Independent Metropolis--Hastings based estimation
• Optimisation of the Population Monte Carlo algorithm: Application to constraining isocurvature models with cosmic microwave background data
• On computational tools for Bayesian data analysis