Skip to contents

Given a sample of exceedances, compute the first four L-moments and use either the first two to obtain the scale and shape (default), or else use L-skewness and L-scale to compute the scale and shape of the generalized Pareto distribution

Usage

gpd.lmom(xdat, thresh, sorted = FALSE, Lskew = FALSE)

Arguments

xdat

[numeric] vector of observations

thresh

[numeric] optional threshold argument

sorted

[logical] if TRUE, observations are sorted in increasing order

Lskew

[logical]; if TRUE, shape is obtained from L-skewness rather than first two moments.

Value

a vector of length two with the scale and shape estimates