test_that(
  "amr_logistic fits a binary logistic regression",
  {

    dat <- data.frame(
      AMR = c(
        0, 1, 0, 1,
        0, 0, 1, 1,
        1, 0, 1, 0
      ),
      age = c(
        2, 5, 3, 7,
        4, 6, 8, 1,
        9, 3, 5, 7
      ),
      treatment = c(
        0, 1, 1, 0,
        0, 1, 0, 1,
        1, 0, 1, 0
      )
    )

    model <- amr_logistic(
      data = dat,
      formula = AMR ~ age + treatment
    )

    expect_s3_class(
      model,
      "glm"
    )

    expect_equal(
      family(model)$family,
      "binomial"
    )
  }
)