R/quaray.R

In wasquith/lmomco: L-Moments, Censored L-Moments, Trimmed L-Moments, L-Comoments, and Many Distributions

"quaray" <-
function(f,para,paracheck=TRUE) {
   if(! check.fs(f)) return()
   if(paracheck == TRUE) {
     if(! are.parray.valid(para)) return()
   }
   U <- para$para[1]
   A <- para$para[2]

   x <- U + sqrt(-2*A^2 * log(1-f))
   names(x) <- NULL
   return(x)
}