On the insufficiency of arbitrarily precise covariance matrices: non-Gaussian weak lensing likelihoods

Abstract

We investigate whether a Gaussian likelihood, as routinely assumed in the analysis of cosmological data, is supported by simulated survey data. We define test statistics, based on a novel method that first destroys Gaussian correlations in a dataset, and then measures the non-Gaussian correlations that remain. This procedure flags pairs of datapoints which depend on each other in a non-Gaussian fashion, and thereby identifies where the assumption of a Gaussian likelihood breaks down. Using this diagnostic, we find that non-Gaussian correlations in the CFHTLenS cosmic shear correlation functions are significant. With a simple exclusion of the most contaminated datapoints, the posterior for $s_8$ is shifted without broadening, but we find no significant reduction in the tension with $s_8$ derived from Planck Cosmic Microwave Background data. However, we also show that the one-point distributions of the correlation statistics are noticeably skewed, such that sound weak lensing data sets are intrinsically likely to lead to a systematically low lensing amplitude being inferred. The detected non-Gaussianities get larger with increasing angular scale such that for future wide-angle surveys such as Euclid or LSST, with their very small statistical errors, the large-scale modes are expected to be increasingly affected. The shifts in posteriors may then not be negligible and we recommend that these diagnostic tests be run as part of future analyses.

Replacement to match accepted MNRAS version

References in corpus (2)

• Why your model parameter confidences might be too optimistic -- unbiased estimation of the inverse covariance matrix
• Quantifying lost information due to covariance matrix estimation in parameter inference

Cited by in corpus (22)

• Scientific Synergy Between LSST and Euclid
• Cosmic shear with small scales: DES-Y3, KiDS-1000 and HSC-DR1
• Quijote-PNG: Simulations of primordial non-Gaussianity and the information content of the matter field power spectrum and bispectrum
• Beware of commonly used approximations II: estimating systematic biases in the best-fit parameters
• Dark Energy Survey internal consistency tests of the joint cosmological probes analysis with posterior predictive distributions
• KiDS-SBI: Simulation-based inference analysis of KiDS-1000 cosmic shear
• EFTofLSS meets simulation-based inference: ${\sigma }_8$ from biased tracers
• Probing massive neutrinos with the Minkowski functionals of the galaxy distribution
• Sufficiency of a Gaussian power spectrum likelihood for accurate cosmology from upcoming weak lensing surveys
• A simulation-based inference pipeline for cosmic shear with the Kilo-Degree Survey
• Almanac: Weak Lensing power spectra and map inference on the masked sphere
• Quantification of high dimensional non-Gaussianities and its implication to Fisher analysis in cosmology
• Testing the accuracy of likelihoods for cluster abundance cosmology
• Bayesian Control Variates for optimal covariance estimation with pairs of simulations and surrogates
• KiDS-Legacy: Covariance validation and the unified OneCovariance framework for projected large-scale structure observables
• The impact of signal-to-noise, redshift, and angular range on the bias of weak lensing 2-point functions
• Hybrid summary statistics: neural weak lensing inference beyond the power spectrum
• Weak Gravitational Lensing
• A comparison of shrinkage estimators of the cosmological precision matrix
• Galaxy-Point Spread Function correlations as a probe of weak-lensing systematics with UNIONS data
• Non-Gaussian likelihood of weak lensing power spectra
• The Kaiser-Rocket effect: three decades and counting