lmomgpaRC {lmomco} | R Documentation |
B-type L-moments of the Generalized Pareto Distribution with Right-Tail Censoring
Description
This function computes the “B”-type L-moments of the Generalized Pareto distribution given the parameters (,
, and
) from
pargpaRC
and the right-tail censoring fraction . The B-type L-moments in terms of the parameters are
where and
is the right-tail censor fraction or the probability
that
is less than the quantile at
nonexceedance probability: (
). In other words, if
, then there is no right-tail censoring. Finally, the
RC
in the function name is to denote R
ight-tail C
ensoring.
Usage
lmomgpaRC(para)
Arguments
para |
The parameters of the distribution. Note that if the |
Value
An R list
is returned.
lambdas |
Vector of the L-moments. First element is
|
ratios |
Vector of the L-moment ratios. Second element is
|
trim |
Level of symmetrical trimming used in the computation, which is |
leftrim |
Level of left-tail trimming used in the computation, which is |
rightrim |
Level of right-tail trimming used in the computation, which is |
source |
An attribute identifying the computational source of the L-moments: “lmomgpaRC”. |
message |
For clarity, this function adds the unusual message to an L-moment object that the |
zeta |
The censoring fraction. Assumed equal to unity if not present in the |
Author(s)
W.H. Asquith
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, v. 52, pp. 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
pargpa
, pargpaRC
, lmomgpa
, cdfgpa
, pdfgpa
, quagpa
Examples
para <- vec2par(c(1500,160,.3),type="gpa") # build a GPA parameter set
lmorph(lmomgpa(para))
lmomgpaRC(para) # zeta = 1 is internally assumed if not available
# The previous two commands should output the same parameter values from
# independent code bases.
# Now assume that we have the sample parameters, but the zeta is nonunity.
para$zeta = .8
lmomgpaRC(para) # The B-type L-moments.