Tutorials › Statistical Theory I - Probability and Distribution Theory › Regular Variation and the Fréchet Tail Index

Reliability and Extreme-Value Laws · Tutorial 691 of 994

Regular Variation and the Fréchet Tail Index

Learn how a power-law survival tail yields a Fréchet limit for maxima and how its exponent determines the GEV shape parameter.

Intermediate 10 min read

What You'll Learn

  • Define regular variation for a survival tail and distinguish its exponent from its slowly varying factor.
  • Derive the limiting distribution of normalized maxima under a tail-based normalization.
  • Explain why the limit is Fréchet and identify its tail exponent.
  • Relate the regular-variation exponent to the positive GEV shape parameter.
  • Work with exact Pareto tails, shifted Pareto tails, and distributions with an atom at zero.

Why Tail Decay Determines Extreme-Value Shape

For a sample of independent observations, the maximum is usually far beyond the range of a typical observation. Its distribution therefore depends strongly on how quickly the survival probability decreases in the far tail. A power-law tail decreases more slowly than an exponential tail, and this slow decay produces a characteristic limit for normalized maxima: the Fréchet distribution.

Regular variation gives a precise way to describe this power-law behavior, including tails that are not exact powers. Its exponent also identifies the shape parameter of the corresponding generalized extreme-value (GEV) limit. The central calculation is that multiplying a large threshold by \(x\) multiplies the tail probability, asymptotically, by \(x^{-\alpha}\).

Regular Variation of a Survival Tail

Let \(X\) be a nonnegative random variable with distribution function \(F\). Its survival function is \(\bar F(x)=P(X>x)=1-F(x)\). The tail is regularly varying with index \(-\alpha\), for \(\alpha>0\), if it is eventually positive and

$$ \lim_{t\to\infty}\frac{\bar F(tx)}{\bar F(t)}=x^{-\alpha} \qquad\text{for every }x>0. $$

Equivalently, for sufficiently large \(t\), one can write \(\bar F(t)=t^{-\alpha}L(t)\), where \(L\) is a positive slowly varying function: \(L(tx)/L(t)\to1\) for each fixed \(x>0\). The factor \(L\) allows departures from an exact power law while preserving the same asymptotic tail exponent. Regular variation is a statement about ratios at large thresholds, not an assertion that the tail equals a power at every threshold.

Definition: The number \(\alpha>0\) is the tail exponent (or Fréchet tail index) when \(\bar F\) is regularly varying with index \(-\alpha\). A smaller \(\alpha\) corresponds to a heavier tail: at a fixed large threshold multiplier \(x>1\), the limiting tail ratio \(x^{-\alpha}\) is larger.

Worked Example: An Exact Pareto Tail

Suppose \(X\) has a Pareto distribution with scale \(c>0\) and exponent \(\alpha>0\), so that \(P(X>x)=(c/x)^\alpha\) for \(x\ge c\). For any fixed \(u>0\), once \(t\) is sufficiently large that both \(t\) and \(tu\) are at least \(c\),

$$ \frac{\bar F(tu)}{\bar F(t)} = \frac{(c/(tu))^\alpha}{(c/t)^\alpha} =u^{-\alpha}. $$

Thus the survival tail is regularly varying with index \(-\alpha\), with slowly varying factor \(L(t)=c^\alpha\), a constant. In particular, the tail exponent is \(\alpha\), regardless of the scale \(c\).

From Regular Variation to a Fréchet Limit

Let \(X_1,X_2,\ldots\) be independent and identically distributed nonnegative random variables with survival function \(\bar F\), regularly varying with index \(-\alpha\), where \(\alpha>0\). Write \(M_n=\max(X_1,\ldots,X_n)\). Choose positive constants \(a_n\) such that

$$ n\bar F(a_n)\longrightarrow 1. $$

Such a tail-based normalization is available for an unbounded regularly varying tail; one may use generalized tail quantiles. The condition states that a threshold of size \(a_n\) has exceedance probability approximately \(1/n\). Since \(\bar F(a_n)\to0\), this normalization also requires \(a_n\to\infty\).

Theorem: Under the conditions above, for every \(x>0\), $$ P(M_n/a_n\le x)\longrightarrow \exp(-x^{-\alpha}). $$ The limiting distribution, with distribution function \(G_\alpha(x)=\exp(-x^{-\alpha})\) for \(x>0\) and \(G_\alpha(x)=0\) for \(x\le0\), is the unit Fréchet distribution with tail exponent \(\alpha\).

To prove the result, fix \(x>0\). Regular variation, applied at \(t=a_n\), gives

$$ n\bar F(a_nx) = n\bar F(a_n)\frac{\bar F(a_nx)}{\bar F(a_n)} \longrightarrow 1\cdot x^{-\alpha} =x^{-\alpha}. $$

By independence, the probability that every observation is at most \(a_nx\) is

$$ P(M_n/a_n\le x) = \bigl(1-\bar F(a_nx)\bigr)^n. $$

Set \(q_n=\bar F(a_nx)\). The preceding limit shows \(nq_n\to x^{-\alpha}\), so \(q_n\to0\). For \(0\le q<1\), \(\log(1-q)/q\to-1\) as \(q\to0\). Consequently,

$$ n\log(1-q_n) =(nq_n)\frac{\log(1-q_n)}{q_n} \longrightarrow -x^{-\alpha}, $$

and exponentiating proves the claimed limit for \(x>0\). For \(x<0\), nonnegativity gives \(P(M_n/a_n\le x)=0\). At \(x=0\), the probability is \(P(X=0)^n\); regular variation ensures \(P(X>0)>0\), so this probability tends to zero. These values agree with \(G_\alpha\), completing the convergence in distribution. \(\square\)

Worked Example: Normalized Maxima from a Pareto Sample

Let \(X_1,\ldots,X_n\) be independent with Pareto survival \(\bar F(x)=(2/x)^3\) for \(x\ge2\), and let \(M_n\) be their maximum. Choose \(a_n=2n^{1/3}\). This threshold is at least 2, and

$$ n\bar F(a_n) =n\left(\frac{2}{2n^{1/3}}\right)^3 =n\cdot\frac{1}{n} =1. $$

For fixed \(x>0\), \(a_nx\ge2\) for all sufficiently large \(n\). Therefore

$$ P(M_n/a_n\le x) = \left[1-\left(\frac{2}{a_nx}\right)^3\right]^n = \left(1-\frac{x^{-3}}{n}\right)^n \longrightarrow \exp(-x^{-3}). $$

The limit has tail exponent \(3\). The scale \(2\) affects the normalization \(a_n\), but not the limiting exponent.

The Tail Exponent and the GEV Shape Parameter

In the usual GEV parameterization, the distribution function is

$$ H_{\xi,\mu,\sigma}(z) = \exp\left\{-\left[1+\xi\left(\frac{z-\mu}{\sigma}\right)\right]^{-1/\xi}\right\}, \qquad 1+\xi\left(\frac{z-\mu}{\sigma}\right)>0, $$

where \(\mu\) is a location parameter, \(\sigma>0\) is a scale parameter, and \(\xi\ne0\) is the shape parameter. For a Fréchet limit, \(\xi>0\). The unit Fréchet variable \(Y\) with \(P(Y\le y)=\exp(-y^{-\alpha})\) for \(y>0\) can be transformed by \(Z=(Y-1)/\xi\), where \(\xi=1/\alpha\). Then, on \(1+\xi z>0\),

$$ P(Z\le z) = P(Y\le1+\xi z) = \exp\left\{-(1+\xi z)^{-\alpha}\right\} = \exp\left\{-(1+\xi z)^{-1/\xi}\right\}. $$

This is a GEV distribution with \(\mu=0\), \(\sigma=1\), and \(\xi=1/\alpha\). The affine transformation changes location and scale conventions, not the shape. The same relationship is visible directly in the far tail: for \(\xi>0\), the GEV survival probability is asymptotic to a constant times \(z^{-1/\xi}\) as \(z\to\infty\). Matching that power with a regularly varying survival tail gives \(\alpha=1/\xi\).

Key relationship: A regularly varying survival tail with exponent \(\alpha>0\) leads to a Fréchet maximum limit with GEV shape parameter \(\xi=1/\alpha>0\). A heavy tail (small \(\alpha\)) therefore corresponds to a larger positive GEV shape parameter.

Worked Example: Reading the GEV Shape from a Tail

Suppose a model has survival tail \(\bar F(t)=t^{-2}L(t)\), where \(L\) is slowly varying. The regular-variation exponent is \(\alpha=2\). The corresponding Fréchet limit has distribution function \(\exp(-x^{-2})\) for \(x>0\), and its GEV shape parameter is

$$ \xi=\frac{1}{\alpha}=\frac12. $$

Conversely, a positive GEV shape value \(\xi=0.4\) corresponds to tail exponent \(\alpha=1/\xi=2.5\). These calculations identify the exponent and shape only; the tail's slowly varying factor and the distribution's scale determine the appropriate normalization of the maxima.

Atoms, Normalization, and a Common Pitfall

Regular variation concerns the tail at large positive thresholds. It does not require a continuous distribution, nor does it forbid an atom at zero. The maximum-limit proof uses the survival probability at \(a_nx\) for \(x>0\); the value at zero must be handled separately if one wants to state convergence there.

Worked Example: A Pareto Tail with an Atom at Zero

Let \(P(X=0)=0.2\), and with probability \(0.8\) let \(X\) have the Pareto survival \((1/x)^2\) for \(x\ge1\). For \(x\ge1\),

$$ \bar F(x)=0.8x^{-2}. $$

Thus the tail is regularly varying with exponent \(\alpha=2\). For sufficiently large \(n\), take \(a_n=(0.8n)^{1/2}\), which is at least 1. Then \(n\bar F(a_n)=1\), and the theorem gives the unit Fréchet limit \(\exp(-x^{-2})\) for \(x>0\). At zero, however, the finite-sample probability is

$$ P(M_n/a_n\le0)=P(X_1=0,\ldots,X_n=0)=0.2^n, $$

which is positive for every finite \(n\), but tends to zero. The limit is zero at \(x=0\), as required.

A common error is to claim that nonnegativity alone makes \(P(M_n/a_n\le0)=0\). Nonnegativity only implies that this event occurs when all observations equal zero; it does not rule out an atom there. The correction is to write its probability as \(P(X=0)^n\) and show that it tends to zero. Regular variation ensures \(P(X>0)>0\), so \(P(X=0)<1\).

Another important distinction is between the tail exponent and the normalization. The exponent \(\alpha\) determines the shape of the Fréchet limit, but it does not by itself determine the finite-sample scale \(a_n\). That scale depends on the full tail, including its slowly varying factor. A sound analysis therefore identifies both the tail ratio and a normalization satisfying \(n\bar F(a_n)\to1\).

Check Your Understanding

Use the definitions and limit argument above to answer the following questions.

  1. If \(\bar F(t)=t^{-4}L(t)\) with \(L\) slowly varying, what is the Fréchet tail exponent and the corresponding positive GEV shape parameter?
  2. Why does \(n\bar F(a_n)\to1\) imply \(a_n\to\infty\) for an unbounded regularly varying tail?
  3. For fixed \(x>0\), what limit does \(n\bar F(a_nx)\) have, and which definition gives this result?
  4. Why does the maximum-limit proof use independence when expressing \(P(M_n/a_n\le x)\) as a power?
  5. If \(P(X=0)=r\) with \(0\le r<1\), what is \(P(M_n/a_n\le0)\), and what is its limit?