Rate function

From HandWiki
Revision as of 15:53, 6 February 2024 by ScienceGen (talk | contribs) (link)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Short description: Probability function

In mathematics — specifically, in large deviations theory — a rate function is a function used to quantify the probabilities of rare events. Such functions are used to formulate large deviation principles. A large deviation principle quantifies the asymptotic probability of rare events for a sequence of probabilities.

A rate function is also called a Cramér function, after the Swedish probabilist Harald Cramér.

Definitions

Rate function An extended real-valued function [math]\displaystyle{ I: X \to [0, +\infty] }[/math] defined on a Hausdorff topological space [math]\displaystyle{ X }[/math] is said to be a rate function if it is not identically [math]\displaystyle{ +\infty }[/math] and is lower semi-continuous i.e. all the sub-level sets

[math]\displaystyle{ \{ x \in X \mid I(x) \leq c \} \mbox{ for } c \geq 0 }[/math]

are closed in [math]\displaystyle{ X }[/math]. If, furthermore, they are compact, then [math]\displaystyle{ I }[/math] is said to be a good rate function.

A family of probability measures [math]\displaystyle{ (\mu_{\delta})_{\delta \gt 0} }[/math] on [math]\displaystyle{ X }[/math] is said to satisfy the large deviation principle with rate function [math]\displaystyle{ I: X \to [0, +\infty) }[/math] (and rate [math]\displaystyle{ 1/\delta }[/math]) if, for every closed set [math]\displaystyle{ F \subseteq X }[/math] and every open set [math]\displaystyle{ G \subseteq X }[/math],

[math]\displaystyle{ \limsup_{\delta \downarrow 0} \delta \log \mu_{\delta} (F) \leq - \inf_{x \in F} I(x), \quad \mbox{(U)} }[/math]
[math]\displaystyle{ \liminf_{\delta \downarrow 0} \delta \log \mu_{\delta} (G) \geq - \inf_{x \in G} I(x). \quad \mbox{(L)} }[/math]

If the upper bound (U) holds only for compact (instead of closed) sets [math]\displaystyle{ F }[/math], then [math]\displaystyle{ (\mu_{\delta})_{\delta \gt 0} }[/math] is said to satisfy the weak large deviations principle (with rate [math]\displaystyle{ 1/\delta }[/math] and weak rate function [math]\displaystyle{ I }[/math]).

Remarks

The role of the open and closed sets in the large deviation principle is similar to their role in the weak convergence of probability measures: recall that [math]\displaystyle{ (\mu_{\delta})_{\delta \gt 0} }[/math] is said to converge weakly to [math]\displaystyle{ \mu }[/math] if, for every closed set [math]\displaystyle{ F \subseteq X }[/math] and every open set [math]\displaystyle{ G \subseteq X }[/math],

[math]\displaystyle{ \limsup_{\delta \downarrow 0} \mu_\delta (F) \leq \mu(F), }[/math]
[math]\displaystyle{ \liminf_{\delta \downarrow 0} \mu_\delta (G) \geq \mu(G). }[/math]

There is some variation in the nomenclature used in the literature: for example, den Hollander (2000) uses simply "rate function" where this article — following Dembo & Zeitouni (1998) — uses "good rate function", and "weak rate function". Rassoul-Agha & Seppäläinen (2015) uses the term "tight rate function" instead of "good rate function" due to the connection with exponential tightness of a family of measures. Regardless of the nomenclature used for rate functions, examination of whether the upper bound inequality (U) is supposed to hold for closed or compact sets tells one whether the large deviation principle in use is strong or weak.

Properties

Uniqueness

A natural question to ask, given the somewhat abstract setting of the general framework above, is whether the rate function is unique. This turns out to be the case: given a sequence of probability measures (μδ)δ>0 on X satisfying the large deviation principle for two rate functions I and J, it follows that I(x) = J(x) for all x ∈ X.

Exponential tightness

It is possible to convert a weak large deviation principle into a strong one if the measures converge sufficiently quickly. If the upper bound holds for compact sets F and the sequence of measures (μδ)δ>0 is exponentially tight, then the upper bound also holds for closed sets F. In other words, exponential tightness enables one to convert a weak large deviation principle into a strong one.

Continuity

Naïvely, one might try to replace the two inequalities (U) and (L) by the single requirement that, for all Borel sets S ⊆ X,

[math]\displaystyle{ \lim_{\delta \downarrow 0} \delta \log \mu_{\delta} (S) = - \inf_{x \in S} I(x). \quad \mbox{(E)} }[/math]

The equality (E) is far too restrictive, since many interesting examples satisfy (U) and (L) but not (E). For example, the measure μδ might be non-atomic for all δ, so the equality (E) could hold for S = {x} only if I were identically +∞, which is not permitted in the definition. However, the inequalities (U) and (L) do imply the equality (E) for so-called I-continuous sets S ⊆ X, those for which

[math]\displaystyle{ I \big( \stackrel{\circ}{S} \big) = I \big( \bar{S} \big), }[/math]

where [math]\displaystyle{ \stackrel{\circ}{S} }[/math] and [math]\displaystyle{ \bar{S} }[/math] denote the interior and closure of S in X respectively. In many examples, many sets/events of interest are I-continuous. For example, if I is a continuous function, then all sets S such that

[math]\displaystyle{ S \subseteq \bar{\stackrel{\circ}{S}} }[/math]

are I-continuous; all open sets, for example, satisfy this containment.

Transformation of large deviation principles

Given a large deviation principle on one space, it is often of interest to be able to construct a large deviation principle on another space. There are several results in this area:

History and basic development

The notion of a rate function emerged in the 1930s with the Swedish mathematician Harald Cramér's study of a sequence of i.i.d. random variables (Zi)i∈[math]\displaystyle{ \mathbb{N} }[/math]. Namely, among some considerations of scaling, Cramér studied the behavior of the distribution of the average [math]\displaystyle{ X_n = \frac 1 n \sum_{i=1}^n Z_i }[/math] as n→∞.[1] He found that the tails of the distribution of Xn decay exponentially as e(x) where the factor λ(x) in the exponent is the Legendre–Fenchel transform (a.k.a. the convex conjugate) of the cumulant-generating function [math]\displaystyle{ \Psi_Z(t)=\log \operatorname E e^{tZ}. }[/math] For this reason this particular function λ(x) is sometimes called the Cramér function. The rate function defined above in this article is a broad generalization of this notion of Cramér's, defined more abstractly on a probability space, rather than the state space of a random variable.

See also


References

  1. Cramér, Harald (1938). "Sur un nouveau théorème-limite de la théorie des probabilités" (in fr). Colloque consacré à la théorie des probabilités, Part 3, Actualités scientifiques et industrielles 731: 5–23. 
  • Dembo, Amir; Zeitouni, Ofer (1998). Large deviations techniques and applications. Applications of Mathematics (New York) 38 (Second ed.). New York: Springer-Verlag. p. xvi+396. ISBN 0-387-98406-2.  MR1619036
  • den Hollander, Frank (2000). Large deviations. Fields Institute Monographs 14. Providence, RI: American Mathematical Society. p. x+143. ISBN 0-8218-1989-5.  MR1739680
  • Rassoul-Agha, Firas; Seppäläinen, Timo (2015). A course on large deviations with an introduction to Gibbs measures. Graduate Studies in Mathematics 162. Providence, RI: American Mathematical Society. p. xiv+318. ISBN 978-0-8218-7578-0.  MR3309619