Constrained correlation functions

Free access article

Issue		A&A Volume 504, Number 3, September IV 2009


Page(s)		705 - 717
Section		Cosmology (including clusters of galaxies)
DOI		http://dx.doi.org/10.1051/0004-6361/200912424
Published online		16 July 2009

A&A 504, 705-717 (2009)
DOI: 10.1051/0004-6361/200912424

Constrained correlation functions

P. Schneider and J. Hartlap

Argelander-Institut für Astronomie, Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany
e-mail: [peter;hartlap]@astro.uni-bonn.de

Received 5 May 2009 / Accepted 7 July 2009

Abstract
Measurements of correlation functions and their comparison with theoretical models are often employed in natural sciences, including astrophysics and cosmology, to determine best-fitting model parameters and their confidence regions. Due to a lack of better descriptions, the likelihood function of the correlation function is often assumed to be a multi-variate Gaussian. Using different methods, we show that correlation functions have to satisfy constraint relations, owing to the non-negativity of the power spectrum of the underlying random process. Specifically, for any statistically homogeneous and (for more than one spatial dimension) isotropic random field with correlation function $\xi(x)$ , we derive inequalities for the correlation coefficients $r_n\equiv \xi(n x)/\xi(0)$ (for integer n) of the form $r_{n{\rm l}}\le r_n\le r_{n{\rm u}}$ , where the lower and upper bounds on r_n depend on the r_j, with j<n, or more explicitly $\Xi_{n-}\left\{ \xi(0),\xi(x), \xi(2x),\dots,\xi([n-1]x) \right\}\le \xi(n x) \le \Xi_{n+}\left\{ \xi(0),\xi(x), \xi(2x),\dots,\xi([n-1]x) \right\}.$ Explicit expressions for the bounds are obtained for arbitrary n. We show that these constraint equations very significantly limit the set of possible correlation functions. For one particular example of a fiducial cosmic shear survey, we show that the Gaussian likelihood ellipsoid has a significant spill-over into the region of correlation functions forbidden by the aforementioned constraints, rendering the resulting best-fitting model parameters and their error region questionable, and indicating the need for a better description of the likelihood function. We conduct some simple numerical experiments which explicitly demonstrate the failure of a Gaussian description for the likelihood of $\xi$ . Instead, the shape of the likelihood function of the correlation coefficients appears to follow approximately the shape of the bounds on the r_n, even if the Gaussian ellipsoid lies well within the allowed region. Therefore, we define a non-linear and coupled transformation of the r_n, based on these bounds. Some numerical experiments then indicate that a Gaussian is a much better description of the likelihood in these transformed variables than of the original correlation coefficients – in particular, the full probability distribution then lies explicitly in the allowed region. For more than one spatial dimension of the random field, the explicit expressions of the bounds on the r_n are not optimal. We outline a geometrical method how tighter bounds may be obtained in principle. We illustrate this method for a few simple cases; a more general treatment awaits future work.

Key words: methods: statistical -- cosmological parameters -- cosmology: miscellaneous

What is OpenURL?

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access. An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.

If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
You can define your own OpenURL resolver with your EDPS Account. In this case your choice will be given priority over that of your library.
You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see http://www.openly.com/openurlref/). You should disable this module if you wish to use the OpenURL server that you or your library have defined.