Categories &

Functions List

Function Reference: evinv

statistics: x = evinv (p)
statistics: x = evinv (p, mu)
statistics: x = evinv (p, mu, sigma)
statistics: [x, xlo, xup] = evinv (p, mu, sigma, pcov)
statistics: [x, xlo, xup] = evinv (p, mu, sigma, pcov, alpha)

Inverse of the extreme value cumulative distribution function (iCDF).

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

Default values are mu = 0 and sigma = 1.

When called with three output arguments, i.e. [x, xlo, xup], evinv computes the confidence bounds for x when the input parameters mu and sigma 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 minima. For modeling maxima, use the alternative Gumbel iCDF, gumbelinv.

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

See also: evcdf, evpdf, evrnd, evfit, evlike, evstat, gumbelinv

Source Code: evinv

Example: 1

 

 ## Plot various iCDFs from the extreme value distribution
 p = 0.001:0.001:0.999;
 x1 = evinv (p, 0.5, 2);
 x2 = evinv (p, 1.0, 2);
 x3 = evinv (p, 1.5, 3);
 x4 = evinv (p, 3.0, 4);
 plot (p, x1, "-b", p, x2, "-g", p, x3, "-r", p, x4, "-c")
 grid on
 ylim ([-10, 10])
 legend ({"μ = 0.5, σ = 2", "μ = 1.0, σ = 2", ...
          "μ = 1.5, σ = 3", "μ = 3.0, σ = 4"}, "location", "northwest")
 title ("Extreme value iCDF")
 xlabel ("probability")
 ylabel ("values in x")

                    
plotted figure