Misplaced Pages

Gumbel distribution: Difference between revisions

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
Browse history interactively← Previous editNext edit →Content deleted Content addedVisualWikitext
Revision as of 18:02, 22 October 2013 editGareth Jones (talk | contribs)Autopatrolled, Extended confirmed users8,840 edits Application: complete ref← Previous edit Revision as of 04:01, 23 October 2013 edit undoAManWithNoPlan (talk | contribs)Extended confirmed users95,771 edits | Assisted by Citation bot r454Next edit →
Line 60: Line 60:
] ]


When the cdf of ''Y'' is the converse of the Gumbel standard cumulative distribution, <math>P(Y \leq y) = 1 - F(y)</math>, then ''Y'' has a density function that is a ]:<ref>Willemse, W. J. and Kaas, R., "Rational reconstruction of frailty-based mortality models by a generalisation of Gompertz’ law of mortality", ''Insurance: Mathematics and Economics'', 40 (3) (2007), 468&ndash;484.</ref> however, ''Y'' does not have a ] since the Gompertz distribution is restricted to positive values while ''Y'' can take both positive and negative values. When the cdf of ''Y'' is the converse of the Gumbel standard cumulative distribution, <math>P(Y \leq y) = 1 - F(y)</math>, then ''Y'' has a density function that is a ]:<ref>{{Cite journal|doi=10.1016/j.insmatheco.2006.07.003|title=Rational reconstruction of frailty-based mortality models by a generalisation of Gompertz' law of mortality|year=2007|last1=Willemse|first1=W.J.|last2=Kaas|first2=R.|journal=Insurance: Mathematics and Economics|volume=40|issue=3|pages=468}}</ref> however, ''Y'' does not have a ] since the Gompertz distribution is restricted to positive values while ''Y'' can take both positive and negative values.


Theory related to the ] provides a multivariate version of the Gumbel distribution. Theory related to the ] provides a multivariate version of the Gumbel distribution.
Line 72: Line 72:
] with ] of a cumulative Gumbel distribution to maximum one-day October rainfalls.]] ] with ] of a cumulative Gumbel distribution to maximum one-day October rainfalls.]]


Gumbel has shown that the maximum value (or last ]) in a sample of a ] following an ] approaches the Gumbel distribution closer with increasing sample size.<ref>Gumbel, E.J. 1954. "Statistical theory of extreme values and some practical applications". ''Applied Mathematics Series'', 33. U.S. Department of Commerce, National Bureau of Standards.</ref> Gumbel has shown that the maximum value (or last ]) in a sample of a ] following an ] approaches the Gumbel distribution closer with increasing sample size.<ref>{{cite book|last=Gumbel|first= E. J. |year=1954|asin=B0007DSHG4| title=Statistical theory of extreme values and some practical applications|series=Applied Mathematics Series|volume= 33|edition=1st|url=http://books.google.com/books/about/Statistical_theory_of_extreme_values_and.html?id=SNpJAAAAMAAJ|lccn=QA 1 .U58 no.33|publisher= U.S. Department of Commerce, National Bureau of Standards}}</ref>


In ], therefore, the Gumbel distribution is used to analyze such variables as monthly and annual maximum values of daily rainfall and river discharge volumes,<ref>{{cite book|last=Ritzema (ed.)|first=H.P.|title=Frequency and Regression Analysis|year=1994|publisher=Chapter 6 in: Drainage Principles and Applications, Publication 16, International Institute for Land Reclamation and Improvement (ILRI), Wageningen, The Netherlands|pages=175–224|url=http://www.waterlog.info/pdf/freqtxt.pdf|isbn=90-70754-33-9}}</ref> and also to describe droughts.<ref>{{cite doi|10.1016/j.jhydrol.2010.04.035}}</ref> In ], therefore, the Gumbel distribution is used to analyze such variables as monthly and annual maximum values of daily rainfall and river discharge volumes,<ref>{{cite book|editor-last=Ritzema |editor-first=first=H. P.|first1=R. J.|last1=Oosterbaan|chapter=Chapter 6 Frequency and Regression Analysis|year=1994|title= Drainage Principles and Applications, Publication 16|publisher= International Institute for Land Reclamation and Improvement (ILRI)|location= Wageningen, The Netherlands|pages=175–224|url=http://www.waterlog.info/pdf/freqtxt.pdf|isbn=90-70754-33-9}}</ref> and also to describe droughts.<ref>{{cite journal|doi=10.1016/j.jhydrol.2010.04.035|title=An extreme value analysis of UK drought and projections of change in the future|year=2010|last1=Burke|first1=Eleanor J.|last2=Perry|first2=Richard H.J.|last3=Brown|first3=Simon J.|journal=Journal of Hydrology|volume=388|pages=131}}</ref>
Gumbel has also shown that the ] <big>''r'' / (''n''+1)</big> for the probability of an event - where <big>''r''</big> is the rank number of the observed value in the data series and <big>''n''</big> is the total number of observations - is an ] of the ] around the ] of the distribution. Therefore, this estimator is often used as a ]. Gumbel has also shown that the ] <big>''r'' / (''n''+1)</big> for the probability of an event - where <big>''r''</big> is the rank number of the observed value in the data series and <big>''n''</big> is the total number of observations - is an ] of the ] around the ] of the distribution. Therefore, this estimator is often used as a ].
Line 80: Line 80:
The blue picture illustrates an example of fitting the Gumbel distribution to ranked maximum one-day October rainfalls showing also the 90% ] based on the ]. The rainfall data are represented by the ] ''r'' / (''n''+1) as part of the ]. The blue picture illustrates an example of fitting the Gumbel distribution to ranked maximum one-day October rainfalls showing also the 90% ] based on the ]. The rainfall data are represented by the ] ''r'' / (''n''+1) as part of the ].


In ], the Gumbel distribution approximates the number of terms in a ]<ref>{{cite doi|10.1215/S0012-7094-41-00826-8}}</ref> as well as the trend-adjusted sizes of record ] and record gaps between ].<ref>Kourbatov, A. "Maximal gaps between prime k-tuples: a statistical approach". ''Journal of Integer Sequences'', 16 (2013). Article 13.5.2.</ref> In ], the Gumbel distribution approximates the number of terms in a ]<ref>{{cite journal|doi=10.1215/S0012-7094-41-00826-8|title=Of a positive integer|year=1941|last1=Erdös|first1=Paul|last2=Lehner|first2=Joseph|journal=Duke Mathematical Journal|volume=8|issue=2|pages=335}}</ref> as well as the trend-adjusted sizes of record ] and record gaps between ].<ref>{{cite journal|arxiv=1301.2242|last=Kourbatov|first= A.|title=Maximal gaps between prime k-tuples: a statistical approach|journal=Journal of Integer Sequences|volume= 16 |year=2013}} Article 13.5.2.</ref>


==See also== ==See also==

Revision as of 04:01, 23 October 2013

Gumbel
Probability density functionProbability distribution function
Cumulative distribution functionCumulative distribution function
Parameters μ {\displaystyle \mu \!} location (real)
β > 0 {\displaystyle \beta >0\!} scale (real)
Support x ( ; + ) {\displaystyle x\in (-\infty ;+\infty )\!}
PDF 1 β e ( z + e z ) {\displaystyle {\frac {1}{\beta }}e^{-(z+e^{-z})}\!}
where z = x μ β {\displaystyle z={\frac {x-\mu }{\beta }}\!}
CDF e e ( x μ ) / β {\displaystyle e^{-e^{-(x-\mu )/\beta }}\!}
Mean μ + β γ {\displaystyle \mu +\beta \,\gamma \!}
Median μ β ln ( ln ( 2 ) ) {\displaystyle \mu -\beta \,\ln(\ln(2))\!}
Mode μ {\displaystyle \mu \!}
Variance π 2 6 β 2 {\displaystyle {\frac {\pi ^{2}}{6}}\,\beta ^{2}\!}
Skewness 12 6 ζ ( 3 ) π 3 1.14 {\displaystyle {\frac {12{\sqrt {6}}\,\zeta (3)}{\pi ^{3}}}\approx 1.14\!}
Excess kurtosis 12 5 {\displaystyle {\frac {12}{5}}}
Entropy ln ( β ) + γ + 1 {\displaystyle \ln(\beta )+\gamma +1\!}
MGF Γ ( 1 β t ) e μ t {\displaystyle \Gamma (1-\beta \,t)\,e^{\mu \,t}\!}
CF Γ ( 1 i β t ) e i μ t {\displaystyle \Gamma (1-i\,\beta \,t)\,e^{i\,\mu \,t}\!}

In probability theory and statistics, the Gumbel distribution is used to model the distribution of the maximum (or the minimum) of a number of samples of various distributions. Such a distribution might be used to represent the distribution of the maximum level of a river in a particular year if there was a list of maximum values for the past ten years. It is useful in predicting the chance that an extreme earthquake, flood or other natural disaster will occur.

The potential applicability of the Gumbel distribution to represent the distribution of maxima relates to extreme value theory which indicates that it is likely to be useful if the distribution of the underlying sample data is of the normal or exponential type.

The Gumbel distribution is a particular case of the generalized extreme value distribution (also known as the Fisher-Tippett distribution). It is also known as the log-Weibull distribution and the double exponential distribution (a term that is alternatively sometimes used to refer to the Laplace distribution). It is often incorrectly labelled as Gompertz distribution.

In the latent variable formulation of the multinomial logit model — common in discrete choice theory — the errors of the latent variables follow a Gumbel distribution. This is useful because the difference of two Gumbel-distributed random variables has a logistic distribution.

The Gumbel distribution is named after Emil Julius Gumbel (1891–1966).

Properties

This section includes a list of references, related reading, or external links, but its sources remain unclear because it lacks inline citations. Please help improve this section by introducing more precise citations. (February 2012) (Learn how and when to remove this message)

The cumulative distribution function of the Gumbel distribution is

F ( x ; μ , β ) = e e ( x μ ) / β . {\displaystyle F(x;\mu ,\beta )=e^{-e^{-(x-\mu )/\beta }}.\,}

The mode is μ, while the median is μ β ln ( ln 2 ) , {\displaystyle \mu -\beta \ln \left(\ln 2\right),} and the mean is given by

E ( X ) = μ + γ β , {\displaystyle \operatorname {E} (X)=\mu +\gamma \beta ,}

where γ {\displaystyle \gamma } = Euler–Mascheroni constant 0.5772. {\displaystyle \approx 0.5772.} The standard deviation is β π / 6 . {\displaystyle \beta \pi /{\sqrt {6}}.}

Standard Gumbel distribution

The standard Gumbel distribution is the case where μ = 0 {\displaystyle \mu =0} and β = 1 {\displaystyle \beta =1} with cumulative distribution function

F ( x ) = e e ( x ) {\displaystyle F(x)=e^{-e^{(-x)}}\,}

and probability density function

f ( x ) = e ( x + e x ) . {\displaystyle f(x)=e^{-(x+e^{-x})}.}

In this case the mode is 0, the median is ln ( ln ( 2 ) ) 0.3665 {\displaystyle -\ln(\ln(2))\approx 0.3665} , the mean is γ {\displaystyle \gamma } , and the standard deviation is π / 6 1.2825. {\displaystyle \pi /{\sqrt {6}}\approx 1.2825.}

Quantile function and Generating Gumbel variates

Since the quantile function(inverse cumulative distribution function), Q ( p ) {\displaystyle Q(p)} , of a Gumbel distribution is given by

Q ( p ) = μ β ln ( ln ( p ) ) , {\displaystyle Q(p)=\mu -\beta \ln(-\ln(p)),}

the variate Q ( U ) {\displaystyle Q(U)} has a Gumbel distribution with parameters μ {\displaystyle \mu } and β {\displaystyle \beta } when the random variate U {\displaystyle U} is drawn from the uniform distribution on the interval ( 0 , 1 ) {\displaystyle (0,1)} .

Related distributions

A piece of graph paper that incorporates the Gumbel distribution.

When the cdf of Y is the converse of the Gumbel standard cumulative distribution, P ( Y y ) = 1 F ( y ) {\displaystyle P(Y\leq y)=1-F(y)} , then Y has a density function that is a Gompertz function: however, Y does not have a Gompertz distribution since the Gompertz distribution is restricted to positive values while Y can take both positive and negative values.

Theory related to the generalized multivariate log-gamma distribution provides a multivariate version of the Gumbel distribution.

Graphic paper

In pre-software times graphic paper was used to picture the Gumbel distribution (see illustration). The paper is based on linearization of the cumulative distribution function F {\displaystyle F}  :

l n [ l n ( F ) ] = ( x μ ) / β {\displaystyle -ln=(x-\mu )/\beta }

In the paper the horizontal axis is constructed at a double log scale. The vertical axis is linear. By plotting F {\displaystyle F} on the horizontal axis of the paper and the x {\displaystyle x} -variable on the vertical axis, the distribution is represented by a straight line with a slope 1 / β {\displaystyle /\beta } . When distribution fitting software like CumFreq became available, the task of plotting the distribution was made easier, as is demonstrated in the section below.

Application

Distribution fitting with confidence belt of a cumulative Gumbel distribution to maximum one-day October rainfalls.

Gumbel has shown that the maximum value (or last order statistic) in a sample of a random variable following an exponential distribution approaches the Gumbel distribution closer with increasing sample size.

In hydrology, therefore, the Gumbel distribution is used to analyze such variables as monthly and annual maximum values of daily rainfall and river discharge volumes, and also to describe droughts.

Gumbel has also shown that the estimator r / (n+1) for the probability of an event - where r is the rank number of the observed value in the data series and n is the total number of observations - is an unbiased estimator of the cumulative probability around the mode of the distribution. Therefore, this estimator is often used as a plotting position.

The blue picture illustrates an example of fitting the Gumbel distribution to ranked maximum one-day October rainfalls showing also the 90% confidence belt based on the binomial distribution. The rainfall data are represented by the plotting position r / (n+1) as part of the cumulative frequency analysis.

In number theory, the Gumbel distribution approximates the number of terms in a partition of an integer as well as the trend-adjusted sizes of record prime gaps and record gaps between prime constellations.

See also

External links

References

  1. Willemse, W.J.; Kaas, R. (2007). "Rational reconstruction of frailty-based mortality models by a generalisation of Gompertz' law of mortality". Insurance: Mathematics and Economics. 40 (3): 468. doi:10.1016/j.insmatheco.2006.07.003.
  2. Gumbel, E. J. (1954). Statistical theory of extreme values and some practical applications. Applied Mathematics Series. Vol. 33 (1st ed.). U.S. Department of Commerce, National Bureau of Standards. ASIN B0007DSHG4. LCCN 1 .U58 no.33 QA 1 .U58 no.33. {{cite book}}: Check |lccn= value (help)
  3. Oosterbaan, R. J. (1994). "Chapter 6 Frequency and Regression Analysis". In Ritzema, first=H. P. (ed.). Drainage Principles and Applications, Publication 16 (PDF). Wageningen, The Netherlands: International Institute for Land Reclamation and Improvement (ILRI). pp. 175–224. ISBN 90-70754-33-9. {{cite book}}: Missing pipe in: |editor-first= (help)
  4. Burke, Eleanor J.; Perry, Richard H.J.; Brown, Simon J. (2010). "An extreme value analysis of UK drought and projections of change in the future". Journal of Hydrology. 388: 131. doi:10.1016/j.jhydrol.2010.04.035.
  5. Erdös, Paul; Lehner, Joseph (1941). "Of a positive integer". Duke Mathematical Journal. 8 (2): 335. doi:10.1215/S0012-7094-41-00826-8.
  6. Kourbatov, A. (2013). "Maximal gaps between prime k-tuples: a statistical approach". Journal of Integer Sequences. 16. arXiv:1301.2242. Article 13.5.2.
Probability distributions (list)
Discrete
univariate
with finite
support
with infinite
support
Continuous
univariate
supported on a
bounded interval
supported on a
semi-infinite
interval
supported
on the whole
real line
with support
whose type varies
Mixed
univariate
continuous-
discrete
Multivariate
(joint)
Directional
Univariate (circular) directional
Circular uniform
Univariate von Mises
Wrapped normal
Wrapped Cauchy
Wrapped exponential
Wrapped asymmetric Laplace
Wrapped Lévy
Bivariate (spherical)
Kent
Bivariate (toroidal)
Bivariate von Mises
Multivariate
von Mises–Fisher
Bingham
Degenerate
and singular
Degenerate
Dirac delta function
Singular
Cantor
Families
Categories: