parwak
are consistent with the corresponding distribution,
otherwise a list would not have been returned. However, other
functions (cdfwak
, quawak
, and
lmomwak
require consistent parameters to return the cumulative
probability (nonexceedance), quantile, and L-moments of the distribution,
respectively. These functions internally use the are.parwak.valid
function. The FORTRAN source code of Hosking provides the basis
for the function.are.parwak.valid(para,nowarn=FALSE)
parwak
.TRUE
then options(warn=-1)
is made and restored on return. This switch is to permit calls in which warnings are not desired as the user knows how to handle the returned value---say in awak
consistent.wak
consistent.Hosking, J.R.M. and Wallis, J.R., 1997, Regional frequency analysis---An approach based on L-moments: Cambridge University Press.
is.wak
para <- parwak(lmom.ub(c(123,34,4,654,37,78)))
if(are.parwak.valid(para)) Q <- quawak(0.5,para)
Run the code above in your browser using DataLab