gevfit

Generalized extreme value parameter estimates

collapse all in page

Syntax

pHat = gevfit(x)
[pHat,pCI] = gevfit(x)
[pHat,pCI] = gevfit(x,alpha)
[___] = gevfit(x,alpha,options)

Description

pHat = gevfit(x) returns the maximum likelihood estimates (MLEs) of the generalized extreme value distribution parameters (shape, scale, and location), given the sample data in x.

[pHat,pCI] = gevfit(x) also returns the 95% confidence intervals for the parameter estimates.

[pHat,pCI] = gevfit(x,alpha) specifies the confidence level for the confidence intervals to be 100(1 – alpha)%.

example

[___] = gevfit(x,alpha,options) specifies optimization options for the iterative algorithm gevfit uses to compute the MLEs. Create options by using the function statset.

You can specify [] for alpha to use its default value of 0.05.

Examples

collapse all

Open Live Script

Generate 100 random numbers from the generalized extreme value distribution with the shape parameter k=1, scale parameter sigma=4, and location parameter mu=3.

rng(0,"twister") % For reproducibility
k=1;
sigma = 4;
mu = 3;
x = gevrnd(k,sigma,mu,100,1);

Find the maximum likelihood estimates and the 99% confidence intervals.

[pHat,pCI] = gevfit(x,0.01)
pHat = 1×3

    1.0816    4.3981    3.2821


pCI = 2×3

    0.7102    2.9074    1.9684
    1.4531    6.6530    4.5957

pHat contains the MLE values for the shape, scale, and location parameters, in that order. pCI contains the 99% confidence intervals of each parameter. The values in the first row are the lower bounds, and the values in the second row are the upper bounds.

Input Arguments

collapse all

Sample data, specified as a numeric vector.

Data Types: single | double

Significance level for the confidence intervals, specified as a scalar in the range [0,1]. The confidence level is 100(1 – alpha)%, where alpha is the probability that the confidence intervals do not contain the true value. You can specify [] for alpha to use its default value of 0.05.

Data Types: single | double

Optimization options, specified as a structure. options determines the control parameters for the iterative algorithm used by gevfit to compute MLEs for censored data.

Create options by using the function statset or by creating a structure array containing the fields and values described in this table.

Field NameValueDefault Value
Display

Amount of information displayed by the algorithm

  • "off" — Displays no information

  • "final" — Displays the final output

  • "notify" — Displays output only if the function does not converge

"off"
MaxFunEvals

Maximum number of objective function evaluations allowed, specified as a positive integer

400
MaxIter

Maximum number of iterations allowed, specified as a positive integer

200
TolBnd

Lower bound of the standard deviation parameter estimate, specified as a positive scalar.

The bounds for the mean and standard deviation parameter estimates are [–Inf,Inf] and [TolBnd,Inf], respectively.

1e-6
TolFun

Termination tolerance for the objective function value, specified as a positive scalar

1e-6
TolX

Termination tolerance for the parameters, specified as a positive scalar

1e-6
OutputFcn

Specify one or more user-defined functions that an optimization function calls at each iteration, either as a function handle or as a cell array of function handles. For more information, see Optimization Solver Output Functions.

[]

You can also enter statset("gevfit") in the Command Window to see the names and default values of the fields included in the options structure.

Example: statset(Display="final",MaxIter=1000) specifies to display the final information of the iterative algorithm results, and change the maximum number of iterations allowed to 1000.

Data Types: struct

Output Arguments

collapse all

Estimates of the parameters k (shape), sigma (scale), and mu (location) of the generalized extreme value distribution, returned as a 1-by-3 numeric row vector.

Confidence intervals for the parameters of the generalized extreme value distribution, returned as a 2-by-3 numeric matrix containing the lower and upper bounds of the 100(1 – alpha)% confidence interval.

The first and second rows correspond to the lower and upper bounds of the confidence intervals, respectively.

Alternative Functionality

gevfit is a function specific to the generalized extreme value distribution. Statistics and Machine Learning Toolbox™ also offers the generic functions mle, fitdist, and paramci and the Distribution Fitter app, which support various probability distributions.

  • mle returns MLEs and the confidence intervals of MLEs for the parameters of various probability distributions. You can specify the probability distribution name or a custom probability density function.

  • Create a GeneralizedExtremeValueDistribution probability distribution object by fitting the distribution to data using the fitdist function or the Distribution Fitter app. The object properties k, sigma and mu store the parameter estimates. To obtain the confidence intervals for the parameter estimates, pass the object to paramci.

References

[1] Embrechts, P., C. Klüppelberg, and T. Mikosch. Modelling Extremal Events for Insurance and Finance. New York: Springer, 1997.

[2] Kotz, S., and S. Nadarajah. Extreme Value Distributions: Theory and Applications. London: Imperial College Press, 2000.

Extended Capabilities

expand all

Version History

Introduced before R2006a

See Also

| | | | | |

Topics