skip_on_cran()
skip_on_os(c("mac", "solaris"))
skip_if_not_installed("lme4")
skip_if_not_installed("glmmTMB")
skip_if_not_installed("marginaleffects")
# glmmTMB ----
test_that("ggpredict, glmmTMB prediction intervals random effects", {
data(iris)
m <- glmmTMB::glmmTMB(Sepal.Length ~ Sepal.Width + (1 | Species), data = iris)
out1 <- ggpredict(m, "Species", type = "random")
out2 <- ggpredict(m, "Species", type = "random", interval = "confidence")
out3 <- ggpredict(m, "Species", type = "random", verbose = FALSE)
out4 <- ggpredict(m, c("Sepal.Width", "Species"), type = "random")
expect_equal(out1$conf.low, c(3.85055, 5.29633, 5.78545), tolerance = 1e-3)
expect_equal(out2$conf.low, c(4.57423, 6.02764, 6.52735), tolerance = 1e-3)
expect_message(print(out1), regex = "prediction")
expect_no_message(print(out3))
expect_no_message(print(hypothesis_test(out1)))
expect_message(
print(hypothesis_test(out4)),
regex = "Intervals"
)
expect_no_message(
print(print(hypothesis_test(out1, verbose = FALSE)))
)
})
# lmer ----
test_that("ggpredict, lmer prediction intervals random effects", {
data(iris)
m <- lme4::lmer(Sepal.Length ~ Sepal.Width + (1 | Species), data = iris)
out1 <- ggpredict(m, "Species", type = "random")
# not sure why this works on devel again - maybe fix in Matrix?
if (getRversion() > "4.4.1") {
expect_equal(out1$conf.low, c(3.27188912101159, 4.72096548348341, 5.20851196398241), tolerance = 1e-4)
} else {
expect_null(out1$conf.low)
}
out1 <- ggpredict(m, c("Sepal.Width [2:4]", "Species"), type = "random")
out2 <- ggpredict(m, c("Sepal.Width [2:4]", "Species"), type = "random", interval = "confidence")
expect_equal(
out1$conf.low,
c(
2.41206, 3.86114, 4.34869, 3.22614, 4.67521, 5.16276, 4.00984,
5.45891, 5.94646
),
tolerance = 1e-3
)
expect_equal(
out2$conf.low,
c(
2.69634, 4.14542, 4.63296, 3.51457, 4.96365, 5.45119, 4.29495,
5.74403, 6.23158
),
tolerance = 1e-3
)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.