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

Crèmes <- cbind.data.frame(Crèmes[,1:2],Text1²=I(Crèmes[,1]^2),Text2²=I(Crèmes[,2]^2),Text1_Text2=I(Crèmes[,1]*Crèmes[,2]),Crèmes)
for (i in 6:7) Crèmes[,i]=as.factor(Crèmes[,i])
levels(Crèmes[,6]) <- paste("Text1=",levels(Crèmes[,6]),sep="")
levels(Crèmes[,7]) <- paste("Text2=",levels(Crèmes[,7]),sep="")

library(FactoMineR)
PCA(Crèmes, quali.sup=6:7,quanti.sup = 1:5)