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é)