context("or_glm")
library("mgcv")
library("MASS")
test_that("correct level count of indicator variable for glm", {
data("data_glm")
fit_glm <- glm(admit ~ rank, data = data_glm, family = "binomial") # fit model
out <- or_glm(data = data_glm, model = fit_glm)
expect_length(out$predictor, length(levels(data_glm$rank)) - 1)
})
test_that("or_glm works with glmmPQL", {
data(bacteria)
fit_glmmpql <- glmmPQL(y ~ trt + week,
random = ~ 1 | ID,
family = binomial, data = bacteria,
verbose = FALSE
)
# Apply function
out <- expect_warning(or_glm(data = bacteria, model = fit_glmmpql, incr = list(week = 5)))
expect_length(out, 5)
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.