str(warpbreaks)
m.lm <- lm(breaks~wool*tension, data=warpbreaks)
m.glm <- glm(breaks~wool*tension, data=warpbreaks, family=gaussian)
anova(m.lm)
anova(m.glm, test="F")
str(swiss)
m.lm <- lm(Infant.Mortality~., data=swiss)
m.glm <- glm(Infant.Mortality~., data=swiss, family=gaussian)
anova(m.lm)
anova(m.glm, test="F")