Bias on Tensor-to-Scalar Ratio Inference With Estimated Covariance Matrices

Abstract

We investigate simulation-based bandpower covariance matrices commonly used in cosmological parameter inferences such as the estimation of the tensor-to-scalar ratio $r$. We find that upper limits on $r$ can be biased low by tens of percent. The underestimation of the upper limit is most severe when the number of simulation realizations is similar to the number of observables. Convergence of the covariance-matrix estimation can require a number of simulations an order of magnitude larger than the number of observables, which could mean $\mathcal{O}(10000)$ simulations. This is found to be caused by an additional scatter in the posterior probability of $r$ due to Monte Carlo noise in the estimated bandpower covariance matrix, in particular, by spurious non-zero off-diagonal elements. We show that matrix conditioning can be a viable mitigation strategy in the case that legitimate covariance assumptions can be made.

arXiv comment9 pages, 10 figures, accepted for publication in MNRAS

References in corpus (28)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological analysis of the DR12 galaxy sample
• The Simons Observatory: Science goals and forecasts
• Early Dark Energy Can Resolve The Hubble Tension
• Signature of Gravity Waves in Polarization of the Microwave Background
• BICEP / Keck XIII: Improved Constraints on Primordial Gravitational Waves using Planck, WMAP, and BICEP/Keck Observations through the 2018 Observing Season
• Why your model parameter confidences might be too optimistic -- unbiased estimation of the inverse covariance matrix
• A Probe of Primordial Gravity Waves and Vorticity
• Planck 2018 results. VIII. Gravitational lensing
• Cobaya: Code for Bayesian Analysis of hierarchical physical models
• A unified pseudo-$C_{\ell }$ framework
• The Quest for B Modes from Inflationary Gravitational Waves
• Planck intermediate results. LVII. Joint Planck LFI and HFI data processing
• Measurements of the E-Mode Polarization and Temperature-E-Mode Correlation of the CMB from SPT-3G 2018 Data
• Improved limits on the tensor-to-scalar ratio using BICEP and Planck
• Likelihood Analysis of CMB Temperature and Polarization Power Spectra
• Parameter inference with estimated covariance matrices
• Neutrino mass from Cosmology
• The Effect of Covariance Estimator Error on Cosmological Parameter Constraints
• Putting the Precision in Precision Cosmology: How accurate should your data covariance matrix be?
• Planck constraints on the tensor-to-scalar ratio
• Measurements of B-mode Polarization of the Cosmic Microwave Background from 500 Square Degrees of SPTpol Data
• A Measurement of the Degree Scale CMB B-mode Angular Power Spectrum with POLARBEAR
• Large-scale CMB temperature and polarization cross-spectra likelihoods
• Cosmology with the CMB temperature-polarization correlation
• MCMC generation of cosmological fields far beyond Gaussianity