Accelerating Approximate Bayesian Computation with Quantile Regression: Application to Cosmological Redshift Distributions

Abstract

Approximate Bayesian Computation (ABC) is a method to obtain a posterior distribution without a likelihood function, using simulations and a set of distance metrics. For that reason, it has recently been gaining popularity as an analysis tool in cosmology and astrophysics. Its drawback, however, is a slow convergence rate. We propose a novel method, which we call qABC, to accelerate ABC with Quantile Regression. In this method, we create a model of quantiles of distance measure as a function of input parameters. This model is trained on a small number of simulations and estimates which regions of the prior space are likely to be accepted into the posterior. Other regions are then immediately rejected. This procedure is then repeated as more simulations are available. We apply it to the practical problem of estimation of redshift distribution of cosmological samples, using forward modelling developed in previous work. The qABC method converges to nearly same posterior as the basic ABC. It uses, however, only 20\% of the number of simulations compared to basic ABC, achieving a fivefold gain in execution time for our problem. For other problems the acceleration rate may vary; it depends on how close the prior is to the final posterior. We discuss possible improvements and extensions to this method.

10 pages, 5 figures, prepared for submission to JCAP

References in corpus (18)

• The First Release COSMOS Optical and Near-IR Data and Catalog
• The Cosmic Evolution Survey (COSMOS): Subaru Observations of the HST COSMOS Field
• Estimating conditional quantiles with the help of the pinball loss
• The Third Gravitational Lensing Accuracy Testing (GREAT3) Challenge Handbook
• Lensing substructure quantification in RXJ1131-1231: A 2 keV lower bound on dark matter thermal relic mass
• Likelihood-Free Cosmological Inference with Type Ia Supernovae: Approximate Bayesian Computation for a Complete Treatment of Uncertainty
• Approximate Bayesian Computation for Forward Modeling in Cosmology
• cosmoabc: Likelihood-free inference via Population Monte Carlo Approximate Bayesian Computation
• astroABC: An Approximate Bayesian Computation Sequential Monte Carlo sampler for cosmological parameter estimation
• Accelerating ABC methods using Gaussian processes
• Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection
• Bayesian Modeling with Gaussian Processes using the GPstuff Toolbox
• The redshift distribution of cosmological samples: a forward modeling approach
• liquidSVM: A Fast and Versatile SVM package
• The $z<1.2$ optical luminosity function from a sample of $\sim 410000$ galaxies in bootes
• A new approach for obtaining cosmological constraints from Type Ia Supernovae using Approximate Bayesian Computation
• Simulation-based model selection for dynamical systems in systems and population biology
• Semi-automatic selection of summary statistics for ABC model choice

Cited by in corpus (21)

• Fast likelihood-free cosmology with neural density estimators and active learning
• Massive optimal data compression and density estimation for scalable, likelihood-free inference in cosmology
• Likelihood-free inference with neural compression of DES SV weak lensing map statistics
• Nuisance hardened data compression for fast likelihood-free inference
• Automatic physical inference with information maximising neural networks
• Generalized massive optimal data compression
• Bayesian optimisation for likelihood-free cosmological inference
• Simulation-Based Inference of Reionization Parameters From 3D Tomographic 21 cm Lightcone Images
• Inference of the optical depth to reionization from low multipole temperature and polarisation Planck data
• Accelerated Bayesian SED Modeling using Amortized Neural Posterior Estimation
• The sum of the masses of the Milky Way and M31: a likelihood-free inference approach
• Likelihood Non-Gaussianity in Large-Scale Structure Analyses
• Measurement of the B-band Galaxy Luminosity Function with Approximate Bayesian Computation
• Exploring Nonlinear Electrodynamics Theories: Shadows of Regular Black Holes and Horizonless Ultra-Compact Objects
• The PAU Survey: Measurement of Narrow-band galaxy properties with Approximate Bayesian Computation
• 21st Century Statistical and Computational Challenges in Astrophysics
• DIGS: Deep Inference of Galaxy Spectra with Neural Posterior Estimation
• Likelihood-free Forward Modeling for Cluster Weak Lensing and Cosmology
• Catalog-free modeling of galaxy types in deep images: Massive dimensional reduction with neural networks
• Manifold regularization based on Nystr{ö}m type subsampling
• Observational distinguishability of the Kerr and Kerr-Hayward metrics to EHT