parkap
.
The quantile function of the distribution is$$x(F) = \xi + \frac{\alpha}{\kappa}\left(1-{\left(\frac{1-F^h}{h}\right)}^\kappa\right) \mbox{,}$$
where $x(F)$ is the quantile for nonexceedance probability $f$, $\xi$ is a location parameter, $\alpha$ is a scale parameter, $\kappa$ is a shape parameter, and $h$ is another shape parameter.
quakap(f, para, paracheck=TRUE)
parkap
or similar.Hosking, J.R.M., 1996, FORTRAN routines for use with the method of L-moments: Version 3, IBM Research Report RC20525, T.J. Watson Research Center, Yorktown Heights, New York.
Hosking, J.R.M. and Wallis, J.R., 1997, Regional frequency analysis---An approach based on L-moments: Cambridge University Press.
cdfkap
, parkap
lmr <- lmom.ub(c(123,34,4,654,37,78,21,32,231,23))
quakap(0.5,parkap(lmr))
Run the code above in your browser using DataLab