test_that("blr_test_lr throws an error when the model is not an object of class glm", { model1 <- glm( honcomp ~ 1, data = hsb2, family = binomial(link = "logit") ) model2 <- lm(write ~ 1, data = hsb2) expect_error(blr_test_lr(model1, model2)) expect_error(blr_test_lr(model2, model1)) }) test_that("output is as expected when reduced model is not specified", { model <- glm( honcomp ~ female + read + science, data = hsb2, family = binomial(link = "logit") ) actual <- round(blr_test_lr(model)$test_result[['lr_ratio']], 4) expected <- 71.0525 expect_equal(actual, expected) })