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.

9 pages, 10 figures, accepted for publication in MNRAS

References in corpus (10)

• Array Programming with NumPy
• 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
• Cobaya: Code for Bayesian Analysis of hierarchical physical models
• 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
• Neutrino mass from Cosmology
• Planck constraints on the tensor-to-scalar ratio
• MCMC generation of cosmological fields far beyond Gaussianity