-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathACP.R
39 lines (29 loc) · 936 Bytes
/
ACP.R
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
33
34
35
36
37
38
39
# Analyse en Composantes Principales
## Étape 1 - lecture du jeu de données :
### > Importer le jeu de données
attach(data)
rownames(data) = data$id
data = data[,-1]
## Étape 2 - description des données :
summary(data)
## Étape 3 - analyse du jeu de données :
### > Charger la librairie ade4
### Calcul de l'ACP :
acp = dudi.pca(data)
### > Sélectionner le nombre d'axes qui permet de représenter environ 70 % du jeu de données
### Choix du nombre d'axes :
inertia.dudi(acp)
acp = dudi.pca(data)
### Coordonnées dans le nouveau référentiel :
acp$co
acp$li
### Graphiques :
s.corcircle(acp$co,xax = 1,yax = 2)
s.label(acp$li,xax=1,yax=2)
s.arrow(3*acp$co,xax=1,yax=2,add.plot = TRUE)
## Étape 4 - étude des contributions :
inertia.dudi(acp,col=T,row=T)$col.abs
inertia.dudi(acp,col=T,row=T)$row.abs
## Étape 5 - étude des qualités :
inertia.dudi(acp,col=T,row=T)$col.rel
inertia.dudi(acp,col=T,row=T)$row.rel