Orange <- read.table("https://husson.github.io/livre_exo/chap9_ex9.csv", header = TRUE, sep=";", encoding = "latin1", stringsAsFactors = TRUE)

summary(lm(Appréciation ~Sucrée,data=Orange,subset=Sexe=="F"))
summary(lm(Appréciation ~Sucrée,data=Orange,subset=Sexe=="H"))

library(FactoMineR)
AovSum(Appréciation ~Sucrée+Sexe+Sucrée:Sexe,data=Orange)