Parameter inference with estimated covariance matrices
read the original abstract
When inferring parameters from a Gaussian-distributed data set by computing a likelihood, a covariance matrix is needed that describes the data errors and their correlations. If the covariance matrix is not known a priori, it may be estimated and thereby becomes a random object with some intrinsic uncertainty itself. We show how to infer parameters in the presence of such an estimated covariance matrix, by marginalising over the true covariance matrix, conditioned on its estimated value. This leads to a likelihood function that is no longer Gaussian, but rather an adapted version of a multivariate t-distribution, which has the same numerical complexity as the multivariate Gaussian. As expected, marginalisation over the true covariance matrix improves inference when compared with Hartlap et al.'s method, which uses an unbiased estimate of the inverse covariance matrix but still assumes that the likelihood is Gaussian.
This paper has not been read by Pith yet.
Forward citations
Cited by 4 Pith papers
-
If at First You Don't Succeed, Trispectrum: I. Estimating the Matter Power Spectrum Covariance with Higher-Order Statistics
Estimators from squeezed bispectrum and collapsed trispectrum recover unbiased small-scale matter power spectrum covariance at the percent level using 25 Quijote simulations.
-
Probing the limits of cosmological information from the Lyman-$\alpha$ forest 2-point correlation functions
Using idealized synthetic data, knowing the true continuum in Lyα forest auto- and cross-correlations reduces uncertainties on the AP parameter and Ω_m by ~10%, with extension to 240 h^{-1}Mpc scales adding up to ~15%...
-
DESI 2024 V: Full-Shape Galaxy Clustering from Galaxies and Quasars
DESI DR1 full-shape galaxy clustering constrains Omega_m = 0.296 ± 0.010, H0 = 68.63 ± 0.79 km/s/Mpc, and sigma_8 = 0.841 ± 0.034, consistent with LambdaCDM and Planck.
-
Modeling and measuring the anisotropic halo 3-point correlation function: a coordinated study
Novel implementation of anisotropic 3PCF model and estimator tested on 298 halo catalogs shows degeneracy breaking between f and b1 in 3PCF-only analysis but limited added value in joint 2PCF+3PCF due to tree-level mo...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.