R/f2f.R

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

"f2f" <- function(f, pp=NA, xlo=NULL, type=c("ge", "gt")) {
   if(! check.fs(f)) return(FALSE)
   if(! is.null(xlo)) pp <- xlo$pp
   if(is.na(pp)) {
      warning("pp can not be NA")
      return(FALSE)
   } else {
      if(pp < 0 || pp > 1) {
        print("pp argument is not a valid nonexceedance probability")
        return(FALSE)
      }
   }
   type <- match.arg(type)
   if(type == "gt") {
     f <- f[f >  pp] # see identifical subsetting in f2flo
   } else {
     f <- f[f >= pp] # see identifical subsetting in f2flo
   }
   return(f)
}