Precision matrix expansion - efficient use of numerical simulations in estimating errors on cosmological parameters

Abstract

Computing the inverse covariance matrix (or precision matrix) of large data vectors is crucial in weak lensing (and multi-probe) analyses of the large scale structure of the universe. Analytically computed covariances are noise-free and hence straightforward to invert, however the model approximations might be insufficient for the statistical precision of future cosmological data. Estimating covariances from numerical simulations improves on these approximations, but the sample covariance estimator is inherently noisy, which introduces uncertainties in the error bars on cosmological parameters and also additional scatter in their best fit values. For future surveys, reducing both effects to an acceptable level requires an unfeasibly large number of simulations. In this paper we describe a way to expand the true precision matrix around a covariance model and show how to estimate the leading order terms of this expansion from simulations. This is especially powerful if the covariance matrix is the sum of two contributions, $\mathbf{C}=\mathbf{A}+\mathbf{B}$, where $\mathbf{A}$ is well understood analytically and can be turned off in simulations (e.g. shape-noise for cosmic shear) to yield a direct estimate of $\mathbf{B}$. We test our method in mock experiments resembling tomographic weak lensing data vectors from the Dark Energy Survey (DES) and the Large Synoptic Survey Telecope (LSST). For DES we find that $400$ N-body simulations are sufficient to achive negligible statistical uncertainties on parameter constraints. For LSST this is achieved with $2400$ simulations. The standard covariance estimator would require >$10^5$ simulations to reach a similar precision. We extend our analysis to a DES multi-probe case finding a similar performance.

14 pages, submitted to mnras

References in corpus (4)

• Why your model parameter confidences might be too optimistic -- unbiased estimation of the inverse covariance matrix
• Shrinkage Estimation of the Power Spectrum Covariance Matrix
• Quantifying lost information due to covariance matrix estimation in parameter inference
• Non-linear shrinkage estimation of large-scale structure covariance

Cited by in corpus (29)

• Weak lensing for precision cosmology
• Cosmological constraints from BOSS with analytic covariance matrices
• Dark Energy Survey Year 3 Results: Covariance Modelling and its Impact on Parameter Estimation and Quality of Fit
• 2D-FFTLog: Efficient computation of real space covariance matrices for galaxy clustering and weak lensing
• Fewer Mocks and Less Noise: Reducing the Dimensionality of Cosmological Observables with Subspace Projections
• Nuw CDM cosmology from the weak lensing convergence PDF
• Dark Energy Survey Year 3 results: cosmology with moments of weak lensing mass maps -- validation on simulations
• Massive data compression for parameter-dependent covariance matrices
• Cosmic Shear Covariance Matrix in $w$CDM: Cosmology Matters
• Primordial non-Gaussianity without tails -- how to measure fNL with the bulk of the density PDF
• Mock galaxy shape catalogs in the Subaru Hyper Suprime-Cam Survey
• Cosmology with the Wide-Field Infrared Survey Telescope -- Synergies with the Rubin Observatory Legacy Survey of Space and Time
• Dark Energy Survey Year 1 Results: galaxy mock catalogues for BAO
• The integrated 3-point correlation function of cosmic shear
• RascalC: A Jackknife Approach to Estimating Single and Multi-Tracer Galaxy Covariance Matrices
• BAO from angular clustering: optimization and mitigation of theoretical systematics
• Cosmological inference from emulator based halo model I: Validation tests with HSC and SDSS mock catalogs
• Non-Gaussianity in the Weak Lensing Correlation Function Likelihood -- Implications for Cosmological Parameter Biases
• A Bayesian method for combining theoretical and simulated covariance matrices for large-scale structure surveys
• RedMaPPer: Evolution and Mass Dependence of the Conditional Luminosity Functions of Red Galaxies in Galaxy Clusters
• It takes two to know one: Computing accurate one-point PDF covariances from effective two-point PDF models
• Modal compression of the redshift-space galaxy bispectrum
• Bayesian Control Variates for optimal covariance estimation with pairs of simulations and surrogates
• Covariances of density probability distribution functions. Lessons from hierarchical models
• Simulation-based inference has its own Dodelson-Schneider effect (but it knows that it does)
• Euclid: Fast two-point correlation function covariance through linear construction
• A comparison of shrinkage estimators of the cosmological precision matrix
• Super sample covariance and the volume scaling of galaxy survey covariance matrices
• Matching Bayesian and frequentist coverage probabilities when using an approximate data covariance matrix