Parameter estimation for specific distributions by the method of L-moments
Computes the parameters of a probability distribution as a function of the L-moments. The following distributions are recognized:
pelexp |
exponential | |
pelgam |
gamma | |
pelgev |
generalized extreme-value | |
pelglo |
generalized logistic | |
pelgpa |
generalized Pareto | |
pelgno |
generalized normal | |
pelgum |
Gumbel (extreme-value type I) | |
pelkap |
kappa | |
pelln3 |
three-parameter lognormal | |
pelnor |
normal | |
pelpe3 |
Pearson type III | |
pelwak |
Wakeby | |
pelwei |
Weibull | |
pelexp(lmom) pelgam(lmom) pelgev(lmom) pelglo(lmom) pelgno(lmom) pelgpa(lmom, bound = NULL) pelgum(lmom) pelkap(lmom) pelln3(lmom, bound = NULL) pelnor(lmom) pelpe3(lmom) pelwak(lmom, bound = NULL, verbose = FALSE) pelwei(lmom, bound = NULL)
lmom |
Numeric vector containing the L-moments of the distribution or of a data sample. |
bound |
Lower bound of the distribution. If |
verbose |
Logical: whether to print a message when not all parameters of the distribution can be computed. |
Numerical methods and accuracy are as described in
Hosking (1996, pp. 10–11).
Exception:
if pelwak
is unable to fit a Wakeby distribution using all 5 L-moments,
it instead fits a generalized Pareto distribution to the first 3 L-moments.
(The corresponding routine in the LMOMENTS Fortran package
would attempt to fit a Wakeby distribution with lower bound zero.)
A numeric vector containing the parameters of the distribution.
J. R. M. Hosking jrmhosking@gmail.com
Hosking, J. R. M. (1996). Fortran routines for use with the method of L-moments, Version 3. Research Report RC20525, IBM Research Division, Yorktown Heights, N.Y.
pelp
for parameter estimation of a general distribution
specified by its cumulative distribution function or quantile function.
lmrexp
, etc., to compute the L-moments
of a distribution given its parameters.
For individual distributions, see their cumulative distribution functions:
# Sample L-moments of Ozone from the airquality data data(airquality) lmom <- samlmu(airquality$Ozone) # Fit a GEV distribution pelgev(lmom)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.