1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32
| # Chargement des données
library(MASS)
# Analyse discriminante factorielle
Dataset.lda <- lda(QUAL ~ ., data=Dataset)
print(Dataset.lda)
plot(Dataset.lda)
#Calcul des projections sur les axes factoriels
Dataset.pred <- predict(Dataset.lda, data=Dataset)
# Calcul des moyennes conditionnelles sur les axes
m <- matrix (rep(0,6), nrow=3, ncol = 2)
for (i in 1:3){
for (j in 1:2){
m[i,j] <- mean (Dataset.pred$x[unclass(Dataset$QUAL)--i,j])
}
}
#graphique - carte factorielle avec les moyennes conditionnelles des groupes (asp=1 pour que le positionnement relatif sur les axes soit respecté)
plot (Dataset.pred$x[,1], Dataset.pred$x[,2], main="Carte factorielle", xlab="Axe.1", ylab="Axe.2", pch=21, bg=c("red", "green", "blue")[unclass(Dataset$QUAL)], asp=1)
abline (0, 0, h=0)
abline (a=0, b=0, v=0)
text(m[,1],m[,2],labels=levels(Dataset$QUAL), col=c("red","green","blue")) |
Partager