ATTandCACE

ATTandCACE.R

library(DeclareDesign)
library(tidyverse)

design <-
  declare_population(N = 50,
                     type = sample(c("complier", "nevertaker"), N, replace = TRUE)) +
  declare_potential_outcomes(D ~ Z * (type == "complier")) +
  declare_potential_outcomes(Y ~ 5 * Z + rnorm(N)) +
  declare_assignment(prob = 0.5) +
  declare_estimand(CACE = mean(Y_Z_1 - Y_Z_0),
                   subset = type == "complier") +
  declare_estimand(ATT = mean(Y_Z_1 - Y_Z_0),
                   subset = D == 1)

# Different in a single draw
draw_estimands(design)

# Same in expectation
simulate_design(design, sims = 500) %>%
  group_by(estimand_label) %>%
  summarise(mean_estimand = mean(estimand))