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)