Categories &

Functions List

Function Reference: gumbelinv

statistics: x = gumbelinv (p)
statistics: x = gumbelinv (p, mu)
statistics: x = gumbelinv (p, mu, beta)
statistics: [x, xlo, xup] = gumbelinv (p, mu, beta, pcov)
statistics: [x, xlo, xup] = gumbelinv (p, mu, beta, pcov, alpha)

Inverse of the Gumbel cumulative distribution function (iCDF).

For each element of p, compute the quantile (the inverse of the CDF) of the Gumbel distribution (also known as the extreme value or the type I generalized extreme value distribution) with location parameter mu and scale parameter beta. The size of x is the common size of p, mu and beta. A scalar input functions as a constant matrix of the same size as the other inputs.

Default values are mu = 0 and beta = 1.

When called with three output arguments, i.e. [x, xlo, xup], gumbelinv computes the confidence bounds for x when the input parameters mu and beta are estimates. In such case, pcov, a 2×2 matrix containing the covariance matrix of the estimated parameters, is necessary. Optionally, alpha, which has a default value of 0.05, specifies the 100 * (1 - alpha) percent confidence bounds. xlo and xup are arrays of the same size as x containing the lower and upper confidence bounds.

The Gumbel distribution is used to model the distribution of the maximum (or the minimum) of a number of samples of various distributions. This version is suitable for modeling maxima. For modeling minima, use the alternative extreme value iCDF, evinv.

Further information about the Gumbel distribution can be found at https://en.wikipedia.org/wiki/Gumbel_distribution

Input arguments must be double or single; integer, logical, and character arrays are rejected. MATLAB accepts a character array and evaluates it at the character codes, which Octave deliberately does not, since a character array is an integer type and integers are refused too.

See also: gumbelcdf, gumbelpdf, gumbelrnd, gumbelfit, gumbellike, gumbelstat, evinv

Source Code: gumbelinv

Plot various iCDFs from the Gumbel distribution

 p = 0.001:0.001:0.999;
 x1 = gumbelinv (p, 0.5, 2);
 x2 = gumbelinv (p, 1.0, 2);
 x3 = gumbelinv (p, 1.5, 3);
 x4 = gumbelinv (p, 3.0, 4);
 plot (p, x1, '-b', p, x2, '-g', p, x3, '-r', p, x4, '-c')
 grid on
 ylim ([-5, 20])
 legend ({'μ = 0.5, β = 2', 'μ = 1.0, β = 2', ...
          'μ = 1.5, β = 3', 'μ = 3.0, β = 4'}, 'location', 'northwest')
 title ('Gumbel iCDF')
 xlabel ('probability')
 ylabel ('values in x')
plotted figure