R/calibrate.R

In melt: Multiple Empirical Likelihood Tests

#' Calibrate critical value and p-value
#'
#' Calibrate critical value and p-value in [elt()].
#'
#' @param alpha A single numeric.
#' @param statistic A single numeric.
#' @param calibrate A single character.
#' @param p A single integer.
#' @param par A numeric vector.
#' @param object An object that inherits from class \linkS4class{EL}.
#' @param control An object of class of \linkS4class{ControlEL}.
#' @return A numeric vector of length two for the critical value and the
#'   p-value.
#' @noRd
calibrate <- function(calibrate, alpha, statistic, p, par, object, control) {
  switch(calibrate,
    "chisq" = {
      c(
        cv = qchisq(1 - alpha, df = p),
        pval = pchisq(statistic, df = p, lower.tail = FALSE)
      )
    },
    "boot" = {
      seed <- if (is.null(control@seed)) {
        sample.int(.Machine$integer.max, 1L)
      } else {
        control@seed
      }
      compute_bootstrap_calibration(
        alpha, statistic, control@b, seed, control@nthreads,
        getMethodEL(object), getData(object), par, coef(object),
        control@maxit_l, control@tol_l, control@th, getWeights(object)
      )
    },
    "f" = {
      n <- nrow(getData(object))
      c(
        cv = qf(1 - alpha, df1 = p, df2 = n - p) * (p * (n - 1)) / (n - p),
        pval = pf(statistic * (n - p) / (p * (n - 1)),
          df1 = p, df2 = n - p, lower.tail = FALSE
        )
      )
    }
  )
}