evinv
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 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
## 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") |