Lmoments provides the estimate of L-moments of a sample or regional L-moments of a region.
1 2 3 4 5
vector representing a data-sample (or data from many samples defined with
array that defines the data subdivision among sites
The estimation of L-moments is based on a sample of size n, arranged in ascending order. Let x(1:n) <= x(2:n) <= ... <= x(n:n) be the ordered sample. An unbiased estimator of the probability weighted moments βr is:
br = 1/n sum[j from r+1 to n](x(j:n) (j-1)(j-2)...(j-r)/(n-1)/(n-2)/.../(n-r))
The sample L-moments are defined by:
l1 = b0
l2 = 2b1 - b0
l3 = 6b2 - 6b1 + b0
l4 = 20b3 - 30b2 + 12b1 - b0
and in general
l(r+1) = sum[k from 0 to r](b_k (-1)^(r-k) (r+k)! / (k!)^2 / (r-k)!)
where r=0, 1, ..., n-1.
The sample L-moment ratios are defined by
tr = lr / l2
and the sample L-CV by
t = l2 / l1
Sample regional L-CV, L-skewness and L-kurtosis coefficients are defined as
t^R = sum[i from 1 to k](ni t^(i)) / sum[i from 1 to k](ni)
t3^R = sum[i from 1 to k](ni t3^(i)) / sum[i from 1 to k](ni)
t4^R = sum[i from 1 to k](ni t4^(i)) / sum[i from 1 to k](ni)
Lmoments gives the L-moments (l_1, l_2, t, t_3, t_4),
regionalLmoments gives the regional weighted L-moments (l_1^R, l_2^R, t^R, t_3^R, t_4^R),
LCV gives the coefficient of L-variation,
LCA gives the L-skewness and
Lkur gives the L-kurtosis of
For information on the package and the Author, and for all the references, see
1 2 3 4 5 6 7 8 9 10 11 12 13 14
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.