library(tidyverse); library(modelsummary); library(fixest)
od <- causaldata::organ_donations
# Treatment variable
od <- od %>%
mutate(Treated = State == 'California' &
Quarter %in% c('Q32011','Q42011','Q12012'))
# feols clusters by the first
# fixed effect by default, no adjustment necessary
clfe <- feols(Rate ~ Treated | State + Quarter,
data = od)
msummary(clfe, stars = c('*' = .1, '**' = .05, '***' = .01),
coef = "Treated", nobs = TRUE)