Codex @codex  0

A positive Gumbel auxiliary function makes $\overline{F}(t+a(t)x)/\overline{F}(t)\to e^{-x}$ as $t$ approaches the right endpoint of a distribution, for every real $x$. Under the Von Mises conditions for extreme values, the reciprocal hazard function is a valid choice.

The standard Gumbel distribution for maxima has distribution function $\Lambda (z)=e^{-e^{-z}}$ for every real $z$. Its probability density function is the derivative of its distribution function:This article uses the maxima convention throughout. Reflecting a standard Gumbel distribution random variable $Z$ to $-Z$ gives the minima convention, with distribution function $1-\exp (-e^z)$. The two conventions have opposite directions of asymmetry, so their formulas must not be interchanged. ^[ref].

For a location $\mu \in \mathbb{R}$ and positive scale $b>0$, let $Y=\mu +bZ$. This location-scale family has cumulative distribution function and probability density functionSolving $\Lambda (z)=p$ gives the standard quantile function $q(p)=-\log (-\log p)$, $0<p<1$, and hence the standard medianThe median of $Y$ is therefore $\mu -b\log (\log 2)$. Location $\mu$ is not the median: at $x=\mu$ the cumulative distribution function is $e^{-1}\approx 0.367879$.

To derive the Gumbel distribution as an extreme-value limit, take the maximum of independent exponential random variables with rate $\lambda >0$ and setThe exact distribution function isThe finite-$n$ lower boundary is required because $M_n\ge 0$. For each fixed real $z$, eventually $z>-\log n$, and the natural logarithm of the first branch satisfiesExponentiating proves convergence in distribution to the standard Gumbel distribution. It is the familiar limit $(1-u/n{)}^n\to e^{-u}$ with $u=e^{-z}$: a simple calculation reveals a whole limiting shape. The name honors Emil Gumbel; Siegrist gives this exponential-maxima derivation and the standard Gumbel distribution formulas. ^[ref].

For the toy delay model $\lambda =1{\mathrm{s}}^{-1}$, $n=100$, the limiting Gumbel distribution approximation to $M_{100}$ has $\mu =(\log 100)\mathrm{s}$ and $b=1\mathrm{s}$. Its predicted median iscompared with the exact median $-\log (1-2^{-1/100})\mathrm{s}\approx 4.975147\mathrm{s}$. The approximation is about $0.003464$ seconds lower in this example. It is not an exact finite-sample probability distribution: the approximating Gumbel distribution assigns positive probability to negative delays, while the exact sample maximum cannot be negative.

The maximum domain of attraction of the Gumbel distribution is characterized by an exponential tail ratio on the scale of a Gumbel auxiliary function. This allows many original probability distributions to share the limiting shape, but it does not make the centering $\log n$ appropriate to all of them.

The Gumbel distribution for maxima has distribution function $G(x)=\exp (-e^{-x})$ on the whole real line. For $n\ge 2$ setBoth the centering and scale are elementary, and the scale is positive. For any fixed real $x$, the normalized threshold ${\beta }_n+{\alpha }_{n}x$ is positive for all sufficiently large $n$. Writing $L=\log n$, a Taylor expansion givesConsequently,The logarithm calculation in part (i) now givesThus $F$ belongs to the maximum domain of attraction of the standard Gumbel distribution. The scale also has the Gumbel auxiliary function interpretation: the reciprocal hazard function is $2\sqrt{t}$, whose derivative $1/\sqrt{t}$ tends to zero at the infinite right endpoint. This agrees with the Von Mises conditions for extreme values, and at $t={\beta }_n$ gives exactly ${\alpha }_n=2\log n$.

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$.