lmomlap {lmomco}R Documentation

L-moments of the Laplace Distribution

Description

This function estimates the L-moments of the Laplace distribution given the parameters (ξ\xi and α\alpha) from parlap. The L-moments in terms of the parameters are λ1=ξ\lambda_1 = \xi, λ2=3α/4\lambda_2 = 3\alpha/4, τ3=0\tau_3 = 0, τ4=17/22\tau_4 = 17/22, τ5=0\tau_5 = 0, and τ6=31/360\tau_6 = 31/360.

For rr odd and r3r \ge 3, λr=0\lambda_r = 0, and for rr even and r4r \ge 4, the L-moments using the hypergeometric function 2F1(){}_2F_1() are

λr=2αr(r1)[12F1(r,r1,1,1/2)]\mbox,\lambda_r = \frac{2\alpha}{r(r-1)}[1 - {}_2F_1(-r, r-1, 1, 1/2)]\mbox{,}

where 2F1(a,b,c,z){}_2F_1(a, b, c, z) is defined as

2F1(a,b,c,z)=n=0(a)n(b)n(c)nznn!\mbox,{}_2F_1(a, b, c, z) = \sum_{n=0}^\infty \frac{(a)_n(b)_n}{(c)_n}\frac{z^n}{n!}\mbox{,}

where (x)n(x)_n is the rising Pochhammer symbol, which is defined by

(x)n=1\mbox for n=0\mbox,and(x)_n = 1 \mbox{\ for\ } n = 0\mbox{, and}

(x)n=x(x+1)(x+n1)\mbox for n>0\mbox.(x)_n = x(x+1)\cdots(x+n-1) \mbox{\ for\ } n > 0\mbox{.}

Usage

lmomlap(para)

Arguments

para

The parameters of the distribution.

Value

An R list is returned.

lambdas

Vector of the L-moments. First element is λ1\lambda_1, second element is λ2\lambda_2, and so on.

ratios

Vector of the L-moment ratios. Second element is τ\tau, third element is τ3\tau_3 and so on.

trim

Level of symmetrical trimming used in the computation, which is 0.

leftrim

Level of left-tail trimming used in the computation, which is NULL.

rightrim

Level of right-tail trimming used in the computation, which is NULL.

source

An attribute identifying the computational source of the L-moments: “lmomlap”.

Author(s)

W.H. Asquith

References

Hosking, J.R.M., 1986, The theory of probability weighted moments: IBM Research Report RC12210, T.J. Watson Research Center, Yorktown Heights, New York.

See Also

parlap, cdflap, pdflap, qualap

Examples

lmr <- lmoms(c(123,34,4,654,37,78))
lmr
lmomlap(parlap(lmr))

[Package lmomco version 2.5.1 Index]