lmdummy <- lm(colhe ~ dummy1 + dummy2 + dummy3 , data = colheitaDummy)
## avalie o modelo 
summary(lmdummy)
anova(lmdummy)

• ajuste o modelo normal de anova

lmAnova <- lm(colhe~solo, data=colheita)

## avalie o modelo
summary(lmAnova)
anova(lmAnova)

• compare os coeficientes dos dois modelos