This function computes the quantiles of the Wakeby distribution given
parameters (\(\xi\), \(\alpha\), \(\beta\), \(\gamma\), and \(\delta\)) computed by parwak
. The quantile function is
$$x(F) = \xi+\frac{\alpha}{\beta}(1-(1-F)^\beta)-
\frac{\gamma}{\delta}(1-(1-F))^{-\delta} \mbox{,}$$
where \(x(F)\) is the quantile for nonexceedance probability \(F\), \(\xi\) is a location parameter, \(\alpha\) and \(\beta\) are scale parameters, and \(\gamma\) and \(\delta\) are shape parameters. The five returned parameters from parwak
in order are \(\xi\), \(\alpha\), \(\beta\), \(\gamma\), and \(\delta\).
quawak(f, wakpara, paracheck=TRUE)
Quantile value for nonexceedance probability \(F\).
Nonexceedance probability (\(0 \le F \le 1\)).
The parameters from parwak
or vec2par
.
A logical controlling whether the parameters are checked for validity. Overriding of this check might be extremely important and needed for use of the quantile function in the context of TL-moments with nonzero trimming.
W.H. Asquith
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., 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.
cdfwak
, pdfwak
, lmomwak
, parwak
lmr <- lmoms(c(123,34,4,654,37,78))
quawak(0.5,parwak(lmr))
Run the code above in your browser using DataLab