R/internal_functions.R

In MOsemiind: Marshall-Olkin Shock Models with Semi-Independent Time

computeCopula <- function(type, theta) {
  if (type == "Clayton") {
    return(copula::claytonCopula(param = theta, dim = 3))
  } else if (type == "Gumbel") {
    if (theta < 1) {
      stop("Error: For a Gumbel copula, theta must be greater than or equal to 1.", call. = FALSE)
    }
    return(copula::gumbelCopula(param = theta, dim = 3))
  } else {
    stop("Invalid copula type")
  }
}



computeDependenceMeasures <- function(copula) {
  list(kendall = 0.5, spearman = 0.6)  # Placeholder values
}