Likelihood non-Gaussianity in large-scale structure analyses. Issue 2 (26th February 2019)

• Record Type: Journal Article
• Authors: Hahn, ChangHoon; Beutler, Florian; Sinha, Manodeep; Berlind, Andreas; Ho, Shirley; Hogg, David W
• Is Part Of: Monthly notices of the Royal Astronomical Society. Volume 485:Issue 2(2019)
• Page Start: 2956
• Page End: 2969
• Publication Date: 2019-02-26
• Subjects: methods: data analysis -- methods: statistical -- galaxies: statistics -- cosmology: observations -- cosmological parameters -- large-scale structure of Universe; Astronomy -- Periodicals; Periodicals; 520.5
• DOI: 10.1093/mnras/stz558
• Languages: English
• ISSNs: 0035-8711

Abstract

ABSTRACT: Standard present-day large-scale structure (LSS) analyses make a major assumption in their Bayesian parameter inference – that the likelihood has a Gaussian form. For summary statistics currently used in LSS, this assumption, even if the underlying density field is Gaussian, cannot be correct in detail. We investigate the impact of this assumption on two recent LSS analyses: the Beutler et al. power spectrum multipole ( P ℓ ) analysis and the Sinha et al. group multiplicity function (ζ) analysis. Using non-parametric divergence estimators on mock catalogues originally constructed for covariance matrix estimation, we identify significant non-Gaussianity in both the P ℓ and ζ likelihoods. We then use Gaussian mixture density estimation and independent component analysis on the same mocks to construct likelihood estimates that approximate the true likelihood better than the Gaussian pseudo -likelihood. Using these likelihood estimates, we accurately estimate the true posterior probability distribution of the Beutler et al. and Sinha et al. parameters. Likelihood non-Gaussianity shifts the f σ8 constraint by −0.44σ, but otherwise does not significantly impact the overall parameter constraints of Beutler et al. For the ζ analysis, using the pseudo-likelihood significantly underestimates the uncertainties and biases the constraints of the Sinha et al. halo occupation parameters. For$\log \, M_1$ and α, the posteriors are shifted by +0.43σ and −0.51σ and broadened by$42{{\ \rm per\ cent}}$ and$66{{\ \rm per\ cent}}$, respectively. The divergence and likelihood estimation methods we present provide a straightforward framework for quantifying the impact of likelihood non-Gaussianity and deriving more accurate parameter constraints.