如何执行具有健壮的标准错误的TukeyHSD？

Question

我目前在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不允许使用群集健壮的标准错误。正确的做法是什么？