Generalized Pareto cumulative distribution function
p = gpcdf(x,k,sigma,theta)
p = gpcdf(x,k,sigma,theta,'upper')
p = gpcdf(x,k,sigma,theta)
returns the cdf of the generalized Pareto (GP) distribution with the tail index (shape) parameterk
, scale parametersigma
, and threshold (location) parameter,theta
, evaluated at the values inx
. The size ofp
is the common size of the input arguments. A scalar input functions as a constant matrix of the same size as the other inputs.
p = gpcdf(x,k,sigma,theta,'upper')
returns the complement of the cdf of the generalized Pareto (GP) distribution, using an algorithm that more accurately computes the extreme upper tail probabilities.
Default values fork
,sigma
, andtheta
are 0, 1, and 0, respectively.
Whenk = 0
andtheta = 0
, the GP is equivalent to the exponential distribution. Whenk > 0
andtheta = sigma/k
, the GP is equivalent to a Pareto distribution with a scale parameter equal tosigma/k
and a shape parameter equal to1/k
. The mean of the GP is not finite whenk
≥1
, and the variance is not finite whenk
≥1/2
. Whenk
≥0
, the GP has positive density for
x > theta
, or, when
k < 0
,
.
[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.