Cacher is the code snippet organizer for pro developers

#M1: random match firm + same day
treat = which(df$IsBirthday == 1)
sic_i = df$SICCD[treat[110]]%/%100
date_i = df$date[treat[110]]

sample <- df %>%
filter(SIC2 == sic_i & date == date_i)

m1 <- glm(IsBirthday ~ Ln_ME_lag1 + Q_lag1 + Leverage_lag1 + ROA_lag1 + PPE_lag1 + Age + logTenure,
family=binomial(link = "probit"), data = sample)

#M3: random match firm + random match day