Espace des individus ℝᵖ Prioritaire
Le nuage des n individus dans l'espace des variables, le principe d'ajustement, et la matrice à diagonaliser — cœur du calcul de l'ACP.
1. Le nuage des individus dans ℝᵖ
Chaque individu \(i\) est représenté par le vecteur de ses \(p\) mesures : un point dans \(\mathbb{R}^p\). L'ensemble des \(n\) points forme un nuage dont la forme traduit les relations entre individus. Deux individus proches ont des profils similaires pour toutes les variables.
Interprétation : deux individus sont proches s'ils ont des valeurs similaires pour toutes les \(p\) variables.
2. Centrage — translation au centre de gravité
Le principe d'ajustement consiste à maximiser la variance des projections, c'est-à-dire la dispersion autour de la moyenne. On translate donc l'origine au centre de gravité \(G\) du nuage. Le tableau centré \(X\) a pour terme général :
rᵢⱼ : donnée brute · r̄ⱼ : moyenne de la variable \(j\) · xᵢⱼ : donnée centréePropriété clé : \(\sum_{i,i'} d^2(i,i') = 2n \sum_i d^2(i,G)\) — maximiser les distances entre paires est équivalent à maximiser les distances au centre de gravité.
3. ACP normée — réduction des variables
Si les variables ont des unités différentes (cm, kg, €...), celles à grande variance dominent l'analyse. Pour donner un rôle identique à chaque variable, on réduit en divisant par l'écart-type :
Conséquence : la matrice \(C = X^TX\) devient la matrice de corrélation.
4. Visualisation — centrage et nuage
Données brutes R
Données centrées X (origine = G)
La croix jaune = centre de gravité G. Après centrage, G est à l'origine. La forme du nuage ne change pas.
5. La matrice à diagonaliser C = X'X
On cherche le vecteur unitaire \(u\) qui maximise la somme des carrés des projections \(\sum_i (OH_i)^2 = u^T X^T X u\) sous la contrainte \(u^T u = 1\). D'après le théorème des extrema de formes quadratiques sous contrainte quadratique, la solution est le vecteur propre de \(X^TX\) associé à la plus grande valeur propre.
C : matrice de taille \(p \times p\) (matrice de corrélation en ACP normée)uα : \(\alpha\)-ième vecteur propre = direction du \(\alpha\)-ième axe principalλα : \(\alpha\)-ième valeur propre = variance des projections sur cet axeLes vecteurs propres \(u_1, u_2, \dots, u_p\) sont orthogonaux et de norme 1.
Le terme général \(C_{jj'}\) de la matrice \(C = X^TX\) est :
6. Exemple numérique complet
Appliquons la démarche sur un petit dataset — 4 individus, 2 variables.
| Individu | X₁ | X₂ |
|---|---|---|
| i₁ | 2 | 4 |
| i₂ | 4 | 6 |
| i₃ | 6 | 8 |
| i₄ | 8 | 2 |
\(\bar{x}_1 = \frac{2+4+6+8}{4} = 5\) · \(\bar{x}_2 = \frac{4+6+8+2}{4} = 5\)
\(s_1^2 = \frac{(2-5)^2+(4-5)^2+(6-5)^2+(8-5)^2}{4} = \frac{9+1+1+9}{4} = 5\) → \(s_1 = \sqrt{5}\)
\(s_2^2 = \frac{(4-5)^2+(6-5)^2+(8-5)^2+(2-5)^2}{4} = 5\) → \(s_2 = \sqrt{5}\)
| Individu | x₁ = (r₁−5)/(√5·2) | x₂ = (r₂−5)/(√5·2) |
|---|---|---|
| i₁ | −3/(2√5) | −1/(2√5) |
| i₂ | −1/(2√5) | 1/(2√5) |
| i₃ | 1/(2√5) | 3/(2√5) |
| i₄ | 3/(2√5) | −3/(2√5) |
\(C_{11} = \sum x_{i1}^2 = \frac{9+1+1+9}{20} = 1\) · \(C_{22} = 1\) (vérification : diagonale = 1 ✓)
\(C_{12} = \sum x_{i1}x_{i2} = \frac{(-3)(-1)+(- 1)(1)+(1)(3)+(3)(-3)}{20} = \frac{3-1+3-9}{20} = \frac{-4}{20} = -0.2\)
\(C = \begin{pmatrix}1 & -0.2 \\ -0.2 & 1\end{pmatrix}\)
\(\det(C - \lambda I) = (1-\lambda)^2 - 0.04 = 0\)
\((1-\lambda) = \pm 0.2\) → \(\lambda_1 = 1.2\) et \(\lambda_2 = 0.8\)
Vérification : \(\lambda_1 + \lambda_2 = 2 = p\) ✓ (ACP normée, somme = p)
Le premier axe capte \(\frac{1.2}{2} = 60\%\) de l'inertie totale. Le second \(40\%\). Les variables sont légèrement anticorrélées (\(r = -0.2\)).
7. Calculateur interactif — matrice de corrélation
Modifie les données brutes (5 individus × 2 variables). La matrice de corrélation se calcule en temps réel.
Données brutes R
| X₁ | X₂ |
|---|
📌 À retenir / À l'examen
- Chaque individu = un point dans \(\mathbb{R}^p\). La distance entre deux individus = distance euclidienne.
- Centrage : \(x_{ij} = r_{ij} - \bar{r}_j\) — translate l'origine au centre de gravité G
- ACP normée : \(x_{ij} = (r_{ij}-\bar{r}_j)/(s_j\sqrt{n})\) — toutes les variables ont la même norme
- Matrice à diagonaliser : \(C = X^TX\) — en ACP normée, c'est la matrice de corrélation
- Terme général : \(C_{jj'} = \text{cor}(j,j')\) — la diagonale vaut toujours 1
- Les axes principaux \(u_\alpha\) = vecteurs propres de \(C\), la variance sur l'axe = valeur propre \(\lambda_\alpha\)
- Propriété fondamentale : \(\sum_{i,i'} d^2(i,i') = 2n \sum_i d^2(i,G)\)