# NOT RUN {
# Simulate example data
simulated_data <- generate_simulated_data(n = 100)
# Fit by outcome modeling
estimate <- gapclosing(
data = simulated_data,
outcome_formula = formula(outcome ~ treatment * category + confounder),
treatment_name = "treatment",
category_name = "category",
counterfactual_assignments = 1
)
summary(estimate)
# Convert to a data frame
estimate.df <- as.data.frame(estimate)
# }
Run the code above in your browser using DataLab