R/expexpcp.R

In ERPeq: Probabilistic Hazard Assessment

#' @title Probabilistic estimation of earthquake recurrence interval using exponentiated exponential distribution
#' @description Computes the probability of an earthquake within a specified time "r" and elapsed time "te".
#'
#' @param fit Fit is the fitexpexp object. See ?fitexpexp for details.
#' @param r The specified time in which the probability of an earthquake is desired to be calculated.
#' @param te Elapsed time since the last earthquake
#' @return A numeric value
#' @references Pasari, S. and Dikshit, O. (2014). Impact of three-parameter Weibull models in probabilistic assessment of earthquake hazards. Pure and Applied Geophysics, 171, 1251-1281.
#' @importFrom stats ks.test optim
#' @importFrom VGAM dpareto ppareto rpareto
#' @import methods
#' @export
#' @examples
#' fit=fitexpexp(c(1,1),data=data_earthquake_7)
#' expexpcp(fit,r=2,te=5)

expexpcp=function(fit,r,te){
  if (missingArg(fit) == TRUE)
    stop("fit object is missing. The fit object must be fitexpexp object. See ?fitexpexp for details.")
  if (missingArg(r) == TRUE)
    stop("r object is missing.")
  if (missingArg(te) == TRUE)
    stop("te object is missing.")

  GG=cdfeexp(fit$mle,r+te)
  G=cdfeexp(fit$mle,te)
  f=(GG-G)/(1-G)
  return(f)
}