I do not understand why aov and anova deliver different results if I consider the same linear model. I first do this with aov function:
res.aov &
