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

X <- model.matrix(~ Text1 + Text2 + I(Text1^2) + I(Text2^2) + I(Text1*Text2), Crèmes)
round(X,3)
round(cor(X),3)

## Avec toutes les interactions :
summary(lm(Visc.t0 ~ Text1 + Text2 + I(Text1^2) + I(Text2^2) + I(Text1*Text2), data = Crèmes))

TabCoeff <- TabProba <- matrix(NA,6,6)
for (i in 1:6){
  res <- summary(lm(Crèmes[,i+2] ~ Text1 + Text2 + I(Text1^2) + I(Text2^2) + I(Text1*Text2), data = Crèmes))
  TabProba[,i] <- res$coeff[,4]
  TabCoeff[,i] <- res$coeff[,1]
}
round(TabProba,3)
round(TabCoeff,3)

par(mfrow=c(1,2))
plot(Visc.t0 ~ Text1, xlab="Texturant 1", ylab="Viscosité à t0", data=Crèmes,pch=20)
text(Crèmes[,c(1,3)],rownames(Crèmes),pos=4)
plot(Visc.t60 ~ Text2, xlab="Texturant 1", ylab="Viscosité à t60", data=Crèmes,pch=20)
text(Crèmes[,c(1,6)],rownames(Crèmes),pos=4)

par(mfrow=c(1,2))
plot(Visc.t0 ~ Text2, xlab="Texturant 2", ylab="Viscosité à t0", data=Crèmes,pch=20)
text(Crèmes[,c(2,3)],rownames(Crèmes),pos=4)
plot(Visc.t60 ~ Text2, xlab="Texturant 2", ylab="Viscosité à t60", data=Crèmes,pch=20)
text(Crèmes[,c(2,6)],rownames(Crèmes),pos=4)