Data analysis in cosmology requires reliable covariance matrices. Covariance matrices derived from numerical simulations often require a very large number of realizations to be accurate. When a theoretical model for the covariance matrix exists, the parameters of the model can often be fit with many fewer simulations. We establish a rigorous Bayesian method for performing such a fit, but show that using the maximum posterior point is often sufficient. We show how a model covariance matrix can be tested by examining the appropriate χ2 distributions from simulations. We demonstrate our method on two examples. First, we measure the two-point correlation function of halos from a large set of 10000 mock halo catalogs. We build a model covariance with 2 free parameters, which we fit using our procedure. The resulting best-fit model covariance obtained from just 100 simulation realizations proves to be as reliable as the numerical covariance matrix built from the full 10000 set. We also test our method on a setup where the covariance matrix is large by measuring the halo bispectrum for thousands of triangles for the same set of mocks. We build a block diagonal model covariance with 2 free parameters as an improvement over the diagonal Gaussian covariance. Our model covariance passes the χ2 test only partially in this case, signaling that the model is insufficient even using free parameters, but significantly improves over the Gaussian one.

Fitting covariance matrix models to simulations

Alessandra Fumagalli
;
Alexandro Saro;Emiliano Sefusatti;Pierluigi Monaco;
2022-01-01

Abstract

Data analysis in cosmology requires reliable covariance matrices. Covariance matrices derived from numerical simulations often require a very large number of realizations to be accurate. When a theoretical model for the covariance matrix exists, the parameters of the model can often be fit with many fewer simulations. We establish a rigorous Bayesian method for performing such a fit, but show that using the maximum posterior point is often sufficient. We show how a model covariance matrix can be tested by examining the appropriate χ2 distributions from simulations. We demonstrate our method on two examples. First, we measure the two-point correlation function of halos from a large set of 10000 mock halo catalogs. We build a model covariance with 2 free parameters, which we fit using our procedure. The resulting best-fit model covariance obtained from just 100 simulation realizations proves to be as reliable as the numerical covariance matrix built from the full 10000 set. We also test our method on a setup where the covariance matrix is large by measuring the halo bispectrum for thousands of triangles for the same set of mocks. We build a block diagonal model covariance with 2 free parameters as an improvement over the diagonal Gaussian covariance. Our model covariance passes the χ2 test only partially in this case, signaling that the model is insufficient even using free parameters, but significantly improves over the Gaussian one.
File in questo prodotto:
File Dimensione Formato  
3021243(1)_organized.pdf

Accesso chiuso

Tipologia: Documento in Versione Editoriale
Licenza: Copyright Editore
Dimensione 9.14 MB
Formato Adobe PDF
9.14 MB Adobe PDF   Visualizza/Apri   Richiedi una copia
2206.05191.pdf

accesso aperto

Tipologia: Altro materiale allegato
Licenza: Creative commons
Dimensione 1.37 MB
Formato Adobe PDF
1.37 MB Adobe PDF Visualizza/Apri
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11368/3044402
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 2
  • ???jsp.display-item.citation.isi??? 3
social impact