I’ve often used linear regression to test if mean values differ between groups by dummy coding my categorical variable, which I think is basically the same thing (or at least I get the same results) as using ANOVA. I have used lm() function in R for doing this.
Previously, if my data did not meet the assumptions of linear regression, I’ve used data transformations. Sometimes this works better and sometimes not so good. As far as I’m concerned, I could use generalized linear models to compare group means for data that follow e.g. Poisson or negative binomial distributions without the need to transform data.
The problem is, when I fit the model and get the model summary (using glm() function in R), I don’t see the p-value for the full model – which I get in the last line of model summary when I am fitting linear models using lm() function. Model summary - when using glm() - gives me only the p and Z-values for each coefficient which I can use for pairwise comparison.
The main idea why I would like to get the p-value for full model is that I could use glm() as a substitute for ANOVA for data that does not meet its assumptions.
All help is much appreciated!
?anova.glm
[r] mixed models p-values