# NOT RUN {
dat <- fabricate(
N = 10,
U = rnorm(N),
potential_outcomes(Y ~ 0.1 * Z + U)
)
fabricate(
data = dat,
Z = rbinom(N, 1, prob = 0.5),
Y = reveal_outcomes(Y ~ Z)
)
fabricate(
N = 10,
U = rnorm(N),
potential_outcomes(Y ~ 0.1 * Z1 + 0.3 * Z2 + 0.5 * Z1 * Z2 + U,
conditions = list(Z1 = c(0, 1),
Z2 = c(0, 1))),
Z1 = rbinom(N, 1, prob = 0.5),
Z2 = rbinom(N, 1, prob = 0.5),
Y = reveal_outcomes(Y ~ Z1 + Z2)
)
# }
Run the code above in your browser using DataLab