Creating accurate and low-noise covariance matrices represents a formidable challenge in modem-day cosmology. We present a formalism to compress arbitrary observables into a small number of bins by projection into a model-specific subspace that minimizes the prior-averaged log-likelihood error. The lower dimensionality leads to a dramatic reduction in covariance matrix noise, significantly reducing the number of mocks that need to be computed. Given a theory model, a set of priors, and a simple model of the covariance, our method works by using singular value decompositions to construct a basis for the observable that is close to Euclidean; by restricting to the first few basis vectors, we can capture almost all the constraining power in a lower-dimensional subspace. Unlike conventional approaches, the method can be tailored for specific analyses and captures nonlinearities that are not present in the Fisher matrix, ensuring that the full likelihood can be reproduced. The procedure is validated with full-shape analyses of power spectra from Baryon Oscillation Spectroscopic Survey (BOSS) DR12 mock catalogs, showing that the 96-bin power spectra can be replaced by 12 subspace coefficients without biasing the output cosmology; this allows for accurate parameter inference using only similar to 100 mocks. Such decompositions facilitate accurate testing of power spectrum covariances; for the largest BOSS data chunk, we find the following: (a) analytic covariances provide accurate models (with or without trispectrum terms); and (b) using the sample covariance from the MultiDark-Patchy mocks incurs a similar to 0.5 sigma shift in Omega(m) unless the subspace projection is applied. The method is easily extended to higher order statistics; the similar to 2000-bin bispectrum can be compressed into only similar to 10 coefficients, allowing for accurate analyses using few mocks and without having to increase the bin sizes.

Fewer mocks and less noise: Reducing the dimensionality of cosmological observables with subspace projections / Oliver H. E. Philcox; Mikhail M. Ivanov; Matias Zaldarriaga; Marko Simonovic; Marcel Schmittfull. - In: PHYSICAL REVIEW D. - ISSN 2470-0010. - STAMPA. - 103:(2021), pp. 43508-43527. [10.1103/physrevd.103.043508]

Fewer mocks and less noise: Reducing the dimensionality of cosmological observables with subspace projections

Marko Simonovic;
2021

Abstract

Creating accurate and low-noise covariance matrices represents a formidable challenge in modem-day cosmology. We present a formalism to compress arbitrary observables into a small number of bins by projection into a model-specific subspace that minimizes the prior-averaged log-likelihood error. The lower dimensionality leads to a dramatic reduction in covariance matrix noise, significantly reducing the number of mocks that need to be computed. Given a theory model, a set of priors, and a simple model of the covariance, our method works by using singular value decompositions to construct a basis for the observable that is close to Euclidean; by restricting to the first few basis vectors, we can capture almost all the constraining power in a lower-dimensional subspace. Unlike conventional approaches, the method can be tailored for specific analyses and captures nonlinearities that are not present in the Fisher matrix, ensuring that the full likelihood can be reproduced. The procedure is validated with full-shape analyses of power spectra from Baryon Oscillation Spectroscopic Survey (BOSS) DR12 mock catalogs, showing that the 96-bin power spectra can be replaced by 12 subspace coefficients without biasing the output cosmology; this allows for accurate parameter inference using only similar to 100 mocks. Such decompositions facilitate accurate testing of power spectrum covariances; for the largest BOSS data chunk, we find the following: (a) analytic covariances provide accurate models (with or without trispectrum terms); and (b) using the sample covariance from the MultiDark-Patchy mocks incurs a similar to 0.5 sigma shift in Omega(m) unless the subspace projection is applied. The method is easily extended to higher order statistics; the similar to 2000-bin bispectrum can be compressed into only similar to 10 coefficients, allowing for accurate analyses using few mocks and without having to increase the bin sizes.
2021
103
43508
43527
Oliver H. E. Philcox; Mikhail M. Ivanov; Matias Zaldarriaga; Marko Simonovic; Marcel Schmittfull
File in questo prodotto:
File Dimensione Formato  
PhysRevD103043508.pdf

accesso aperto

Tipologia: Pdf editoriale (Version of record)
Licenza: Open Access
Dimensione 1.13 MB
Formato Adobe PDF
1.13 MB Adobe PDF

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

Utilizza questo identificatore per citare o creare un link a questa risorsa: https://hdl.handle.net/2158/1345820
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 32
  • ???jsp.display-item.citation.isi??? 32
social impact