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
gpd.lmom(xdat, thresh, sorted = FALSE, Lskew = FALSE)a vector of length two with the scale and shape estimates
[numeric] vector of observations
[numeric] optional threshold argument
[logical] if TRUE, observations are sorted in increasing order
[logical]; if TRUE, shape is obtained from L-skewness rather than first two moments.