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])
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.