Don <- read.table("https://husson.github.io/livre_exo/chap5_ex3.csv", header = TRUE, sep=";", encoding = "latin1", stringsAsFactors = TRUE) boxplot(Temps~Traitement , data=Don) var.test(Temps~Traitement,data=Don,subset=(Traitement!="B")) summary(aov(Temps~Traitement , data=Don)) ## ou library(FactoMineR) AovSum(Temps~Traitement , data=Don)