Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection

Abstract

Standard approaches to Bayesian parameter inference in large scale structure assume a Gaussian functional form (chi-squared form) for the likelihood. This assumption, in detail, cannot be correct. Likelihood free inferences such as Approximate Bayesian Computation (ABC) relax these restrictions and make inference possible without making any assumptions on the likelihood. Instead ABC relies on a forward generative model of the data and a metric for measuring the distance between the model and data. In this work, we demonstrate that ABC is feasible for LSS parameter inference by using it to constrain parameters of the halo occupation distribution (HOD) model for populating dark matter halos with galaxies. Using specific implementation of ABC supplemented with Population Monte Carlo importance sampling, a generative forward model using HOD, and a distance metric based on galaxy number density, two-point correlation function, and galaxy group multiplicity function, we constrain the HOD parameters of mock observation generated from selected "true" HOD parameters. The parameter constraints we obtain from ABC are consistent with the "true" HOD parameters, demonstrating that ABC can be reliably used for parameter inference in LSS. Furthermore, we compare our ABC constraints to constraints we obtain using a pseudo-likelihood function of Gaussian form with MCMC and find consistent HOD parameter constraints. Ultimately our results suggest that ABC can and should be applied in parameter inference for LSS analyses.

16 pages, 10 figures

References in corpus (10)

• Why your model parameter confidences might be too optimistic -- unbiased estimation of the inverse covariance matrix
• Galaxy Evolution from Halo Occupation Distribution Modeling of DEEP2 and SDSS Galaxy Clustering
• Component separation methods for the Planck mission
• Introducing Decorated HODs: modeling assembly bias in the galaxy-halo connection
• The Cosmic Code Comparison Project
• Nucleosynthesis in a Primordial Supernova: Carbon and Oxygen Abundances in SMSS J031300.36-670839.31
• The Effect of Fiber Collisions on the Galaxy Power Spectrum Multipole
• nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics
• Bayesian inference of cosmic density fields from non-linear, scale-dependent, and stochastic biased tracers
• Satellite Kinematics I: A New Method to Constrain the Halo Mass-Luminosity Relation of Central Galaxies

Cited by in corpus (39)

• Warm dark matter chills out: constraints on the halo mass function and the free-streaming length of dark matter with 8 quadruple-image strong gravitational lenses
• 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
• Probing dark matter structure down to $10^7$ solar masses: flux ratio statistics in gravitational lenses with line of sight halos
• 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
• Hack Weeks as a model for Data Science Education and Collaboration
• Bayesian optimisation for likelihood-free cosmological inference
• Cosmology From CMB Lensing and Delensed EE Power Spectra Using 2019-2020 SPT-3G Polarization Data
• Simulation-Based Inference of Reionization Parameters From 3D Tomographic 21 cm Lightcone Images
• Probing the nature of dark matter by forward modeling flux ratios in strong gravitational lenses
• Inference of the optical depth to reionization from low multipole temperature and polarisation Planck data
• A New Method to Measure the Post-Reionization Ionizing Background from the Joint Distribution of Lyman-$\alpha$ and Lyman-$\beta$ Forest Transmission
• The multiplicity distribution of Kepler's exoplanets
• Implicit Likelihood Inference of Reionization Parameters from the 21 cm Power Spectrum
• Accelerated Bayesian SED Modeling using Amortized Neural Posterior Estimation
• How are galaxies assigned to halos? Searching for assembly bias in the SDSS galaxy clustering
• Machine Learning for Observational Cosmology
• The sum of the masses of the Milky Way and M31: a likelihood-free inference approach
• Likelihood Non-Gaussianity in Large-Scale Structure Analyses
• $\mathrm{S}\mathrm{IM}\mathrm{BIG}$: Mock Challenge for a Forward Modeling Approach to Galaxy Clustering
• Accelerating Approximate Bayesian Computation with Quantile Regression: Application to Cosmological Redshift Distributions
• How proper are Bayesian models in the astronomical literature?
• The PAU Survey: Measurement of Narrow-band galaxy properties with Approximate Bayesian Computation
• Inferring the photometric and size evolution of galaxies from image simulations
• BASS. XXXVI. Constraining the Local Supermassive Black Hole - Halo Connection with BASS DR2 AGN
• CoLFI: Cosmological Likelihood-free Inference with Neural Density Estimators
• Accurate X-ray Timing in the Presence of Systematic Biases With Simulation-Based Inference
• Machine Learning Accelerated Likelihood-Free Event Reconstruction in Dark Matter Direct Detection
• Likelihood-free Forward Modeling for Cluster Weak Lensing and Cosmology
• A Preferential Attachment Model for the Stellar Initial Mass Function
• Pushing the Limits of Detectability: Mixed Dark Matter from Strong Gravitational Lenses
• A forward-modelling method to infer the dark matter particle mass from strong gravitational lenses
• Minimum-entropy constraints on galactic potentials
• Trans-Neptunian Space and the Post-Pluto Paradigm
• Quantifying Weighted Morphological Content of Large-Scale Structures via Simulation-Based Inference
• IQ Collaboratory III: The Empirical Dust Attenuation Framework -- Taking Hydrodynamical Simulations with a Grain of Dust