result_sf <- causens_sf(Z ~ X.1 + X.2 + X.3, "Y",
data = data_sf, c1 = 0.25, c0 = 0.25, bootstrap = TRUE
test_that("summary.causens_sf produces correct output", {
"Estimate     Std.Error    95% C.I.             ",
N = 10000, seed = 123, y_type = "binary", alpha_uz = 1,
result_mc <- causens_monte_carlo("Y", "Z", c("X.1", "X.2", "X.3"),
test_that("summary.monte_carlo_causens produces correct output", {

"Estimate     Std.Error    95% C.I.             "