don <- read.table("https://husson.github.io/livre_exo/chap10_ex9.csv", header = TRUE, sep=";",row.names=1, encoding = "latin1")

données <- cbind.data.frame(Acidité = as.vector(as.matrix(don[1:15,1:2])), Cocktail = as.factor(rep(rownames(don)[1:15],2)))
summary(aov(Acidité ~ Cocktail,data=données))

données2 <- cbind.data.frame(Acidité = apply(don[,1:2],1,mean), don[,3:6])
summary(lm(Acidité ~ -1+Orange+Banane+Mangue+Citron,data=données2[1:15,]))