Gumbel Distribution

DOWNLOAD Mathematica Notebook

There are essentially three types of Fisher-Tippett extreme value distributions. The most common is the type I distribution, which are sometimes referred to as Gumbel types or just Gumbel distributions. These are distributions of an extreme order statistic for a distribution of N elements X_i. In this work, the term "Gumbel distribution" is used to refer to the distribution corresponding to a minimum extreme value distribution (i.e., the distribution of the minimum X^(<1>)).

The Gumbel distribution with location parameter alpha and scale parameter beta is implemented in the Wolfram Language as GumbelDistribution[alpha, beta].

It has probability density function and distribution function

P(x)=1/betaexp[(x-alpha)/beta-exp((x-alpha)/beta)]
(1)
D(x)=1-exp[-exp((x-alpha)/beta)].
(2)

The mean, variance, skewness, and kurtosis are

mu=alpha-gammabeta
(3)
sigma^2=1/6pi^2beta^2
(4)
gamma_1=-(12sqrt(6)zeta(3))/(pi^3)
(5)
gamma_2=(12)/5,
(6)

where gamma is the Euler-Mascheroni constant and zeta(3) is Apéry's constant.

ExtremeValueDistribution

The distribution of X^(<1>) taken from a continuous uniform distribution over the unit interval has probability function

 P_N(x)=Nx^(N-1),
(7)

and distribution function

 D_N(x)=x^N.
(8)

The kth raw moment is given by

 mu_k^'=N/(k+N).
(9)

The first few central moments are

mu_2=N/((N+1)^2(N+2))
(10)
mu_3=-(2N(N-1))/((N+1)^3(N^2+5N+6))
(11)
mu_4=(N(9N^2-3N+6))/((N+1)^4(N+2)(N+3)(N+4)).
(12)

The mean, variance, skewness, and kurtosis excess are therefore given by

mu=N/(N+1)
(13)
sigma^2=N/((N+1)^2(N+2))
(14)
gamma_1=-(2(N-1))/(N+3)sqrt((N+2)/N)
(15)
gamma_2=(6(N^3-N^2-6N+2))/(N(N+3)(N+4)).
(16)

If X_i are instead taken from a standard normal distribution, then the corresponding cumulative distribution is

F(x)=1/(sqrt(2pi))int_(-infty)^xe^(-t^2/2)dt
(17)
=1/2+Phi(x),
(18)

where Phi(x) is the normal distribution function. The probability distribution of X^(<1>) is then

P(M_n<x)=[F(x)]^n
(19)
=n/(sqrt(2pi))int_(-infty)^x[F(t)]^(n-1)e^(-t^2/2)dt.
(20)

The mean mu(n) and variance sigma^2(n) are expressible in closed form for small n,

mu(1)=0
(21)
mu(2)=1/(sqrt(pi))
(22)
mu(3)=3/(2sqrt(pi))
(23)
mu(4)=3/(2sqrt(pi))[1+2/pisin^(-1)(1/3)]
(24)
mu(5)=5/(4sqrt(pi))[1+6/pisin^(-1)(1/3)]
(25)

and

sigma^2(1)=1
(26)
sigma^2(2)=1-1/pi
(27)
sigma^2(3)=(4pi-9+2sqrt(3))/(4pi)
(28)
sigma^2(4)=1+(sqrt(3))/pi-mu^2(4)
(29)
sigma^2(5)=1+(5sqrt(3))/(4pi)+(5sqrt(3))/(2pi^2)sin^(-1)(1/4)-mu^2(5).
(30)

No exact expression is known for mu(6) or sigma^2(6), but there is an equation connecting them

 mu^2(6)+sigma^2(6)=1+(5sqrt(3))/(4pi)+(15sqrt(3))/(2pi^2)sin^(-1)(1/4).
(31)

Wolfram Web Resources

Mathematica »

The #1 tool for creating Demonstrations and anything technical.

Wolfram|Alpha »

Explore anything with the first computational knowledge engine.

Wolfram Demonstrations Project »

Explore thousands of free applications across science, mathematics, engineering, technology, business, art, finance, social sciences, and more.

Computerbasedmath.org »

Join the initiative for modernizing math education.

Online Integral Calculator »

Solve integrals with Wolfram|Alpha.

Step-by-step Solutions »

Walk through homework problems step-by-step from beginning to end. Hints help you try the next step on your own.

Wolfram Problem Generator »

Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.

Wolfram Education Portal »

Collection of teaching and learning tools built by Wolfram education experts: dynamic textbook, lesson plans, widgets, interactive Demonstrations, and more.

Wolfram Language »

Knowledge-based programming for everyone.