Acidité <- read.table("https://husson.github.io/livre_exo/chap4_ex3.csv", header = TRUE, sep=";", encoding = "latin1", stringsAsFactors = TRUE) by(Acidité[,"Note"],Acidité[,"Jury"],var) var.test(Note ~ Jury, data = Acidité)