Acidité <- read.table("https://husson.github.io/livre_exo/chap5_ex2.csv", header = TRUE, sep=";", encoding = "latin1")
Acidité[,"Bière"] <- as.factor(Acidité[,"Bière"])
boxplot(Note~Bière , data=Acidité, ylab="Note d'acidité",xlab="Bière")
summary(aov(Note~Bière,data=Acidité))

## ou

library(FactoMineR)
AovSum(Note~Bière,data=Acidité)