Likelihood Non-Gaussianity in Large-Scale Structure Analyses

Abstract

Standard present day large-scale structure (LSS) analyses make a major assumption in their Bayesian parameter inference --- that the likelihood has a Gaussian form. For summary statistics currently used in LSS, this assumption, even if the underlying density field is Gaussian, cannot be correct in detail. We investigate the impact of this assumption on two recent LSS analyses: the Beutler et al. (2017) power spectrum multipole ($P_{\ell }$) analysis and the Sinha et al. (2017) group multiplicity function ($\zeta$) analysis. Using non-parametric divergence estimators on mock catalogs originally constructed for covariance matrix estimation, we identify significant non-Gaussianity in both the $P_{\ell }$ and $\zeta$ likelihoods. We then use Gaussian mixture density estimation and Independent Component Analysis on the same mocks to construct likelihood estimates that approximate the true likelihood better than the Gaussian $pseudo$-likelihood. Using these likelihood estimates, we accurately estimate the true posterior probability distribution of the Beutler et al. (2017) and Sinha et al. (2017) parameters. Likelihood non-Gaussianity shifts the $f{\sigma }_8$ constraint by $-0.44\sigma$, but otherwise, does not significantly impact the overall parameter constraints of Beutler et al. (2017). For the $\zeta$ analysis, using the pseudo-likelihood significantly underestimates the uncertainties and biases the constraints of Sinha et al. (2017) halo occupation parameters. For $\log M_1$ and $\alpha$, the posteriors are shifted by $+0.43\sigma$ and $-0.51\sigma$ and broadened by $42\%$ and $66\%$, respectively. The divergence and likelihood estimation methods we present provide a straightforward framework for quantifying the impact of likelihood non-Gaussianity and deriving more accurate parameter constraints.

arXiv comment33 pages, 7 figures

References in corpus (76)

• The cosmological simulation code GADGET-2
• Estimating Mutual Information
• The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological analysis of the DR12 galaxy sample
• Information criteria for astrophysical model selection
• Cosmology and Fundamental Physics with the Euclid Satellite
• Karhunen-Loeve eigenvalue problems in cosmology: how should we tackle large data sets?
• MultiDark simulations: the story of dark matter halo concentrations and density profiles
• The imprints of primordial non-gaussianities on large-scale structure: scale dependent bias and abundance of virialized objects
• Planck 2015 results. XI. CMB power spectra, likelihoods, and robustness of parameters
• Transients from Initial Conditions in Cosmological Simulations
• Planck 2013 results. XX. Cosmology from Sunyaev-Zeldovich cluster counts
• Planck 2013 results. XV. CMB power spectra and likelihood
• Constraints on local primordial non-Gaussianity from large scale structure
• Planck 2015 results. XXIV. Cosmology from Sunyaev-Zeldovich cluster counts
• Statistical Analysis of Galaxy Surveys - I. Robust error estimation for 2-point clustering statistics
• The DESI Experiment, a whitepaper for Snowmass 2013
• Percolation Galaxy Groups and Clusters in the SDSS Redshift Survey: Identification, Catalogs, and the Multiplicity Function
• The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Anisotropic galaxy clustering in Fourier-space
• The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: RSD measurement from the power spectrum and bispectrum of the DR12 BOSS galaxies
• The WiggleZ Dark Energy Survey: Improved Distance Measurements to z = 1 with Reconstruction of the Baryonic Acoustic Feature
• Transients from Initial Conditions: A Perturbative Analysis
• The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: mock galaxy catalogues for the BOSS Final Data Release
• The Bispectrum as a Signature of Gravitational Instability in Redshift-Space
• The WiggleZ Dark Energy Survey: Final data release and cosmological results
• The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Observational systematics and baryon acoustic oscillations in the correlation function
• nbodykit: an open-source, massively parallel toolkit for large-scale structure
• Parameter inference with estimated covariance matrices
• Myths and Truths Concerning Estimation of Power Spectra
• Fast Estimators for Redshift-Space Clustering
• The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Modeling the clustering and halo occupation distribution of BOSS-CMASS galaxies in the Final Data Release
• Data analysis recipes: Fitting a model to data
• The Bispectrum: From Theory to Observations
• Modelling Baryon Acoustic Oscillations with Perturbation Theory and Stochastic Halo Biasing
• Evidence for Low Black Hole Spin and Physically Motivated Accretion Models from Millimeter VLBI Observations of Sagittarius A*
• Constraining primordial non-Gaussianity with future galaxy surveys
• Massive optimal data compression and density estimation for scalable, likelihood-free inference in cosmology
• Extreme deconvolution: Inferring complete distribution functions from noisy, heterogeneous and incomplete observations
• The Clustering of Galaxies in SDSS-III DR9 Baryon Oscillation Spectroscopic Survey: Constraints on Primordial Non-Gaussianity
• Measuring line-of-sight dependent Fourier-space clustering using FFTs
• Galaxy And Mass Assembly: Deconstructing Bimodality - I. Red ones and blue ones
• The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: Cosmological implications of the Fourier space wedges of the final sample
• Detection of Baryon Acoustic Oscillation Features in the Large-Scale 3-Point Correlation Function of SDSS BOSS DR12 CMASS Galaxies
• Improved Primordial Non-Gaussianity Constraints from Measurements of Galaxy Clustering and the Integrated Sachs-Wolfe Effect
• A New Method to Correct for Fiber Collisions in Galaxy Two-Point Statistics
• Automatic physical inference with information maximising neural networks
• Dependence of cosmic shear covariances on cosmology - Impact on parameter estimation
• Breaking the Degeneracies Between Cosmology and Galaxy Bias
• Cosmological Structure Formation with Augmented Lagrangian Perturbation Theory
• Generalized massive optimal data compression
• Application of the Gaussian mixture model in pulsar astronomy -- pulsar classification and candidates ranking for {\it Fermi} 2FGL catalog
• An optimal FFT-based anisotropic power spectrum estimator
• Power Spectra Estimation for Weak Lensing
• P-MaNGA: Full spectral fitting and stellar population maps from prototype observations
• The Effect of Fiber Collisions on the Galaxy Power Spectrum Multipole
• The 3-point function in large scale structure: redshift distortions and galaxy bias
• A Machine Learning Approach for Dynamical Mass Measurements of Galaxy Clusters
• Correlations in the Spatial Power Spectrum Inferred from Angular Clustering: Methods and Application to APM
• Constraining the halo bispectrum in real and redshift space from perturbation theory and non-linear stochastic bias
• Hybrid Estimation of CMB Polarization Power Spectra
• Massive data compression for parameter-dependent covariance matrices
• Perturbative approach to covariance matrix of the matter power spectrum
• Nonparametric Divergence Estimation with Applications to Machine Learning on Distributions
• On the insufficiency of arbitrarily precise covariance matrices: non-Gaussian weak lensing likelihoods
• Dynamical Mass Measurements of Contaminated Galaxy Clusters Using Machine Learning
• The non-Gaussianity of the cosmic shear likelihood - or: How odd is the Chandra Deep Field South?
• Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection
• Large Covariance Matrices: Smooth Models from the 2-Point Correlation Function
• Towards Accurate Modelling of Galaxy Clustering on Small Scales: Testing the Standard $\Lambda \mathrm{CDM}$ + Halo Model
• On estimating cosmology-dependent covariance matrices
• A First Look at creating mock catalogs with machine learning techniques
• Imprint of DESI fiber assignment on the anisotropic power spectrum of emission line galaxies
• SDSS-III Baryon Oscillation Spectroscopic Survey: Analysis of Potential Systematics in Fitting of Baryon Acoustic Feature
• Accelerating Approximate Bayesian Computation with Quantile Regression: Application to Cosmological Redshift Distributions
• Kernels on Sample Sets via Nonparametric Divergence Estimates
• Including parameter dependence in the data and covariance for cosmological inference
• Mixture Models in Astronomy

Cited by in corpus (25)

• The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: measurement of the BAO and growth rate of structure of the luminous red galaxy sample from the anisotropic power spectrum between redshifts 0.6 and 1.0
• The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: measurement of the BAO and growth rate of structure of the emission line galaxy sample from the anisotropic power spectrum between redshift 0.6 and 1.1
• Probing Parity-Violation with the Four-Point Correlation Function of BOSS Galaxies
• Unified galaxy power spectrum measurements from 6dFGS, BOSS, and eBOSS
• Bayesian optimisation for likelihood-free cosmological inference
• Toward a robust inference method for the galaxy bispectrum: likelihood function and model selection
• Weak lensing scattering transform: dark energy and neutrino mass sensitivity
• Accelerated Bayesian SED Modeling using Amortized Neural Posterior Estimation
• Galaxy Clustering Analysis with SimBIG and the Wavelet Scattering Transform
• $\mathrm{S}\mathrm{IM}\mathrm{BIG}$: Mock Challenge for a Forward Modeling Approach to Galaxy Clustering
• The degeneracy between primordial non-Gaussianity and foregrounds in 21cm intensity mapping experiments
• A Parameter-Masked Mock Data Challenge for Beyond-Two-Point Galaxy Clustering Statistics
• Finding eV-scale Light Relics with Cosmological Observables
• EFTofLSS meets simulation-based inference: ${\sigma }_8$ from biased tracers
• Debiasing inference with approximate covariance matrices and other unidentified biases
• A First Detection of the Connected 4-Point Correlation Function of Galaxies Using the BOSS CMASS Sample
• Beyond Gaussian Noise: A Generalized Approach to Likelihood Analysis with non-Gaussian Noise
• A Bayesian method for combining theoretical and simulated covariance matrices for large-scale structure surveys
• A Composite Likelihood Approach for Inference under Photometric Redshift Uncertainty
• Quantification of high dimensional non-Gaussianities and its implication to Fisher analysis in cosmology
• Flow-Based Likelihoods for Non-Gaussian Inference
• Toward Accurate Modeling of Galaxy Clustering on Small Scales: Halo Model Extensions and Lingering Tension
• Constraining Star Formation Histories of Blue Galaxies using the Scatter between Stellar Mass and Halo Mass
• IQ Collaboratory III: The Empirical Dust Attenuation Framework -- Taking Hydrodynamical Simulations with a Grain of Dust
• Non-Gaussian likelihood of weak lensing power spectra