context("mass-negbin")
skip_if_not_installed("modeltests")
library(modeltests)
skip_if_not_installed("MASS")
library(MASS)
fit <- glm.nb(Days ~ Sex/(Age + Eth*Lrn), data = MASS::quine)
test_that("MASS::glm.nb tidier arguments", {
check_arguments(tidy.negbin)
check_arguments(glance.negbin)
})
test_that("tidy.negbin", {
td1 <- tidy(fit)
td2 <- tidy(fit, conf.int = TRUE)
td3 <- tidy(fit, exponentiate = TRUE, conf.int = TRUE)
check_tidy_output(td2)
check_tidy_output(td1)
expect_false(NA %in% td2$conf.low)
expect_false(NA %in% td2$conf.high)
# exponentiate arg check
expect_equal(
as.matrix(td3[, c("estimate", "conf.low", "conf.high")]),
exp(as.matrix(td2[, c("estimate", "conf.low", "conf.high")]))
)
})
test_that("glance.negbin", {
gl <- glance(fit)
check_glance_outputs(gl)
check_dims(gl, 1, 8)
})
Generated by dwww version 1.15 on Thu Jun 20 16:43:45 CEST 2024.