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?