test_that("glm contrast error is captured", {
  ## mm has a factor variable with only one active level
  nb_data <- 50
  mm <- data.frame(
    x = rnorm(nb_data),
    z = rep(factor("a", levels = c("a", "b"))), nb_data
  )
  target <- sample(factor(0:1), nb_data, replace = TRUE)
  ## trigger an error with stats::glm
  expect_error(d_model <- stats::glm(target ~ ., data = mm, family = stats::binomial()))
  ## the error is captured by fit_glm
  expect_no_error(model <- fit_glm(target, mm, 2, covlmc_control()))
  ## returns a constant model
  expect_s3_class(model, "constant_model")
  expect_equal(model$rank, 0)
})