fatal_fe_lm_mod <- lm(fatal_rate ~ beertax + state, data = USFatalities) fatal_fe_lm_mod summary(fatal_fe_lm_mod)