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")) |
