evalue

transform_to_rr <- function(estimate, type) {
  if (type == "or") estimate <- sqrt(estimate)
  if (type == "hr") {
    estimate <- (1 - 0.5^sqrt(estimate)) / (1 - 0.5^sqrt(1 / estimate))
  }
  if (estimate < 1) estimate <- 1 / estimate
    
  estimate
}

evalue <- function(estimate, type = c("rr", "or", "hr")) {
  # validate arguments
  if (!is.numeric(estimate)) stop("`estimate` must be numeric")
  type <- match.arg(type)

  # calculate evalue
  estimate <- transform_to_rr(estimate, type)
  estimate + sqrt(estimate * (estimate - 1))
}