R/mbco_parametric.R

In RMediation: Mediation Analysis Confidence Intervals

mbco_parametric <- function(h0 = NULL,
                            h1 = NULL,
                            R = 10L,
                            alpha = .05,
                            checkHess = "No",
                            checkSE = "No",
                            optim = "NPSOL",
                            precision = 1e-9) {
  OpenMx::mxOption(NULL, "Calculate Hessian", checkHess)
  OpenMx::mxOption(NULL, "Standard Errors", checkSE)
  OpenMx::mxOption(NULL, "Function precision", precision)
  OpenMx::mxOption(NULL, "Default optimizer", optim)

  res <- OpenMx::mxCompare(h1,
                           h0,
                           boot = TRUE,
                           replications = R)
  mbco_chisq <- res$diffLL[2]  # Chi-square
  mbco_pvalue <- res$p[2]   # p-value
  mbco_df <- res$df[2] - res$df[1] # df
  mbcoTest <-
    list(chisq = mbco_chisq,
         df = mbco_df,
         p = mbco_pvalue)
  return(mbcoTest)
}