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