gplike
Generalized Pareto negative loglikelihood
Syntax
nlogL = gplike(params,data)
[nlogL,acov] = gplike(params,data)
Description
nlogL = gplike(params,data)
returns the negative of the loglikelihoodnlogL
for the two-parameter generalized Pareto (GP) distribution, evaluated at parametersparams
.params(1)
is the tail index (shape) parameter,k
, andparams(2)
is the scale parameter.gplike
does not allow a threshold (location) parameter.
[nlogL,acov] = gplike(params,data)
returns the inverse of Fisher's information matrix,acov
. If the input parameter values inparams
are the maximum likelihood estimates, the diagonal elements ofacov
are their asymptotic variances.acov
is based on the observed Fisher's information, not the expected information.
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
,
.
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.