tests/testthat/test-empinf.R

In marklhc/bootmlm: Bootstrap Resampling for Multilevel Models

context("Empirical influence function")

library(lme4)

m1 <- lmer(Yield ~ (1 | Batch), Dyestuff, REML = FALSE)

test_that("empinf_mer() and empinf_merm() give same results", {
  L1 <- empinf_mer(m1, sigma)
  L2 <- empinf_merm(m1, function(x) c(x@theta, sigma(x)))

  expect_length(L1, lme4::ngrps(m1))
  expect_equal(ncol(L2), 2)
  expect_identical(L1, L2[ , 2])
})