Large deviations of Gaussian random functions

Large deviations of Gaussian random functions

A random function – of either one variable (a random process), or two or more variables(a random field) – is called Gaussian if every finite-dimensional distribution is a multivariate normal distribution. Gaussian random fields on the sphere are useful (for example) when analysing

* the anomalies in the cosmic microwave background radiation (see Robert J. Adler, "On excursion sets, tube formulas and maxima of random fields", [http://dx.doi.org/10.1214/aoap/1019737664 The Annals of Applied Probability 2000, Vol. 10, No. 1, 1-74] . (Special invited paper.)] , pp. 8-9);

* brain images obtained by positron emission tomography (see , pp. 9-10).

Sometimes, a value of a Gaussian random function deviates from its expected value by several standard deviations. This is a large deviation. Though rare in a small domain (of space or/and time), large deviations may be quite usual in a large domain.

The probability of such deviations is an old topic in probability theory. [cite journal
last = Cramér
first = Harald
title = Sur un nouveau théorème—limite de la théorie des probabilités
journal = Actualités Scientifiques et Industrielles
volume = 736
year = 1938
pages = 5–23
language = French
] Recent progressFact|date=August 2008 can be characterized as "good" (rather than "pure" or "applied") mathematics: it is more elegant than most "applied" mathematics, and at the same time more useful than some "pure" mathematics.

Basic statement

Let M be the maximal value of a Gaussian random function X on the(two-dimensional) sphere. Assume that the expected value of X is 0 (at every point of the sphere), and the standard deviation of X is 1 (at every point of the sphere). Then, for large a>0, P(M>a) is close to C a exp(-a^2/2) + 2P(xi>a),where xi is distributed N(0,1) (the standard normal distribution), and C is a constant; it does not depend on a, but depends on the correlation function of X (see below). The relative error of the approximation decays exponentially for large a.

The constant C is easy to determine in the important special case described in terms of the directional derivative of X at a given point (of the sphere) in a given direction (tangential to the sphere). The derivative is random, with zero expectation and some standard deviation. The latter may depend on the point and the direction. However, if it does not depend, then it is equal to (pi/2)^{1/4} C^{1/2} (for the sphere of radius 1).

The coefficient 2 before P(xi>a) is in fact the Euler characteristic of the sphere (for the torus it vanishes).

It is assumed that X is twice continuously differentiable (almost surely), and reaches its maximum at a single point (almost surely).

The clue: mean Euler characteristic

The clue to the theory sketched above is, Euler characteristic chi_a of the set {X>a} of all points t (of the sphere) such that X(t)>a. Its expected value (in other words, mean value) E(chi_a) can be calculated explicitly:

: E(chi_a) = C a exp(-a^2/2) + 2 P(xi>a)

(which is far from being trivial, and involves Poincare-Hopf theorem, Gauss-Bonnet theorem, Rice formula etc).

The set {X>a} is the empty set whenever M; in this case chi_a=0. In the other case, when M>a, the set {X>a} is non-empty; its Euler characteristic may take various values, depending on the topology of the set (the number of connected components, and possible holes in these components). However, if a is large and M>a then the set {X>a} is usually a small, slightly deformed disk or ellipse (which is easy to guess, but quite difficult to prove). Thus, its Euler characteristic chi_a is usually equal to 1 (given that M>a). This is why E(chi_a) is close to P(M>a).

See also

* Gaussian process
* Gaussian random field
* Large deviations theory

Further reading

The basic statement given above is a simple special case of a much more general (and difficult) theory stated by Adler [Robert J. Adler, Jonathan E. Taylor, "Random fields and geometry", Springer 2007. ISBN 978-0-387-48112-8] [Robert J. Adler, "Some new random field tools for spatial analysis", [http://arxiv.org/abs/0805.1031 arXiv:0805.1031] .] . For a detailed presentation of this special case see Tsirelson's lectures [ [http://www.tau.ac.il/~tsirel/Courses/Gauss2/syllabus.html Lectures of B. Tsirelson] (especially, Sect. 5).] .


Wikimedia Foundation. 2010.

Игры ⚽ Поможем написать курсовую

Look at other dictionaries:

  • List of mathematics articles (L) — NOTOC L L (complexity) L BFGS L² cohomology L function L game L notation L system L theory L Analyse des Infiniment Petits pour l Intelligence des Lignes Courbes L Hôpital s rule L(R) La Géométrie Labeled graph Labelled enumeration theorem Lack… …   Wikipedia

  • Normal distribution — This article is about the univariate normal distribution. For normally distributed vectors, see Multivariate normal distribution. Probability density function The red line is the standard normal distribution Cumulative distribution function …   Wikipedia

  • Riemann hypothesis — The real part (red) and imaginary part (blue) of the Riemann zeta function along the critical line Re(s) = 1/2. The first non trivial zeros can be seen at Im(s) = ±14.135, ±21.022 and ±25.011 …   Wikipedia

  • Nobel Prizes — ▪ 2009 Introduction Prize for Peace       The 2008 Nobel Prize for Peace was awarded to Martti Ahtisaari, former president (1994–2000) of Finland, for his work over more than 30 years in settling international disputes, many involving ethnic,… …   Universalium

  • List of numerical analysis topics — This is a list of numerical analysis topics, by Wikipedia page. Contents 1 General 2 Error 3 Elementary and special functions 4 Numerical linear algebra …   Wikipedia

  • List of statistics topics — Please add any Wikipedia articles related to statistics that are not already on this list.The Related changes link in the margin of this page (below search) leads to a list of the most recent changes to the articles listed below. To see the most… …   Wikipedia

  • Mean squared error — In statistics, the mean squared error (MSE) of an estimator is one of many ways to quantify the difference between values implied by a kernel density estimator and the true values of the quantity being estimated. MSE is a risk function,… …   Wikipedia

  • Errors and residuals in statistics — For other senses of the word residual , see Residual. In statistics and optimization, statistical errors and residuals are two closely related and easily confused measures of the deviation of a sample from its theoretical value . The error of a… …   Wikipedia

  • Matrix (mathematics) — Specific elements of a matrix are often denoted by a variable with two subscripts. For instance, a2,1 represents the element at the second row and first column of a matrix A. In mathematics, a matrix (plural matrices, or less commonly matrixes)… …   Wikipedia

  • Robust statistics — provides an alternative approach to classical statistical methods. The motivation is to produce estimators that are not unduly affected by small departures from model assumptions. Contents 1 Introduction 2 Examples of robust and non robust… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”