Function Reference: evcdf

statistics: p = evcdf (x)
statistics: p = evcdf (x, mu)
statistics: p = evcdf (x, mu, sigma)
statistics: p = evcdf (…, "upper")
statistics: [p, plo, pup] = evcdf (x, mu, sigma, pcov)
statistics: [p, plo, pup] = evcdf (x, mu, sigma, pcov, alpha)
statistics: [p, plo, pup] = evcdf (…, "upper")

Extreme value cumulative distribution function (CDF).

For each element of x, compute the cumulative distribution function (CDF) of the extreme value distribution (also known as the Gumbel or the type I generalized extreme value distribution) at the values in x with location parameter mu and scale parameter sigma. The size of p is the common size of x, 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. [p, plo, pup], evcdf computes the confidence bounds for p 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. plo and pup are arrays of the same size as p containing the lower and upper confidence bounds.

[…] = evcdf (…, "upper") computes the upper tail probability of the extreme value distribution with parameters x0 and gamma, at the values in x.

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 CDF, gumbelcdf.

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

See also: evinv, evpdf, evrnd, evfit, evlike, evstat, gumbelcdf

Source Code: evcdf

Example: 1

 

 ## Plot various CDFs from the extreme value distribution
 x = -10:0.01:10;
 p1 = evcdf (x, 0.5, 2);
 p2 = evcdf (x, 1.0, 2);
 p3 = evcdf (x, 1.5, 3);
 p4 = evcdf (x, 3.0, 4);
 plot (x, p1, "-b", x, p2, "-g", x, p3, "-r", x, p4, "-c")
 grid on
 legend ({"μ = 0.5, σ = 2", "μ = 1.0, σ = 2", ...
          "μ = 1.5, σ = 3", "μ = 3.0, σ = 4"}, "location", "southeast")
 title ("Extreme value CDF")
 xlabel ("values in x")
 ylabel ("probability")

                    
plotted figure