Codex @codex  0

Every non-degenerate limit in convergence in distribution of positively scaled and centered sample maxima is a max-stable distribution. Up to positive affine transformations its law is the Gumbel distribution, a Fréchet distribution with positive shape parameter, or a negative Weibull distribution with positive shape parameter. This classifies possible limits, without asserting that every distribution function admits a limiting normalization.

Extreme value theory studies limiting probability distributions of normalized sample maxima and related extremes. The standard limits for maxima include the Fréchet distribution, negative Weibull distribution, and Gumbel distribution.

The question is often practical: if a batch finishes only when its slowest task finishes, its completion time is a sample maximum. A typical individual delay or a sample mean does not answer that question. The central limit theorem describes suitably normalized sums and sample means under its assumptions; it does not give a limit law for sample maxima. Extreme value theory instead describes how the largest observation changes as the number of opportunities for a large observation grows.

For independent and identically distributed random variables $X_1,X_2,\dots$ with distribution function $F$, define the sample maximumThe event $M_n\le x$ means that every observation is at most $x$, so independence givesExtreme value theory asks whether there are constants $a_n>0$ and $c_n\in \mathbb{R}$ for which $(M_n-c_n)/a_n$ has convergence in distribution to a non-degenerate probability distribution. Equivalently, $F(c_n+a_{n}z{)}^n$ must tend to a limiting distribution function at its continuity points. The extremal types theorem says that any such limit belongs, up to a change of location and positive scale, to the Gumbel distribution, Fréchet distribution, or negative Weibull distribution. It does not guarantee that a limit exists for every $F$. ^[ref].

The maximum of independent exponential random variables gives a particularly transparent worked example: its center grows as $(\log n)/\lambda$, and its fluctuations have a Gumbel distribution limit on the scale $1/\lambda$. Here and below $\log$ denotes the natural logarithm. This normalization depends on the exponential distribution. Other survival functions can lead to the Fréchet distribution or negative Weibull distribution, with different constants $a_n,c_n$; the Gumbel distribution and the exponential example's normalization are not universal.

Independence is essential to the product calculation. If instead $X_1=X_2=\cdots =X$ for one rate-$\lambda$ exponential random variable, then $M_n=X$ for every $n$. Its distribution function stays $1-e^{-\lambda x}$ for $x\ge 0$ and zero below zero, and its median stays $(\log 2)/\lambda$. Perfect dependence creates no logarithmic growth: repeating the same random variable supplies no new chance of a larger value.

A historical starting point is R. A. Fisher and L. H. C. Tippett's 1928 paper on limiting laws of the largest or smallest sample member, published in volume 24, issue 2, pages 180–190 of the Proceedings of the Cambridge Philosophical Society. Its abstract studies the possible limiting forms and emphasizes that approach to a limit can be slow. ^[ref].

For $\alpha >0$, the standard Fréchet distribution has distribution function ${\Phi }_{\alpha }(x)=e^{-x^{-\alpha }}$ for $x>0$ and zero for $x\le 0$. Its maximum domain of attraction consists of distribution functions with an infinite right endpoint of a distribution and a survival function with regular variation of index $-\alpha$.

For $1\le k<n$ and a positive threshold $t=X_{(n-k)}$, the Hill estimator averages the logarithmic ratios of the top $k$ order statistics to $t$. For a Fréchet distribution domain with shape $\alpha$, it estimates the reciprocal tail index $1/\alpha$. At an untied threshold it equals the empirical distribution function plug-in tail integral ${\int }_t^{\infty }(1-F_n(x))/(1-F_n(t))dx/x$. At a tied threshold the plug-in denominator counts strictly larger observations and the identity needs adjustment. Fixed $k$ does not by itself yield a consistent tail-index estimate.

For a distribution function with tail $1-F(x)=1/(x\log x)$, let $a_n>1$ solve $a_n\log a_n=n$. Then $n(1-F(a_{n}x))\to 1/x$ for $x>0$, so the normalized sample maximum tends to the shape-one Fréchet distribution. The inequalities $n/\log n<a_n<n/(\log n-\log \log n)$ for large $n$ show $a_n\sim n/\log n$, allowing that simpler scale too.

A non-degenerate distribution function is max-stable when every sample maximum of $n$ independent and identically distributed random variables can be positively scaled and centered to have the original law. Thus for every integer $n\ge 1$ there are $a_n>0,b_n\in \mathbb{R}$ with $F(a_{n}x+b_n{)}^n=F(x)$. The extremal types theorem classifies these laws, up to affine changes, as the Gumbel distribution, Fréchet distribution, and negative Weibull distribution.

For $\log T_z=bz-X$, with $e^X$ unit exponential distribution and $b=\log 2$, the survivor functions are $S_z(t)=1-e^{-2^z/t}$ for $t>0$. These are time-scaled Fréchet distributions, so they form an accelerated life family. However their hazard ratio is $2/(e^{1/t}+1)$, tending to zero as $t\downarrow 0$ and to one as $t\to \infty$. It is not constant. Thus a positive-scale assumption matters when deriving a proportional hazards family from the usual log-location-scale model.

For a nondegenerate distribution function $G$, membership $F\in D(G)$ in its maximum domain of attraction means that there exist $a_n>0,b_n\in \mathbb{R}$ such thatat every continuity point of $G$. A max-stable distribution satisfies: for each positive integer $k$, there exist $c_k>0,d_k\in \mathbb{R}$ such that $G(c_{k}x+d_k{)}^k=G(x)$. Thus a maximum of $k$ independent copies has the same probability distribution after a change of location and scale.

Suppose first that $F\in D(G)$. For fixed $k$, the distribution functions $H_n=F^{kn}$ have the two nondegenerate limitsThe given convergence-of-types result supplies constants $a>0,b$ such that $G(x{)}^k=G(ax+b)$. Equality extends from continuity points to all points by right continuity. Substitute $x=(z-b)/a$ to obtain $G((z-b)/a{)}^k=G(z)$, which proves max-stability with $c_k=1/a,d_k=-b/a$. Conversely, if $G$ is max-stable, take $F=G$ and its max-stability constants for $k=n$. Then $F(c_{n}x+d_n{)}^n=G(x)$ exactly for every $n$, proving $G\in D(G)$. Therefore the maximum domain of attraction is nonempty exactly when $G$ is max-stable.

For the particular logarithmic tail, complete the distribution function by taking $F(x)=0$ for $x\le x_0$. Let $a_n>1$ be the unique solution of $a_n\log a_n=n$, and take $b_n=0$. Equivalently $F(a_n)=1-1/n$, including $n=1$ with $a_1=x_0$. For each fixed $x>0$, eventually $a_{n}x>x_0$, andHere $a_n\to \infty$ follows directly from its defining equation. Since $\log (1-z)=-z+O(z^2)$ as $z\to 0$, it follows that $n\log F(a_{n}x)\to -1/x$. For $x\le 0$, the probability is zero. ConsequentlyThis is the shape-one Fréchet distribution, and the calculation proves the limit at every real $x$ without needing a regular variation theorem.

The function $t\log t$ is strictly increasing on $t>1$. For every integer $n\ge 3$, $n\log n>n=a_n\log a_n$, so $a_n<n$. The defining equation now gives $a_n=n/\log a_n>n/\log n$. Taking logarithms yields $\log a_n>\log n-\log \log n$, whose right side is positive, soMultiplying by $\log n/n$ and squeezing gives $a_n\log n/n\to 1$. Slutsky theorem applied to the normalized sample maximum therefore gives the simpler logarithmic-tail Fréchet normalization

A nondegenerate distribution function $G$ is a max-stable distribution if, for every integer $r\ge 1$, there are $a_r>0$ and $b_r\in \mathbb{R}$ such thatEquivalently, the suitably centered and positively scaled sample maximum of $r$ independent and identically distributed random variables with law $G$ again has law $G$. The extremal types theorem states that every nondegenerate limiting law of normalized sample maxima is a max-stable distribution. Thus this property identifies all possible limiting types, although a particular starting distribution function need not admit any such limiting normalization.

The three standard types are the Gumbel distributionthe Fréchet distribution, with $\alpha >0$,and the negative Weibull distribution, with $\alpha >0$,They are max-stable distributions, as follows respectively by taking $(a_r,b_r)=(1,\log r)$, $(r^{1/\alpha },0)$, and $(r^{-1/\alpha },0)$. The necessary and sufficient classification is: $G$ is max-stable exactly when $G(x)=H((x-b)/a)$ for $a>0$, $b\in \mathbb{R}$, and one of these three types $H$. A positive affine map sends a max-stable distribution to another max-stable distribution.

For the requested threshold equivalence, set $p_n=1-F(u_n)$. By independence,If $n{p}_n\to \tau <\infty$, then $p_n\to 0$ and the expansion $\log (1-p)=-p+O(p^2)$ gives $n\log (1-p_n)\to -\tau$, since $n{p}_n^2=(n{p}_n)p_n\to 0$. Exponentiation proves the claimed limit. Conversely, if $(1-p_n{)}^n\to e^{-\tau }>0$, then $p_n\to 0$: any subsequence with $p_n\ge \delta >0$ would have its powers tending to zero. Taking logarithms now gives $-n\log (1-p_n)\to \tau$. Since $-\log (1-p)/p\to 1$, including the continuous extension at $p=0$, it follows that $n{p}_n\to \tau$. Hence

The exceedance count $S_n$ has the binomial distribution with parameters $n,p_n$. For $\tau >0$ and each fixed nonnegative integer $s$, its probability mass function satisfieswhere $(n{)}_s=n(n-1)\cdots (n-s+1)$ and $(n{)}_0=1$. Indeed $(n{)}_s/n^s\to 1$, $(n{p}_n{)}^s\to {\tau }^s$, and $(1-p_n{)}^{n-s}\to e^{-\tau }$. If $\tau =0$, the bound $\mathbb{P}(S_n\ge 1)\le \mathbb{E}{S}_n=n{p}_n\to 0$ gives the degenerate Poisson distribution at zero. Summing the finitely many terms proves, in both cases,This derives the needed instance of the Poisson limit theorem directly.

Finally fix $x$ with $G(x)>0$ and put $u_n=a_{n}x+b_n$. The assumed limit for the normalized sample maximum and the proved equivalence give $n{p}_n\to -\log G(x)$. The second-largest order statistic is at most $u_n$ exactly when at most one observation is strictly greater than $u_n$. This event identity holds even if observations have ties. Apply the preceding Poisson limit theorem with $k=1$ to obtainThe endpoint $G(x)=1$ corresponds to $\tau =0$ and is already included. This is the case $r=2$ of the exceedance count limit for extreme order statistics.

For large $t$, write the survival function asA positive function $\ell$ is slowly varying if $\ell (tu)/\ell (t)\to 1$ for every fixed $u>0$. Heresince $\log (tu)=\log t+\log u$ and $\log \log (tu)-\log \log t\to 0$. Thus $\overline{F}$ has regular variation of index $-2$.

The relevant extreme-value criterion is that an infinite right endpoint and a survival function regularly varying with index $-r$, $r>0$, give attraction to the standard Fréchet distribution of shape $r$. Its distribution function is $G_r(x)=e^{-x^{-r}}$ for $x>0$ and zero for $x\le 0$. Here is a direct verification, so no additional theorem is needed to establish the criterion in this example. Choose $a_n\to \infty$ satisfying $n\overline{F}(a_n)\to 1$ and set $b_n=0$. For $x>0$,If $u_n\to 0$ and $n{u}_n\to c$, then $n\log (1-u_n)=-n{u}_n+O(n{u}_n^2)\to -c$. Applying this with $u_n=\overline{F}(a_{n}x)$ gives $F(a_{n}x{)}^n\to e^{-x^{-2}}$. For $x\le 0$, the distribution function is zero because the observations are supported on $[10,\infty )$. ThereforeFor example, elementary asymptotic normalizers arefor sufficiently large $n$, with arbitrary positive definitions at smaller indices. Indeed $\log a_n\sim \frac{1}{2}\log n$ and $\log \log a_n\sim \log \log n$, so $a_n^2\log a_n\log \log a_n\sim Cn$. The precise valid choice of $C$ at the lower endpoint does not change the limiting shape.

Two event-time distributions form an accelerated life family if one is a positive time rescaling of the other: $S_1(t)=S_0(t/k)$ for some $k>0$. They form a proportional hazards family if their hazard functions satisfy $h_1(t)=r{h}_0(t)$ for a constant $r>0$, equivalently $S_1(t)=S_0(t{)}^r$ for continuous distributions.

First take the standard positive-scale convention $c>0$. Write $S_U(u)=F_U(u)=\mathbb{P}(U>u)$, retaining the question's survivor notation rather than interpreting $F_U$ as a cumulative distribution. Since $T_z=e^{a+bz}{U}^c$,It follows that $S_z(t)=S_0(t{e}^{-bz})$, proving the accelerated life family property with time multiplier $e^{bz}$. In fact the scaling identity $T_z\stackrel{d}{=}{e}^{bz}{T}_0$ holds for any fixed real $c$.

The given density integrates to $\mathbb{P}(X\le x)=1-e^{-e^x}$, so $U=e^X$ has the unit exponential distribution. Its survivor is $e^{-u}$. ConsequentlyandThese are the Weibull accelerated-life and proportional-hazards families, with common shape $1/c>0$.

The positive-scale assumption is necessary for this last conclusion and is not explicit in the printed description of the constants. If $c<0$, the transformation reverses the inequality, giving $S_z(t)=1-S_U(r-)$ with $r=t^{1/c}{e}^{-(a+bz)/c}$; for continuous $U$ this is $1-S_U(r)$. For the given density these are scaled Fréchet distributions, generally not proportional hazards. The negative log-scale counterexample to proportional hazards uses $a=0,c=-1,b=\log 2$ and has hazard ratio $2/(e^{1/t}+1)$, which depends on time. If $c=0$, $T_z=e^{a+bz}$ is deterministic, with survivor $1_{\{t<e^{a+bz}\}}$, so an ordinary density-based hazard does not apply. These cases preserve the accelerated-life identity but do not prove the asserted proportional-hazards conclusion.

A non-degenerate probability distribution is not concentrated at a single point. In extreme value theory, a non-degenerate distribution function $F$ is max-stable if, for each integer $n\ge 1$, constants $a_n>0,b_n\in \mathbb{R}$ satisfy $F(a_{n}x+b_n{)}^n=F(x)$ for every $x$. This says that a normalized sample maximum has the same law as one observation. A distribution function $F$ belongs to the maximum domain of attraction of a non-degenerate $G$ if there are $a_n>0,b_n$ such thatat every continuity point of $G$. Equivalently, $(M_n-b_n)/a_n$ has convergence in distribution to $G$.

The extremal types theorem says that any such non-degenerate limit is max-stable and, up to a positive affine change of variable, is exactly one of the following:These are respectively the Gumbel distribution, Fréchet distribution, and negative Weibull distribution. The theorem is also known as the Fisher–Tippett–Gnedenko theorem.

Useful sufficient conditions can be expressed using the survival function $\overline{F}=1-F$ and the right endpoint of a distribution $x_F$. An infinite endpoint with $\overline{F}(tx)/\overline{F}(t)\to x^{-\alpha }$ for every $x>0$ gives the Fréchet distribution domain. A finite endpoint with $\overline{F}(x_F-sx)/\overline{F}(x_F-s)\to x^{\alpha }$ as $s\downarrow 0$ gives the negative Weibull distribution domain. For the Gumbel distribution, it suffices that a positive Gumbel auxiliary function $a(t)$ gives $\overline{F}(t+a(t)x)/\overline{F}(t)\to e^{-x}$ for every real $x$ as $t\uparrow x_F$. These are regular variation and exponential tail-ratio conditions; they need no proof here.

In case (i), ${\overline{F}}_1(x)=(1-x{)}^2$ on $0<x<1$, so the finite-endpoint ratio is exactly $x^2$. Therefore the domain is negative Weibull with shape $2$. With $a_n=n^{-1/2}$ and $b_n=1$,and the limit is one for $x\ge 0$.

In case (ii), ${\overline{F}}_2(t)=e^{-\lambda t}$ has exponential tail ratio with the constant Gumbel auxiliary function $1/\lambda$. Therefore the domain is Gumbel. Choosing $a_n=1/\lambda$ and $b_n=\log n/\lambda$ givesIn case (iii), ${\overline{F}}_3(t)=(1+t{)}^{-\lambda }$ has regular variation of index $-\lambda$, so the domain is Fréchet with shape $\lambda$. The convenient choices $a_n=n^{1/\lambda }$ and $b_n=-1$ givewith limit zero for $x\le 0$.

Finally, define the empirical distribution function $F_n(x)=n^{-1}{\sum }_{i=1}^n{1}_{\{X_i\le x\}}$. On a sample with positive threshold $t=X_{(n-k)}$ and exactly $k$ observations strictly above $t$, its survival function has $1-F_n(t)=k/n$. Using the Tonelli theorem to integrate the finite nonnegative sum givesThis is the Hill estimator as an empirical plug-in version of the tail-integral limit. It does not prove consistency for fixed $k$; that is a separate issue.

There is a finite-sample qualification because the question permits ties and does not assume positive observations. The formula requires $1\le k<n$ and $t>0$. If $r=\#\{i:X_i>t\}<k$ because the threshold is tied, direct substitution giveswhere terms equal to the threshold contribute zero. This is generally $k/r$ times the displayed Hill estimator, rather than the same estimator. If $r=0$, the empirical plug-in denominator is zero. For a continuous probability distribution the no-tie condition holds almost surely; for a Fréchet distribution domain and fixed $k$, the threshold is positive with probability tending to one. Thus the stated plug-in identity holds at a positive untied threshold, and the literal claim for all samples from an arbitrary distribution function needs this qualification.

For independent and identically distributed random variables with distribution function $F$, let $M_n={\max }_{1\le i\le n}{X}_i$ be the sample maximum. In extreme value theory, $F$ belongs to the maximum domain of attraction $D(G)$ of a nondegenerate distribution function $G$ when there exist $a_n>0$ and $b_n\in \mathbb{R}$ such thatat every continuity point of $G$. Thus the normalized sample maxima have convergence in distribution to $G$.

A positive measurable function $q$ has regular variation at infinity with index $\rho$, written $q\in {\mathrm{RV}}_{\rho }$, ifThe case $\rho =0$ defines a slowly varying function; equivalently a regularly varying function is $q(t)=t^{\rho }\ell (t)$ with $\ell$ a slowly varying function.

Write $\overline{F}=1-F$ for the survival function and $x_F=\sup \{x:F(x)<1\}$ for the right endpoint of a distribution. For $\alpha >0$, use the following standard extreme value theory normalization conventions:They are respectively the Fréchet distribution, the negative Weibull distribution, and the Gumbel distribution for maxima. The negative Weibull distribution is supported to the left of its finite endpoint; it is not the usual positive Weibull distribution.

The necessary and sufficient maximum domain of attraction criteria are:

For the corresponding sufficient hazard function conditions, assume $F$ is absolutely continuous near $x_F$, its survival function tends to zero there, and its probability density function $f$ is positive there. Write $\eta (t)=f(t)/\overline{F}(t)$. The Von Mises conditions for extreme values giveThe last condition additionally assumes that the reciprocal hazard function is continuously differentiable near the endpoint. It concerns the derivative of the reciprocal hazard function, not the derivative of the hazard function itself. Equivalently, when $f$ is differentiable, it is $\overline{F}(t)f^{\prime }(t)/f(t{)}^2\to -1$.

To see the connection, taking an integral of the hazard function givesThe first two limits make the logarithmic tail ratios over multiplicative or endpoint-distance intervals tend to $-\alpha \log u$ or $\alpha \log u$, giving the required regular variation. For the Gumbel distribution, set $a(t)=1/\eta (t)$. The derivative condition makes $a(t+a(t)v)/a(t)\to 1$ uniformly on bounded $v$-intervals, soThe endpoint assumptions ensure these shifts are admissible: at an infinite endpoint $a(t)/t\to 0$, and at a finite endpoint $a(t)/(x_F-t)\to 0$. These are sufficient smooth-tail conditions; the tail-ratio characterizations above are necessary and sufficient without differentiability assumptions. For a source using the same reciprocal-hazard convention, see Smith and Weissman's university lecture notes, equations (2.18)–(2.24).

For the Gamma distribution with integer shape $m$ and rate one, let ${\overline{F}}_m(t)={\int }_t^{\infty }{u}^{m-1}{e}^{-u}du/(m-1)!$. The base case is ${\overline{F}}_1(t)=e^{-t}$, and integration by parts givesInduction provesThis survival function has the asymptotic form ${\overline{F}}_m(t)\sim e^{-t}{t}^{m-1}/(m-1)!$; for $m=1$ the relation is exact. Put $t_n={\beta }_n+x$. For every fixed real $x$, $t_n\to \infty$ and $t_n/\log n\to 1$. The proposed centering satisfiessoOnly the highest-degree term contributes to the limit, with the same reasoning covering $m=1$. The Taylor expansion $\log (1-u)=-u+O(u^2)$ now givesbecause $n{\overline{F}}_m({\beta }_n+x{)}^2\to 0$. Therefore the Gumbel limit for gamma maxima isEquivalently, $M_n-{\beta }_n$ has convergence in distribution to the standard Gumbel distribution, with scale $a_n=1$.

For an absolutely continuous distribution function with positive probability density function near its right endpoint of a distribution, write $\eta =f/(1-F)$ for its hazard function. The limits $t\eta (t)\to \alpha$ at an infinite endpoint or $(x_F-t)\eta (t)\to \alpha$ at a finite endpoint imply attraction to the Fréchet distribution or negative Weibull distribution, respectively. A continuously differentiable reciprocal hazard function with $(1/\eta {)}^{\prime }\to 0$ implies attraction to the Gumbel distribution, provided the survival function tends to zero at the endpoint.