Tutorials › Statistical Theory I - Probability and Distribution Theory › Normalizing Maxima in Standard Distribution Families

Reliability and Extreme-Value Laws · Tutorial 685 of 994

Normalizing Maxima in Standard Distribution Families

Learn to choose and verify maximum normalizations directly from the distribution of one observation.

Intermediate 9 min read

What You'll Learn

  • Obtain the exact CDF of a maximum from the common CDF of independent observations
  • Choose centering and scaling constants for exponential samples
  • Scale Pareto maxima using the sample size and tail index
  • Normalize uniform maxima near their finite upper endpoint
  • Verify each proposed limit by evaluating the normalized maximum CDF
  • Distinguish a reverse-Weibull normalization from the equivalent endpoint-gap limit

Why Normalize a Maximum?

For independent observations, the maximum usually grows with the sample size or moves closer to a finite upper endpoint. A useful limit therefore requires more than taking the maximum itself: we must choose a centering constant \(b_n\) and a positive scale \(a_n\), then study \((M_n-b_n)/a_n\). The choices depend on the tail of the individual observations.

This tutorial works out those choices for exponential, Pareto, and uniform samples. In each case, the verification comes from the exact CDF of the maximum, not from an informal argument about what large observations are likely to do.

Definition: Let \(X_1,\ldots,X_n\) be independent and identically distributed random variables with CDF \(F\), and define \(M_n=\max(X_1,\ldots,X_n)\). The constants \(b_n\) and \(a_n>0\) are called a centering and a scaling sequence when the normalized maximum \((M_n-b_n)/a_n\) has a nondegenerate limiting distribution as \(n\) grows.

The Maximum CDF and a Useful Limit

For any real \(t\), the event \(M_n\leq t\) occurs exactly when every observation is at most \(t\). Independence therefore gives the basic identity.

$$ P(M_n\leq t) =P(X_1\leq t,\ldots,X_n\leq t) =F(t)^n. $$

Thus, to study a normalized maximum at a fixed value \(x\), evaluate \(F(b_n+a_nx)^n\). The following tail calculation will be used repeatedly.

Theorem: Suppose \(u_n\in[0,1]\) and \(n u_n\to \tau\), where \(0\leq \tau<\infty\). Then \((1-u_n)^n\to e^{-\tau}\).

Proof. If \(\tau=0\), then \(0\leq (1-u_n)^n\leq 1\), and Bernoulli’s inequality gives \((1-u_n)^n\geq 1-nu_n\) for \(u_n\in[0,1]\). Both bounds tend to 1. If \(\tau>0\), then \(u_n\to0\), and for all sufficiently large \(n\), \(u_n<1\). Taking logarithms,

$$ n\log(1-u_n) =nu_n\frac{\log(1-u_n)}{u_n} \longrightarrow \tau(-1)=-\tau, $$

because \(\log(1-u)/u\to-1\) as \(u\downarrow0\). Exponentiating proves the result. \(\square\)

In an extreme-value calculation, \(u_n\) is often the probability that one observation exceeds the threshold \(b_n+a_nx\). If that probability is approximately \(\tau(x)/n\), the maximum CDF is approximately \((1-\tau(x)/n)^n\), which tends to \(e^{-\tau(x)}\).

Exponential Samples: Center by the Logarithmic Growth

Let each observation have an exponential distribution with rate \(\lambda>0\), so that \(F(t)=1-e^{-\lambda t}\) for \(t\geq0\). The upper-tail probability at \(t\) is \(e^{-\lambda t}\). To make the expected number of observations above a threshold roughly one, set \(n e^{-\lambda t}\) near 1. Solving gives the centering \(b_n=(\log n)/\lambda\). A shift of size \(1/\lambda\) changes the tail by a fixed multiplicative factor, suggesting \(a_n=1/\lambda\).

Worked Example: Exponential Maxima

Suppose \(X_1,\ldots,X_n\) are independent exponential lifetimes with rate \(\lambda=2\). Find a centering and scale, and verify the limiting CDF.

Take \(b_n=(\log n)/2\) and \(a_n=1/2\). For fixed \(x\), the threshold is \(b_n+a_nx=(\log n+x)/2\), which is nonnegative for all sufficiently large \(n\). The one-observation survival probability there is

$$ e^{-2(b_n+a_nx)} =e^{-(\log n+x)} =\frac{e^{-x}}{n}. $$

Using the exact maximum CDF,

$$ P\left(\frac{M_n-b_n}{a_n}\leq x\right) =\left(1-\frac{e^{-x}}{n}\right)^n \longrightarrow \exp(-e^{-x}). $$

For a finite-sample check, take \(n=100\) and \(x=1\). The threshold is \((\log 100+1)/2\), the survival probability for one observation is \(e^{-1}/100\), and the exact CDF is \((1-e^{-1}/100)^{100}\), approximately \(0.692\). The limiting value is \(\exp(-e^{-1})\), also approximately \(0.692\).

The resulting limit is the standard Gumbel CDF \(G(x)=\exp(-e^{-x})\). The rate \(\lambda\) affects the centering and scale, but not this limiting CDF after normalization.

Pareto Samples: Scale by the Tail Index

Consider a Pareto distribution with minimum value \(x_m>0\) and shape parameter \(\alpha>0\). Its CDF is \(F(t)=1-(x_m/t)^\alpha\) for \(t\geq x_m\). The tail probability decreases as a power rather than exponentially. At a threshold \(t=x_m n^{1/\alpha}x\), for fixed \(x>0\) and sufficiently large \(n\), the exceedance probability is \(x^{-\alpha}/n\). This suggests scaling the maximum by \(x_m n^{1/\alpha}\), with no centering.

Worked Example: Pareto Maxima

Let \(X_1,\ldots,X_n\) have a Pareto distribution with \(x_m=2\) and \(\alpha=3\). Set \(W_n=M_n/(2n^{1/3})\). For fixed \(x>0\), the threshold \(2n^{1/3}x\) is at least 2 for all sufficiently large \(n\), and the survival probability of one observation at that threshold is

$$ \left(\frac{2}{2n^{1/3}x}\right)^3 =\frac{x^{-3}}{n}. $$

It follows that

$$ P(W_n\leq x) =\left(1-\frac{x^{-3}}{n}\right)^n \longrightarrow \exp(-x^{-3}). $$

For \(x\leq0\), \(P(W_n\leq x)=0\), since \(M_n>0\). Hence the limiting CDF is zero for \(x\leq0\) and \(\exp(-x^{-3})\) for \(x>0\), the Fréchet limit with shape \(\alpha=3\). As a finite-sample check, at \(n=100\) and \(x=2\), the one-observation exceedance probability is \(2^{-3}/100=1/800\). The exact CDF is \((799/800)^{100}\), approximately \(0.882\), matching the limit \(\exp(-1/8)\), approximately \(0.882\).

More generally, for any \(x_m>0\) and \(\alpha>0\), the same calculation gives the limit \(\exp(-x^{-\alpha})\) for \(x>0\). The power \(n^{1/\alpha}\) is essential: it balances the Pareto tail so that the probability of exceeding the normalized threshold is of order \(1/n\).

Uniform Samples: Center at the Upper Endpoint

If \(X_i\sim\operatorname{Uniform}(0,1)\), the maximum cannot grow without bound. Instead, it approaches the upper endpoint 1. The natural scale is the distance from that endpoint: thresholds of the form \(1+x/n\) have CDF \(1+x/n\) when \(x\leq0\) and \(1+x/n\geq0\). This leads to centering \(b_n=1\) and scale \(a_n=1/n\) for \((M_n-b_n)/a_n\).

Worked Example: Uniform Maxima

For independent \(X_i\sim\operatorname{Uniform}(0,1)\), define \(Z_n=n(M_n-1)\). If \(x<0\), then for all sufficiently large \(n\), \(0\leq1+x/n\leq1\), so

$$ P(Z_n\leq x) =P(M_n\leq1+x/n) =(1+x/n)^n \longrightarrow e^x. $$

For \(x\geq0\), \(1+x/n\geq1\), so \(P(Z_n\leq x)=1\). Thus the limiting CDF is \(H(x)=e^x\) for \(x<0\) and \(H(x)=1\) for \(x\geq0\), a reverse-Weibull limit with finite upper endpoint 0.

For a numerical check, take \(n=20\) and \(x=-0.5\). The exact CDF is \((1-0.5/20)^{20}=0.975^{20}\), approximately \(0.603\). The limit is \(e^{-0.5}\), approximately \(0.607\).

The endpoint gap has a complementary normalization. Let \(Y_n=n(1-M_n)=-Z_n\). For \(y\geq0\) and \(n\geq y\),

$$ P(Y_n\leq y) =P(M_n\geq1-y/n) =1-(1-y/n)^n \longrightarrow 1-e^{-y}. $$

For \(y<0\), this probability is zero. Therefore the endpoint gap converges to an exponential distribution with rate 1. This is the same limit described with a sign reversal; it does not change the centering and scale for the conventional normalized maximum.

A General Verification Pattern

The three calculations share one structure: choose a threshold \(b_n+a_nx\), find the probability \(u_n(x)=1-F(b_n+a_nx)\) that one observation exceeds it, and then use

$$ P\left(\frac{M_n-b_n}{a_n}\leq x\right) =\bigl(1-u_n(x)\bigr)^n. $$

When \(n u_n(x)\to\tau(x)\), the proved limit gives a maximum CDF of \(\exp(-\tau(x))\). For the exponential model, \(\tau(x)=e^{-x}\); for the Pareto model, \(\tau(x)=x^{-\alpha}\) on \(x>0\); and for the uniform model below the endpoint, \(\tau(x)=-x\) on \(x<0\). The resulting CDFs are evaluated on their full support, with zero or one assigned outside the region where the threshold is inside the distribution’s support.

Common Pitfalls

  • Using the scale in the wrong direction. Writing \(a_n(M_n-b_n)\) instead of \((M_n-b_n)/a_n\) changes the normalization. State the normalized variable explicitly before calculating its CDF.
  • Centering a Pareto maximum by a logarithm. Exponential tails produce logarithmic growth, but the Pareto tail requires the power scale \(n^{1/\alpha}\).
  • Forgetting the finite endpoint. For a uniform sample, thresholds above 1 have CDF 1. Extending the expression \(1+x/n\) as though it were a valid uniform CDF for \(x>0\) is incorrect.
  • Dropping the \(n\)th power. A single observation’s exceedance probability tends to zero, but the maximum depends on all \(n\) observations. The factor \(n\) in \(n u_n(x)\) determines the nondegenerate limit.
Key takeaway: Derive a maximum normalization by making the one-observation exceedance probability at the normalized threshold of order \(1/n\). Then verify the proposed limit using the exact identity \(P(M_n\leq t)=F(t)^n\), including the support boundaries.

Check Your Understanding

Use the maximum CDF identity and the tail behavior of one observation to answer these questions.

  1. For exponential observations with rate \(\lambda\), what centering and scale give the standard Gumbel limit?
  2. For a Pareto distribution with shape \(\alpha\), why does the maximum scale involve \(n^{1/\alpha}\) rather than \(\log n\)?
  3. For \(X_i\sim\operatorname{Uniform}(0,1)\), write the limiting CDF of \(n(M_n-1)\) on both sides of its endpoint.
  4. If the exceedance probability at a normalized threshold is \(u_n\) with \(nu_n\to\tau\), what limit does the maximum CDF have there?
  5. Explain why the endpoint-gap variable \(n(1-M_n)\) has a different-looking limit from \(n(M_n-1)\).
Clinical Biostats