Nothing
skip_on_cran()
skip_if_not_installed("lme4")
data("cbpp", package = "lme4")
set.seed(123)
model <- lme4::glmer(
cbind(incidence, size - incidence) ~ period + (1 | herd),
data = cbpp,
family = binomial(),
nAGQ = 0
)
params <- model_parameters(model, effects = "fixed")
test_that("model_parameters.glmer", {
expect_equal(params$SE, c(0.22758, 0.30329, 0.32351, 0.42445), tolerance = 1e-2)
})
test_that("print model_parameters", {
skip_if_not_installed("withr")
skip_if_not_installed("merDeriv")
withr::local_options(
list(
parameters_exponentiate = TRUE,
parameters_warning_exponentiate = TRUE
)
)
expect_snapshot(params)
suppressMessages({
mp <- model_parameters(model, effects = "all", exponentiate = TRUE)
})
expect_snapshot(mp)
set.seed(123)
model <- lme4::glmer(
cbind(incidence, size - incidence) ~ period + (1 | herd),
data = cbpp,
family = binomial(),
nAGQ = 2
)
mp <- model_parameters(model, effects = "all")
expect_snapshot(mp)
})
test_that("model_parameters.glmer ml1", {
params <- model_parameters(model, ci_method = "ml1", effects = "fixed")
expect_equal(params$SE, c(0.22758, 0.30329, 0.32351, 0.42445), tolerance = 1e-2)
expect_equal(params$df, c(54, 54, 54, 54), tolerance = 1e-2)
})
test_that("model_parameters.glmer betwithin", {
params <- model_parameters(model, ci_method = "betwithin", effects = "fixed")
expect_equal(params$SE, c(0.23009, 0.30433, 0.32476, 0.42632), tolerance = 1e-2)
expect_equal(params$df, c(822, 822, 822, 822), tolerance = 1e-2)
})
set.seed(123)
cbpp$time <- runif(nrow(cbpp), 1, 4)
model <- lme4::glmer(
cbind(incidence, size - incidence) ~ period + time + (1 + time | herd),
data = cbpp,
family = binomial(),
nAGQ = 0
)
test_that("model_parameters.glmer", {
params <- model_parameters(model, effects = "fixed")
expect_equal(params$SE, c(0.66539, 0.36178, 0.36223, 0.45528, 0.2379), tolerance = 1e-2)
})
test_that("model_parameters.glmer ml1", {
params <- model_parameters(model, ci_method = "ml1", effects = "fixed")
expect_equal(params$SE, c(0.66539, 0.36178, 0.36223, 0.45528, 0.2379), tolerance = 1e-2)
expect_equal(params$df, c(53, 53, 53, 53, 53), tolerance = 1e-2)
})
test_that("model_parameters.glmer betwithin", {
params <- model_parameters(model, ci_method = "betwithin", effects = "fixed")
expect_equal(params$SE, c(0.66539, 0.36178, 0.36223, 0.45528, 0.2379), tolerance = 1e-2)
expect_equal(params$df, c(821, 821, 821, 821, 9), tolerance = 1e-2)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.