tests/testthat/test-model_performance.rma.R

In performance: Assessment of Regression Models Performance

skip_if_not_installed("metafor")
skip_if_not_installed("metadat")
data(dat.bcg, package = "metadat")
dat <- metafor::escalc(
  measure = "RR",
  ai = tpos,
  bi = tneg,
  ci = cpos,
  di = cneg,
  data = dat.bcg
)

test_that("model_performance.rma", {
  model <- metafor::rma(yi, vi, data = dat, method = "REML")
  mp <- model_performance(model, estimator = "REML")
  expect_null(mp$R2)
  expect_equal(mp$AIC, 28.40474, tolerance = 1e-3)
  expect_equal(mp$I2, 0.9222139, tolerance = 1e-3)
  expect_identical(
    colnames(mp),
    c(
      "AIC",
      "BIC",
      "I2",
      "H2",
      "TAU2",
      "CochransQ",
      "p_CochransQ",
      "df_error",
      "Omnibus",
      "p_Omnibus"
    )
  )
})

test_that("check_outliers.rma", {
  model <- metafor::rma(yi, vi, data = dat, slab = author, method = "REML")
  out <- check_outliers(model)
  expect_equal(
    out,
    c(
      FALSE,
      FALSE,
      FALSE,
      TRUE,
      FALSE,
      FALSE,
      FALSE,
      TRUE,
      FALSE,
      FALSE,
      FALSE,
      FALSE,
      FALSE
    ),
    ignore_attr = TRUE
  )
  expect_identical(attributes(out)$outlier_var, c("Hart & Sutherland", "TPT Madras"))
  expect_snapshot(out)
})

test_that("model_performance.rma", {
  model <- metafor::rma(yi, vi, mods = cbind(ablat, year), data = dat)
  mp <- model_performance(model, estimator = "REML")
  expect_equal(mp$R2, 0.6463217, tolerance = 1e-3)
  expect_equal(mp$AIC, 24.21375, tolerance = 1e-3)
  expect_equal(mp$I2, 0.719778, tolerance = 1e-3)
  expect_identical(
    colnames(mp),
    c(
      "AIC",
      "BIC",
      "I2",
      "H2",
      "TAU2",
      "CochransQ",
      "p_CochransQ",
      "df_error",
      "Omnibus",
      "p_Omnibus",
      "R2"
    )
  )
})