test_that("bias works for guassian", {
  testthat::expect_snapshot_value(
    bias(residuals(fit_glm), treatment = glm_data$treatment, group_idx = list(seq_len(nrow(glm_data)))),
    style = "deparse",
    tolerance = 1e-4
  )
})