compare_binary takes in a variable e.g. union
and runs bivariate regression of x on treatment (for summary statistics)
Examples
dta <- data.frame(x = rnorm(100), treat = rep(c(0, 1), 50))
compare_binary("x", "treat", dta, report = "diff")
#> [1] 0.17410745 -0.03708638 0.00000000
