test_that("Default maximum iterations", {
  Bids$size.sq <- Bids$size ^ 2
  fit <- glm.hP(formula.mu = numbids ~ leglrest + rearest + finrest +
                 whtknght + bidprem + insthold + size + size.sq + regulatn,
                 formula.gamma = numbids ~ 1, data = Bids)
  expect_equal(fit$maxiter_series, 1000)
})