A quasi-Gaussian approximation for the probability distribution of correlation functions

Abstract

Context. Whenever correlation functions are used for inference about cosmological parameters in the context of a Bayesian analysis, the likelihood function of correlation functions needs to be known. Usually, it is approximated as a multivariate Gaussian, though this is not necessarily a good approximation. Aims. We show how to calculate a better approximation for the probability distribution of correlation functions, which we call "quasi-Gaussian". Methods. Using the exact univariate PDF as well as constraints on correlation functions previously derived, we transform the correlation functions to an unconstrained variable for which the Gaussian approximation is well justified. From this Gaussian in the transformed space, we obtain the quasi-Gaussian PDF. The two approximations for the probability distributions are compared to the "true" distribution as obtained from simulations. Additionally, we test how the new approximation performs when used as likelihood in a toy-model Bayesian analysis. Results. The quasi-Gaussian PDF agrees very well with the PDF obtained from simulations; in particular, it provides a significantly better description than a straightforward copula approach. In a simple toy-model likelihood analysis, it yields noticeably different results than the Gaussian likelihood, indicating its possible impact on cosmological parameter estimation.

16 pages, 14 figures, updated to match version accepted by A&A

References in corpus (5)

• On the assumption of Gaussianity for cosmological two-point statistics and parameter dependent covariance matrices
• Copula Cosmology: Constructing a Likelihood Function
• Precise Estimation of Cosmological Parameters Using a More Accurate Likelihood Function
• Effect of model-dependent covariance matrix for studying Baryon Acoustic Oscillations
• Likelihood of the Power Spectrum in Cosmological Parameter Estimation

Cited by in corpus (8)

• Cosmology with cosmic shear observations: a review
• Non-Gaussianity in the Weak Lensing Correlation Function Likelihood -- Implications for Cosmological Parameter Biases
• CFHTLenS: A Gaussian likelihood is a sufficient approximation for a cosmological analysis of third-order cosmic shear statistics
• Cosmological Inference from Galaxy-Clustering Power Spectrum: Gaussianization and Covariance Decomposition
• Sufficiency of a Gaussian power spectrum likelihood for accurate cosmology from upcoming weak lensing surveys
• Likelihood of the Power Spectrum in Cosmological Parameter Estimation
• Constrained correlation functions from the Millennium Simulation
• Non-Gaussian likelihood of weak lensing power spectra