R/prplaplace.R

In susanca/powdist: Power and Reversal Power Distributions

#' @name ReversalPowerLaplace
#' @examples
#' prplaplace(1, 1, 3, 4)
#' @export
prplaplace <- function(q, lambda = 1, mu = 0, sigma = 1, lower.tail = TRUE, log.p = FALSE){
  p = 1 - ( plaplace((-q+mu)/sigma)**lambda )
  if (lower.tail == FALSE) {
    p = 1 - p
  }
  if (log.p == TRUE) {
    p = log(p)
  }
  return(p)
}