Abstract
We consider the shape of the posterior distribution to be used when fitting cosmological models to power spectra measured from galaxy surveys. At very large scales, Gaussian posterior distributions in the power do not approximate the posterior distribution P_R we expect for a Gaussian density field δ k, even if we vary the covariance matrix according to the model to be tested. We compare alternative posterior distributions with P_R, both mode-by-mode and in terms of expected measurements of primordial non-Gaussianity parametrized by fNL. Marginalising over a Gaussian posterior distribution P_f with fixed covariance matrix yields a posterior mean value of fNL which, for a data set with the characteristics of Euclid, will be underestimated by ▵fNL = 0.4, while for the data release 9 of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey (BOSS DR9; Ahn et al.) it will be underestimated by ▵fNL = 19.1. Adopting a different form of the posterior function means that we do not necessarily require a different covariance matrix for each model to be tested: this dependence is absorbed into the functional form of the posterior. Thus, the computational burden of analysis is significantly reduced.
Author
Kalus, B.; Percival, W. J.; Samushia, L.
Journal
Monthly Notices of the Royal Astronomical Society: Letters
Paper Publication Date
January 2016
Paper Type
Astrostatistics