Poulpe <- read.table("https://husson.github.io/livre_exo/chap4_ex1.csv", header = TRUE, sep=";", encoding = "latin1", stringsAsFactors = TRUE)
by(Don[,"Poids"],Don[,"Sexe"],var)

var.test(Poids ~ Sexe, data = Poulpe) 

1 / var.test(Poids ~ Sexe, data = Poulpe)$statistic   ## On retrouve bien la même valeur de Fobs si on calcule la plus grande variance sur la plus petite