Flash Sale | 50% off

Last chance! 50% off unlimited learning

Sale ends in


lmomco (version 1.6.1)

lmomrevgum: L-moments of the Reverse Gumbel Distribution

Description

This function estimates the L-moments of the Reverse Gumbel distribution given the parameters ($\xi$ and $\alpha$) from parrevgum. The type-B L-moments in terms of the parameters are

λ1B=ξ(0.5722)αα{Ei(log(1ζ))}, λ2B=α{log(2)+Ei(2log(1ζ))Ei(log(1ζ))}, τ3=, τ4=, and τ5=.

where $\zeta$ is the right-tail censoring fraction of the sample or the nonexceedance probability of the right-tail censoring threshold, and $\mathrm{Ei}(x)$ is the exponential integral defined as

Ei(X)=Xx1exdx,

where $\mathrm{Ei}(-\log(1-\zeta)) \rightarrow 0$ as $\zeta \rightarrow 1$ and $\mathrm{Ei}(-\log(1-\zeta))$ can not be evaluated as $\zeta \rightarrow 0$.

Usage

lmomrevgum(para)

Arguments

para
The parameters of the distribution.

Value

  • An R list is returned.
  • L1Arithmetic mean.
  • L2L-scale---analogous to standard deviation.
  • LCVcoefficient of L-variation---analogous to coe. of variation.
  • TAU3The third L-moment ratio or L-skew--analogous to skew.
  • TAU4The fourth L-moment ratio or L-kurtosis---analogous to kurtosis.
  • TAU5The fifth L-moment ratio.
  • L3The third L-moment.
  • L4The fourth L-moment.
  • L5The fifth L-moment.
  • zetaNumber of samples observed (noncensored) divided by the total number of samples.
  • sourceAn attribute identifying the computational source of the L-moments: lmomrevgum.

References

Hosking, J.R.M., 1990, L-moments---Analysis and estimation of distributions using linear combinations of order statistics: Journal of the Royal Statistical Society, Series B, vol. 52, p. 105--124.

Hosking, J.R.M., 1995, The use of L-moments in the analysis of censored data, in Recent Advances in Life-Testing and Reliability, edited by N. Balakrishnan, chapter 29, CRC Press, Boca Raton, Fla., pp. 546--560.

See Also

parrevgum, quarevgum, cdfrevgum

Examples

Run this code
lmr <- lmom.ub(c(123,34,4,654,37,78))
rev.para <- lmom2par(lmr,type='revgum')
lmomrevgum(rev.para)

Run the code above in your browser using DataLab