R/qBrXII.R

In LMoFit: Advanced L-Moment Fitting of Distributions

#' Quantile distribution function of BrXII distribution
#'
#' @param u non-exceedance probability
#' @param RP Return Period "don't use in case u is used"
#' @param para parameters as c(scale, shape1, shape2)
#'
#' @return Quantile value/s using the inverse of the cumulative distribution function.
#' @author Mohanad Zaghloul [aut, cre], Simon Michael Papalexiou [aut, ths], Amin Elshorbagy [aut, ths]
#' @export
#'
#' @examples
#' 
#' x <- qBrXII(u = 0.99, para = c(1, 10, 0.8))
#' x <- qBrXII(RP = 100, para = c(1, 10, 0.8))
#' 
qBrXII <- function(u = NULL, RP = 1/(1 - u), para){
  if (is.null(u) & length(RP) >= 1) {u <- 1 - 1/RP}
  scale <- para[1]; shape_1 <- para[2]; shape_2 <- para[3]
  x <- scale*((((1 - u)^(-1*shape_1*shape_2)) - 1)/shape_2)^(1/shape_1) #BrXII quantile distribution function
  return(x)
}