R/quaben.R

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

"quaben" <- function(f, para=list(para=c(1, 10)), ...) {
  m <- para$para[1]
  if(m > 9) {
    warning("expect failure as first parameter is an integer greater than 9")
  }
  b <- para$para[2]
  if(b != 10) {
    warning("pmfben is only configured for base10 at the moment in the second parameter")
    return(NULL)
  }
  s <- (1*b^(m-1))
  e <- as.integer(paste(rep(9, m), collapse=""))
  x <- cut(f, breaks=c(0, sapply(s:e, function(k) sum(pmfben(s:k, para=para)))),
                          labels=s:e, include.lowest=TRUE)
  as.integer( as.character(x) )
}