Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Hill Estimator

Overview

How heavy is the upper tail, and how sensitive is the answer to where we start fitting it? Hill estimates the exponent of a positive Pareto-type tail from the largest observations. Choosing how many observations to retain is part of the estimate.

Read the Pareto and regular variation pages first. The formula uses logarithms and sorted observations; the extreme-value index ξ\xi is the reciprocal of the Pareto exponent α\alpha. For data analysis, inspect several thresholds and check dependence and measurement limits before interpreting a stable-looking curve. A dedicated threshold-selection page is planned.

Estimator

Let x1,,xnx_1,\dots,x_n be observed right-tail data with xi>0x_i>0 for every ii, and write the ascending realized order statistics as

x1:nxn:n.x_{1:n}\le \cdots \le x_{n:n}.

For 1k<n1\le k<n, the Hill estimate of the extreme-value index is

ξ^k,n=1kj=1klog(xnj+1:nxnk:n).\widehat\xi_{k,n} = \frac1k\sum_{j=1}^{k} \log\left(\frac{x_{n-j+1:n}}{x_{n-k:n}}\right).

For a Pareto-type right tail with exponent α\alpha, the corresponding tail exponent estimate is

α^k,n=1ξ^k,n.\widehat\alpha_{k,n}=\frac{1}{\widehat\xi_{k,n}}.

The tuning parameter kk is the number of upper order statistics used. A Hill stability plot graphs ξ^k,n\widehat\xi_{k,n} or α^k,n\widehat\alpha_{k,n} over a range of kk values. On this page, we name the plotted coordinate explicitly: ξ\xi is the extreme-value index, while α=1/ξ\alpha=1/\xi is the Pareto tail exponent.

Here xi:nx_{i:n} is the iith ascending observation, and kk selects the top kk values above the order-statistic threshold. If all selected values equal that threshold, ξ^=0\widehat\xi=0 and there is no finite reciprocal estimate. Ties can also make the strict exceedance count smaller than kk.

When to use it

Hill is a right-tail estimator for strictly positive observations in a Pareto-type regime, meaning the survival tail is expected to be regularly varying. It is appropriate only after deciding what data transformation makes the tail positive and one-sided.

Small kk means the threshold is very high, so the estimate uses only the most extreme observations and has high variance. Large kk lowers the threshold, which adds data but can mix non-tail observations into the calculation. A stable region is a practical compromise between those two failures.

Report a range of kk values, the corresponding thresholds, and the sensitivity of the estimated exponent. A plateau can guide further checks, but it does not establish regular variation or independence.

Exact Pareto calibration

For a Pareto Type I variable,

P(X>x)=(xmx)α,xxm.\mathbb P(X>x)=\left(\frac{x_m}{x}\right)^\alpha,\qquad x\ge x_m.

Define Y=log(X/xm)Y=\log(X/x_m). For y0y\ge0,

P(Y>y)=P(X>xmey)=eαy,\mathbb P(Y>y) = \mathbb P(X>x_m e^y) =e^{-\alpha y},

so YY is exponential with mean 1/α1/\alpha. Equivalently,

E[log(Xu)X>u]=1α\mathbb E\left[\log\left(\frac{X}{u}\right)\mid X>u\right] =\frac1\alpha

for any Pareto threshold uxmu\ge x_m. The Hill estimator replaces this conditional expectation by the empirical average of log-excesses above the random threshold Xnk:nX_{n-k:n}, or by xnk:nx_{n-k:n} after the sample is realized.

For exact Pareto samples this explains why ξ^k,n\widehat\xi_{k,n} targets 1/α1/\alpha. For a positive iid sample whose right tail is regularly varying with extreme-value index ξ>0\xi>0, the standard consistency result uses an intermediate sequence k=knk=k_n with

kn,knn0.k_n\to\infty,\qquad \frac{k_n}{n}\to0.

Under these assumptions, ξ^kn,n\widehat\xi_{k_n,n} converges in probability to ξ\xi. The proof of that full theorem belongs to extreme-value asymptotics; we use the exact Pareto calibration and cite the general result in Resnick 2007. Hill’s original estimator is described in his 1975 paper, pp. 1163–1174.

Stability diagnostics

A stability diagnostic repeats the same calculation over kk and reports both the threshold and the chosen tail coordinate. The arithmetic is visible even in a fixed four-observation example.

Take the fixed ordered observations (1,2,4,8)(1,2,4,8). The Hill calculation gives:

kkThreshold xnk:nx_{n-k:n}ξ^k,n\widehat\xi_{k,n}α^k,n\widehat\alpha_{k,n}
14log2\log 21/log21/\log 2
223log2/23\log 2/22/(3log2)2/(3\log 2)
312log22\log 21/(2log2)1/(2\log 2)

This small deterministic example shows how the selected threshold changes the estimate. It is an illustration of the arithmetic, with no claim that the four values are a random Pareto sample.

In a distribution with a bounded body and a Pareto tail, increasing kk far enough pulls body observations into the calculation. For empirical data, a plateau is evidence to inspect, not a certificate.

Failure modes

References

Source and adaptation

Adapted from incerto-wiki, content/concepts/methods/hill-estimator.md, revision 9717c9c (2026-09-13 import). Copyright (c) 2023 xshi19. Licensed under MIT. Links and notation were adapted for this site; executable figures and simulations were replaced with static calculations. No upstream execution or formal-proof verification is claimed for this adaptation.

MIT permission notice

MIT License

Copyright (c) 2023 xshi19

Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the “Software”), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions:

The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software.

THE SOFTWARE IS PROVIDED “AS IS”, WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

References
  1. Resnick, S. I. (2007). Heavy-Tail Phenomena: Probabilistic and Statistical Modeling. In Springer Series in Operations Research and Financial Engineering. Springer New York. 10.1007/978-0-387-45024-7
  2. Hill, B. M. (1975). A Simple General Approach to Inference About the Tail of a Distribution. The Annals of Statistics, 3(5). 10.1214/aos/1176343247
  3. Embrechts, P., Klüppelberg, C., & Mikosch, T. (1997). Modelling Extremal Events. Springer Berlin Heidelberg. 10.1007/978-3-642-33483-2