I’m wondering if there is any way to get lm_robust objects to work with ANOVA. Alternatively, how would I go about trying to achieve the following:
First I run the following:
fit1 <- lm_robust(outcome ~ treatment / moderator - 1,
cluster = teamID,
se = "stata",
data = data)
fit2 <- lm_robust(outcome ~ treatment - 1,
cluster = teamID,
se = "stata",
data = data)
Now I want to see if there is an overall effect of the moderator by comparing fit1 with fit2 in an ANOVA.
anova(rescue_events, rescue_events2)
However, I get the error that lm_robust objects don’t work with anova. Is there any way to make them compatible?