design_a <-
declare_model(N = 100,
U = rnorm(N),
Y_Z_0 = U,
Y_Z_1 = U + rnorm(N, mean = 2, sd = 2)) +
declare_inquiry(ATE = mean(Y_Z_1 - Y_Z_0)) +
declare_assignment(Z = complete_ra(N, prob = 0.5)) +
declare_measurement(Y = reveal_outcomes(Y ~ Z)) +
declare_estimator(Y ~ Z, inquiry = "ATE")
design_b <- replace_step(
design_a, step = "assignment",
declare_assignment(Z = complete_ra(N, prob = 0.3)) )
comparison <- compare_diagnoses(design_a, design_b, sims = 40)
Run the code above in your browser using DataLab