В настоящее время я использую команду lm_robust в библиотеке estimatr для генерации устойчивых к кластеру стандартных ошибок.Я также могу использовать felm для этой цели.
Я могу получить все попарные сравнения фактора в настоящее время, например, так:
aov_model <- aov(outcome~ treatment*condition + education + income + sex + birth + race + trust_daily + turk_experience, data = all_data)
res <- TukeyHSD(my_lm, which = c('treatment:condition'))
as.data.frame(tidy(res,my_lm)) %>%
filter(adj.p.value < .05)
Однако aov не допускает стандартных ошибок кластера.Какой правильный способ сделать это?