library(FactoMineR)
pollos <- read.table("http://factominer.free.fr/libra/pollos.csv",header=TRUE,sep=";",dec=".",row.names=1)
pollos <- as.data.frame(t(pollos))
regimen <- as.factor(c(rep("N",6),rep("J16",5),rep("J16R5",8),rep("J16R16",9),rep("J48",6),rep("J48R24",9)))
pollos <- cbind.data.frame(regimen,pollos)
colnames(pollos)[1] <- "Régimen"
res.pca <- PCA(pollos,quali.sup=1)
plot(res.pca, habillage=1, cex=0.7)
plot(res.pca, choix="var", invisible="var")
points(res.pca$var$coord[,1:2], cex=0.5)

dimdesc(res.pca,proba=1e-5)
plot(res.pca, habillage=1, axes=3:4)
plot(res.pca, choix="var", invisible="var", axes=3:4)
points(res.pca$var$coord[,3:4], cex=0.5)

aa <- cbind.data.frame(pollos[,1], res.pca$ind$coord[,1:2])
bb <- coord.ellipse(aa,bary=TRUE)
plot.PCA(res.pca, habillage=1, ellipse=bb)

aa <- cbind.data.frame(pollos[,1], res.pca$ind$coord[,3:4])
bb <- coord.ellipse(aa, bary=TRUE)
plot.PCA(res.pca, habillage=1, ellipse=bb, axes=3:4)