Voir la source de Analyse en composantes principales

{{Voir homonymes|ACP|PCA|KLT (homonymie)}}

L''''analyse en composantes principales''' ('''ACP''' ou '''PCA''' en anglais pour ''{{lang|en|principal component analysis}}''), ou selon le domaine d'application la '''transformation de Karhunen–Loève''' ('''KLT''')<ref>http://tcts.fpms.ac.be/cours/1005-07-08/codage/codage/xcodim2.pdf</ref>,  est une méthode de la famille de l'[[analyse des données]] et plus généralement de la  [[statistique multivariée]], qui consiste à transformer des variables liées entre elles (dites « corrélées » en statistique) en nouvelles variables décorrélées les unes des autres. Ces nouvelles variables sont nommées « composantes principales », ou axes principaux. Elle permet au praticien de réduire le nombre de variables et de rendre l'information moins redondante.

Il s'agit d'une approche à la fois géométrique<ref>Une vidéo d’introduction à l’ACP fondée sur la géométrie est accessible [http://math.agrocampus-ouest.fr/infoglueDeliverLive/enseignement/support2cours/videos ici.]</ref> (les variables étant représentées dans un nouvel espace, selon des directions d'inertie maximale) et statistique (la recherche portant sur des axes indépendants expliquant au mieux la variabilité — la [[Variance (statistiques et probabilités)|variance]] — des données). Lorsqu'on veut ''compresser'' un ensemble de <math>N</math> variables aléatoires, les <math>n</math> ''premiers axes'' de l'analyse en composantes principales sont un ''meilleur choix'', du point de vue de l'inertie ou de la variance.

L'outil mathématique est appliqué dans d'autres domaines que les statistiques et est parfois appelé '''décomposition orthogonale aux valeurs propres''' ou '''POD''' ({{lang-en|''proper orthogonal decomposition''}}).

== Histoire et applications ==
[[Fichier:Karl Pearson line of best fit diagramm from philosophical magazine 1901 2 559-572.jpg|thumb|right|200 px|Extrait de l'article de Pearson de 1901 : la recherche de la « droite du meilleur ajustement ».]]

L'ACP prend sa source dans un article de [[Karl Pearson]] publié en 1901<ref>{{article |langue=en| auteur = Pearson, K. | année = 1901 | titre = On Lines and Planes of Closest Fit to Systems of Points in Space | journal = Philosophical Magazine | volume = 2 | numéro = 6 | pages = 559–572 | url = http://stat.smmu.edu.cn/history/pearson1901.pdf |format=PDF}}</ref>. Le père du [[test du χ²]] y prolonge ses travaux dans le domaine de la [[Régression (statistiques)|régression]] et des [[Corrélation (statistiques)|corrélations]] entre plusieurs variables. Pearson utilise ces corrélations non plus pour expliquer une variable à partir des autres (comme en régression), mais pour décrire et résumer l'information contenue dans ces variables.

Encore connue sous le nom de transformée de Karhunen-Loève ou de transformée de Hotelling, l'ACP a été de nouveau développée et formalisée dans les années 1930 par [[Harold Hotelling]]<ref>{{en}} « {{Langue|en|Analysis of a Complex of Statistical Variables with Principal Components}} », 1933, ''{{Langue|en|Journal of Educational Psychology}}''.</ref>. La puissance mathématique de l'économiste et statisticien américain le conduira aussi à développer l'[[analyse canonique des corrélations]], généralisation des analyses factorielles dont fait partie l'ACP.

Les champs d'application sont aujourd'hui multiples, allant de la [[biologie]] à la [[Recherche en économie|recherche économique et sociale]], et plus récemment le [[traitement d'images]]. L'ACP est majoritairement utilisée pour:
* décrire et visualiser des données ;
* les décorréler ; la nouvelle base est constituée d'axes qui ne sont pas corrélés entre eux ;
* les [[débruitage|débruiter]], en considérant que les axes que l'on décide d'''oublier'' sont des axes ''bruités''.

La méthode de décomposition orthogonales aux valeurs propres est également utilisée :
* pour réduire la taille de modèle d'[[éléments finis]]<ref>{{lien web
 | url = https://scilab.io/introduction-to-model-reduction/
 | titre = introduction to model reduction
 | site = [[Scilab]]
 | date = 2017-02-24 | consulté le = 2018-03-15
}}</ref>.

== Exemples introductifs ==
[[Image:PCA fish.png|thumb|right|Les deux axes d'une ACP sur la photo d'un poisson.]]

'''Premier exemple'''

Dans le cas d'une image, comme dans la figure ci-contre, les [[pixel]]s sont représentés dans un plan et considérés comme une [[variable aléatoire]] à deux dimensions. L'ACP va déterminer les deux axes qui expliquent le mieux la dispersion de l'objet, interprété comme un nuage de points. Elle va aussi les ordonner par inertie expliquée, le second axe étant perpendiculaire au premier.

'''Second exemple'''

Dans une école imaginaire, on n'enseigne que deux matières sur lesquelles les élèves sont notés : le français et les mathématiques. En appliquant l'ACP au tableau de notes, on dégagera probablement en premier axe des valeurs par élève très proches de leur moyenne générale dans les deux matières. C'est cet axe qui résumera au mieux la variabilité des résultats selon les élèves. Mais un professeur voulant pousser l'analyse des résultats, s'intéressa aussi au second axe, qui ordonne les élèves selon l'ampleur de leurs écarts entre les deux notes, et indépendamment du premier axe.

On comprend l'intérêt de la méthode d'ACP quand on étend l'analyse à 10 matières enseignées : la méthode va calculer pour chaque élève 10 nouvelles valeurs, selon 10 axes, chacun étant indépendant des autres. Les derniers axes apporteront très peu d'information sur le plan statistique : ils mettront probablement en évidence quelques élèves au profil singulier. Selon son point de vue d'analyse, le professeur, dans sa pratique quotidienne, veillera donc plus particulièrement à ces élèves qui auront été mis en évidence par les derniers axes de la méthode ACP, et corrigera peut-être une erreur qui se serait glissée dans son tableau de notes, mais à l'inverse, il ne prendra pas en compte ces derniers axes s'il mène une réflexion globale s'intéressant aux caractéristiques pédagogiques majeures, ou autrement dit, principales.
Si on prend pour exemple une classe de {{1re}} S, on a de fortes chances pour avoir comme axe principal un regroupement des matières scientifiques, et comme second axe les matières littéraires. Ces deux variables expliquent les notes obtenues par les élèves de la classe.

La puissance de l'ACP est qu'elle sait aussi prendre en compte des données de nature hétérogène : par exemple un tableau des différents pays du monde avec le PNB par habitant, le taux d'alphabétisation, le taux d'équipement en téléphones portables, le prix moyen du hamburger, etc. Elle permet d'avoir une intuition rapide des effets conjoints entre ces variables.

== Échantillon ==

On applique généralement une ACP sur un ensemble de ''N'' variables aléatoires ''X''<sub>1</sub>, …, ''X<sub>N</sub>'' connues à partir d'un échantillon de <math>K</math> réalisations conjointes de ces variables.

Cet échantillon de ces ''N'' variables aléatoires peut être structuré dans une matrice ''M'', à ''K'' lignes et ''N'' colonnes.

: <math>M=\begin{bmatrix} X_{1,1} & \cdots & X_{1,N} \\ \vdots & \ddots & \vdots \\ X_{K,1} & \cdots & X_{K,N}\end{bmatrix}</math>

Chaque variable aléatoire ''X<sub>n</sub>'', dont ''X''<sub>''1'', n</sub>, …, ''X''<sub>''K'', ''n''</sub> sont des réalisations indépendantes, a une moyenne <math>\bar X_n</math> et un écart type σ<sub>''Xn''</sub>.

=== Poids ===
Si les réalisations (les éléments de la matrice ''M'') sont à probabilités égales alors chaque réalisation (un élément <math>X_{i,j}</math> de la matrice) a la même importance <math>1/K</math> dans le calcul des caractéristiques de l'échantillon.  On peut aussi appliquer un poids <math>p_{i}</math> différent à chaque réalisation conjointe des variables (cas des échantillons redressés, des données regroupées, ...).
Ces poids, qui sont des nombres positifs de somme 1 sont représentés par une matrice diagonale ''D'' de taille ''K'':
: <math>D=\begin{bmatrix} p_{1} & & & 0 \\ & p_{2} & & \\ & & \ddots & \\ 0 & & & p_{K}\end{bmatrix}</math>

Dans le cas le plus courant de poids égaux, <math>D = {1 \over K} I</math> où <math>I</math> est la [[matrice identité]].

=== Transformations de l'échantillon ===

Le vecteur <math>(\bar X_1, \cdots, \bar X_N)</math> est le [[centre de gravité]] du nuage de points ; on le note souvent ''g''.
On a '''<math>g = M^{T}D\tilde{1}</math>''' où <math>\tilde{1}</math> désigne le vecteur de <math>{\mathbb{R}}^K</math> dont toutes les coordonnées sont égales à 1.

La matrice ''M'' est généralement centrée sur le centre de gravité :
: <math>\bar M=\begin{bmatrix} X_{1,1}-\bar X_1 & \cdots & X_{1,N}-\bar X_N \\ \vdots & \ddots & \vdots \\ X_{K,1}-\bar X_1 & \cdots & X_{K,N}-\bar X_N\end{bmatrix} = M - \tilde{1}g^{T}</math>.

Elle peut être aussi '''réduite''' :
: <math>\tilde M=\begin{bmatrix} {X_{1,1}-\bar X_1\over \sigma(X_1)} & \cdots & {X_{1,N}-\bar X_N\over \sigma(X_N)} \\ \vdots & \ddots & \vdots \\ {X_{K,1}-\bar X_1\over \sigma(X_1)} & \cdots & {X_{K,N}-\bar X_N\over \sigma(X_N)}\end{bmatrix}</math>.

Le choix de réduire ou non le nuage de points (i.e. les ''K'' réalisations de la variable aléatoire (''X''<sub>1</sub>, …, ''X<sub>N</sub>'')) est un choix de modèle :
* si l'on ne réduit pas le nuage : une variable à forte variance va « tirer » tout l'effet de l'ACP à elle ;
* si l'on réduit le nuage : une variable qui n'est qu'un bruit va se retrouver avec une variance apparente égale à une variable informative.
Si les variables aléatoires sont dans des unités différentes, la réduction est obligatoire.

=== Calcul des covariances et des corrélations ===

Dans le cas de poids uniformes (<math>D = {1 \over K} I</math>), une fois la matrice <math>M</math> transformée en <math>\bar M</math> ou <math>\tilde M</math>, il suffit de la multiplier par sa [[matrice transposée|transposée]] pour obtenir:

* la [[matrice de variance-covariance]] des ''X''<sub>1</sub>, …, ''X<sub>N</sub>'' si ''M'' n'est pas réduite : <math>{\rm Covariances} = 1/K \cdot {}\bar M \cdot \bar M^{T}</math>;
* la [[matrice de corrélation]] des ''X''<sub>1</sub>, …, ''X<sub>N</sub>'' si ''M'' est réduite : <math>{\rm Corr\acute elations} = 1/K \cdot {}\tilde M \cdot \tilde M^{T}</math>.

Ces deux matrices sont [[matrice carrée|carrées]] (de taille ''N''), [[matrice symétrique|symétriques]], et [[nombre réel|réelles]]. Elles sont donc [[matrice diagonalisable|diagonalisables]] dans une base orthonormée en vertu du [[Théorème spectral#Formalisation algébrique|théorème spectral]].

De façon plus générale, dans le cas de ''poids non uniformes'', la matrice de variance-covariance s'écrit <math>V = MDM^{T} - gg^{T} = \bar M \cdot D \cdot \bar M^{T}</math>.

De plus, si l'on note <math>D_{1/s}</math> la matrice diagonale des inverses des écarts-types:
: <math>D_{1/s} = \begin{bmatrix} 1/s_{1} & & 0 \\ & \ddots & \\ 0 & & 1/s_{N}\end{bmatrix}= \begin{bmatrix} 1/\sigma(X_1) & & 0 \\ & \ddots & \\ 0 & & 1/\sigma(X_N) \end{bmatrix} </math>
alors on a:
: <math>\tilde M = \bar M \cdot D_{1/s}</math>.
La matrice des coefficients de corrélation linéaire entre les ''N'' variables prises deux à deux, notée '''''R''''', s'écrit:
:<math>R = \tilde M \cdot D \cdot \tilde M^{T}=D_{1/s}VD_{1/s}</math>.

== Critère d'inertie ==

Dans la suite de cet article, nous considèrerons que le nuage est transformé (centré et réduit si besoin est). Chaque ''X<sub>n</sub>'' est donc remplacé par <math>X_n-\bar X_n</math> ou <math>(X_n-\bar X_n)/\sigma(X_n)</math>. Nous utiliserons donc la matrice <math>M</math> pour noter <math>\bar M</math> ou <math>\tilde M</math> suivant le cas.

Le principe de l'ACP est de trouver un axe ''u'', issu d'une combinaison linéaire des ''X<sub>n</sub>'', tel que la variance du nuage autour de cet axe soit maximale.

Pour bien comprendre, imaginons que la variance de ''u'' soit égale à la variance du nuage; on aurait alors trouvé une combinaison des ''X<sub>n</sub>'' qui contient toute la ''diversité'' du nuage original (en tout cas toute la part de sa diversité captée par la variance).

Un critère couramment utilisé est la [[Variance (statistiques et probabilités)|variance]] de l'échantillon (on veut maximiser la variance expliquée par le vecteur ''u''). Pour les physiciens, cela a plutôt le sens de maximiser l'inertie expliquée par ''u'' (c'est-à-dire minimiser l'inertie du nuage autour de ''u'').

=== Projection ===

Finalement, nous cherchons le vecteur ''u'' tel que la [[Projection orthogonale|projection]] du nuage sur <math>u</math> ait une variance maximale.
La projection de l'échantillon des ''X'' sur ''u'' s'écrit :
: <math>\pi_u(M) = M \cdot u</math>

la variance empirique de <math>\pi_u(M)</math> vaut donc :

: <math>\pi_u(M)^{T} \cdot 1/K \cdot \pi_u(M) = u^{T}\cdot \underbrace{M^{T}\cdot  1/K \cdot M}_C \cdot u</math>
où ''C'' est la matrice de covariance.

Comme nous avons vu plus haut que ''C'' est [[diagonalisable]] dans une base orthonormée, notons ''P'' le [[Matrice de passage|changement de base]] associé et <math>\Delta = \textrm{Diag}(\lambda_1, \ldots, \lambda_N)</math> la matrice diagonale formée de son spectre :

: <math>\pi_u(M)^{T} \cdot 1/K \cdot \pi_u(M) = u^{T} P^{T} \Delta P u = (Pu)^{T} \Delta \underbrace{(Pu)}_v</math>

Les valeurs <math>(\lambda_1, \ldots, \lambda_N)</math> de la diagonale de <math>\Delta</math> sont rangées en ordre décroissant. Le vecteur unitaire ''u'' qui maximise <math>v^{T} \Delta v</math> est un [[vecteur propre]] de ''C'' associé à la [[valeur propre]] <math>\lambda_1</math> ; on a alors :

: <math> v^{T} \cdot \Delta \cdot v = \lambda_1 </math>
La valeur propre <math>\lambda_1</math> est la variance empirique sur le premier axe de l'ACP.

Il est aussi possible de démontrer ce résultat en maximisant la variance empirique des données projetées sur ''u'' sous la contrainte que ''u'' soit de norme 1 (par un [[multiplicateur de Lagrange]] <math>\alpha</math>) :

: <math> L(u,\alpha) = u^{T} \cdot C \cdot u - \alpha (u^{T} u -1) </math>

{{démonstration|déroulante=oui|

La première condition nécessaire d'extremum local sous contrainte est l'ensemble des relations : pour tout <math>i \in \{1, ..., N\}</math>, 

:<math> \frac{\partial L}{\partial u_i} = 2 \frac{\partial u^T}{\partial u_i} C u - 2 \alpha u_i = 0 </math>

<math> \frac{\partial u^T}{\partial u_i} </math> est simplement la transposée du ''i''-ème vecteur de la base canonique de '''R'''<sup>''N''</sup> (toutes les coordonnées nulles sauf la ''i''-ème égale à 1) donc <math> \frac{\partial u^T}{\partial u_i} C </math> est la ''i''-ème ligne de ''C'', de sorte que l'ensemble des ''N'' équations ci-dessus est équivalent à <math> C u=\alpha u </math> : '''α est une valeur propre de ''C'''''.

De plus, <math> u^T C u = \alpha ||u||^2=\alpha </math> (d'après la seconde condition de Lagrange sur l'annulation des dérivées partielles <math> \frac{\partial L}{\partial \alpha} = - (||u||^2 - 1) = 0 </math>) qui est donc '''maximum lorsque <math> \alpha = \lambda_1 </math>'''.

}}

On continue la recherche du deuxième axe de projection ''w'' sur le même principe en imposant qu'il soit orthogonal à ''u''.

=== Diagonalisation ===

La [[diagonalisation]] de la matrice de corrélation (ou de covariance si on se place dans un modèle non réduit), nous a permis d'écrire que le vecteur qui explique le plus d'inertie du nuage est le premier [[Valeur propre (synthèse)|vecteur propre]]. De même le deuxième vecteur qui explique la plus grande part de l'inertie restante est le deuxième vecteur propre, etc.

Nous avons vu en outre que la variance expliquée par le ''k''-ième vecteur propre vaut λ<sub>''k''</sub>.

Finalement, la question de l'ACP se ramène à un problème de diagonalisation de la matrice de corrélation.

=== Optimisation numérique ===

Numériquement, la matrice ''M'' étant rectangulaire, il peut être plus économique de la [[Décomposition en valeurs singulières|décomposer en valeurs singulières]], puis de recombiner la décomposition obtenue, plutôt que de diagonaliser ''M{{'}} M''.

== ACP et variables qualitatives ==

En ACP, il est fréquent que l’on veuille introduire des variables qualitatives supplémentaires. Par exemple, on a mesuré de nombreuses variables quantitatives sur des plantes. Pour ces plantes, on dispose aussi de variables qualitatives, par exemple l’espèce à laquelle appartient la plante. On soumet ces données à une ACP des variables quantitatives. Lors de l’analyse des résultats, il est naturel de chercher à relier les composantes principales à la variable qualitative ''espèce''.
Pour cela on produit les résultats suivants :
* identification, sur les plans factoriels, des différentes espèces en les représentant par exemple par des couleurs différentes ;
* représentation, sur les plans factoriels, des centres de gravité des plantes appartenant à une même espèce ;
* indication, pour chaque centre de gravité et pour chaque axe, d'une probabilité critique pour juger de la significativité de l’écart entre un centre de gravité et l’origine.

Tous ces résultats constituent ce que l’on appelle ''introduire une variable qualitative en supplémentaire''. Cette procédure est détaillée dans {{Harvsp|texte= Escofier&Pagès 2008|id=Escofier&Pagès}}, {{Harvsp|texte=Husson, Lê & Pagès 2009|id=Husson et al }} et {{Harvsp|texte=Pagès 2013|id=Pagès 2013}}.

Peu de logiciels offrent cette possibilité de façon « automatique ». C’est le cas de SPAD<ref>http://www.coheris.com/produits/analytics/logiciel-data-mining/</ref> qui historiquement, à la suite des travaux de Ludovic Lebart, a été le premier logiciel à le proposer, et du package R FactoMineR<ref>http://factominer.free.fr/</ref>.

== Résultats théoriques ==

Si les sections précédentes ont travaillé sur un échantillon issu de la loi conjointe suivie par ''X''<sub>1</sub>, …, ''X<sub>N</sub>'', que dire de la validité de nos conclusions sur n'importe quel autre échantillon issu de la même loi ?

Plusieurs résultats théoriques permettent de répondre au moins partiellement à cette question, essentiellement en se positionnant par rapport à une [[distribution gaussienne]] comme référence.

== Méthodes voisines et extensions : la famille factorielle ==
L’analyse en composantes principales est la plus connue des méthodes factorielles ; d’autres méthodes factorielles existent pour analyser d’autres types de tableau. À chaque fois, le principe général est le même.

* On considère deux nuages de points, l’un associé aux lignes du tableau analysé et l’autre aux colonnes de ce tableau.
* Ces deux nuages sont liés par des relations de dualité : ils ont la même inertie totale ;
* Chacun de ces nuages est projeté sur ses directions d’inertie maximum.
* D’un nuage à l’autre, les directions d’inertie de même rang sont liées par des relations de dualité (ou de transition) : elles ont la même inertie et les coordonnées des projections sur l’une se déduisent des coordonnées des projections sur l’autre.

=== [[Analyse factorielle des correspondances]] (AFC) ===
Elle s’applique à des tableaux de contingence c'est-à-dire des tableaux croisant deux variables qualitatives. Ce type de tableau est très différent de celui analysé par ACP : en particulier, les lignes et les colonnes jouent des rôles symétriques alors que la distinction entre lignes et colonnes (c'est-à-dire entre individus et variables) est majeure en ACP.

=== [[Analyse des correspondances multiples]] (ACM) ===
Elle s’applique à des tableaux dans lesquels un ensemble d’individus est décrit par un ensemble de variables qualitatives. Ce type de tableau est donc voisin de celui analysé en ACP, les variables quantitatives étant remplacées par des variables qualitatives. L’ACM est souvent vue comme un cas particulier de l’AFC mais ce point de vue est très réducteur. L’ACM possède suffisamment de propriétés spécifiques pour être considérée comme une méthode à part entière.

On peut aussi présenter l’ACM à partir de l’ACP comme cela est fait dans {{Harvsp|texte=Pagès 2013|id=Pagès 2013}}. L’intérêt est de relier entre eux les ressorts de l’ACP et ceux de l’ACM ce qui ouvre la voie au traitement simultané des deux types de variables (cf. AFDM et AFM ci-après)

=== [[Analyse factorielle de données mixtes]] (AFDM) ===
Les données sont constituées par un ensemble d’individus pour lesquels on dispose de plusieurs variables, comme en ACP ou en ACM. Mais, ici, les variables sont aussi bien quantitatives que qualitatives.
L’analyse factorielle de données mixtes traite simultanément les deux types de variables en leur faisant jouer un rôle actif. L’AFDM est décrite dans {{Harvsp|texte=Pagès 2013|id=Pagès 2013}} et {{Harvsp|texte= Escofier&Pagès 2008|id=Escofier&Pagès}}.

=== [[Analyse factorielle multiple]] (AFM) ===
Les données sont, ici encore, constituées par un ensemble d’individus pour lesquels on dispose de plusieurs variables. Mais cette fois, outre qu’elles peuvent être quantitatives et/ou qualitatives, les variables sont structurées en groupes. Ce peut être, par exemple, les différents thèmes d’un questionnaire. L’AFM prend en compte cette structure en groupes dans l’analyse de ces données. L’AFM est décrite en détail dans {{Harvsp|texte=Pagès 2013|id=Pagès 2013}} et {{Harvsp|texte= Escofier&Pagès 2008|id=Escofier&Pagès}}.

== Applications ==

=== Compression ===
L'Analyse en Composantes Principales est usuellement utilisée comme outil de [[Compression statistique|compression]] linéaire. Le principe est alors de ne retenir que les ''n'' premiers vecteurs propres issus de la diagonalisation de la matrice de corrélation (ou covariance), lorsque l'inertie du nuage projeté sur ces ''n'' vecteurs représente ''q<sub>n</sub>'' pourcents de l'inertie du nuage original, on dit qu'on a un taux de compression de 1 - ''q<sub>n</sub>'' pourcents, ou que l'on a compressé à ''q<sub>n</sub>'' pourcents. Un taux de compression usuel est de 20 %.

Les autres méthodes de compressions statistiques habituelles sont:
* l'[[analyse en composantes indépendantes]] ;
* les [[Carte auto adaptative|cartes auto-adaptatives]] (SOM, ''{{lang|en|self organizing maps}}'' en anglais) ; appelées aussi cartes de [[Teuvo Kohonen|Kohonen]] ;
* l'[[analyse en composantes curvilignes]] ;
* la [[compression par ondelettes]].

Il est possible d'utiliser le résultat d'une ACP pour construire une [[classification statistique]] des variables aléatoires ''X''<sub>1</sub>, …, ''X<sub>N</sub>'', en utilisant la [[Distance (mathématiques)|distance]] suivante (''C<sub>n, n{{'}}</sub>'' est la corrélation entre ''X<sub>n</sub>'' et ''X<sub>n{{'}} </sub>''):

: <math>d(X_n,X_{n'})=\sqrt{2\,(1-C_{n,n'})}</math>

=== Analyse de séries dynamiques d'images ===
L'ACP, désignée en général dans le milieu du traitement du signal et de l'analyse d'images plutôt sous son nom de Transformée de Karhunen-Loève (TKL) est utilisée pour analyser les séries dynamiques d'images<ref>Évaluation de la perfusion et de la fonction contractile du myocarde à l’aide de l’analyse de Karhunen-Loève en tomographie d’émission monophotonique myocardique synchronisée à l’ECG par P. Berthout, R. Sabbah, L. Comas, J. Verdenet, O. Blagosklonov, J.-C. Cardot et M. Baud dans ''Médecine nucléaire'' Volume 31, Volume 12, décembre 2007, Pages 638-646.</ref>, c'est-à-dire une succession d'images représentant la cartographie d'une grandeur physique, comme les [[scintigraphie]]s dynamiques en [[médecine nucléaire]], qui permettent d'observer par gamma-caméra le fonctionnement d'organes comme le cœur ou les  [[scintigraphie rénale|reins]].

Dans une série de P images, chaque [[pixel]] est considéré comme un point d'un [[espace affine]] de dimension P dont les coordonnées sont la valeur du pixel pour chacune des P images au cours du temps. Le nuage ainsi formé par tous les points de l'image peut être analysé par l'ACP, (il forme un hyper-ellipsoïde à P dimensions) ce qui permet de déterminer ses axes principaux.

En exprimant tous les points dans le repère orthogonal à P dimensions des axes de l'ACP, on passe ainsi de la série temporelle d'origine (les pixels représentent la valeur en fonction du temps) à une nouvelle série (également de P images) dans l'espace de Karhunen-Loève : c'est la Transformée de Karhunen-Loève, qui est une opération réversible : on parle de « TKL » et de « TKL inverse » ou « TLK{{exp|-1}} ».

La compression est possible car l'information est contenue presque entièrement sur les premiers axes de l'ACP. Mais la notion de « compression » sous-entend que les autres images correspondant aux autres axes sont volontairement ignorées. La TKL étant réversible, la suppression arbitraire des axes les moins énergétiques constitue alors un filtrage permettant de réduire le bruit temporel de la série d'images.

Concrètement, l'application de TKL + suppression des axes les moins significatifs + TKL{{exp|-1}} permet de supprimer le fourmillement apparent (bruit temporel) d'une série animée d'images.

En imagerie médicale fonctionnelle, on améliore ainsi la qualité visuelle de la visualisation scintigraphique du cycle cardiaque moyen.

Par ailleurs, l'analyse de l'importance respective des valeurs propres de l'ACP permet d'approcher le nombre de fonctionnements physiologiques différents. On a ainsi pu montrer que le cœur sain peut être entièrement représenté avec 2 images (2 axes de l'ACP contiennent toute l'information utile), alors que pour certaines pathologies l'information utile s'étale sur 3 images<ref>Baud, Cardot, Verdenet et al, Service de médecine nucléaire, Hôpital Jean-Minjoz, boulevard Fleming, 25030 Besançon cedex, France (nombreuses publications sur plus de 30 ans)</ref>.

=== Analyse d'images multi-spectrales ===
Comme pour l'application précédente, la longueur d'onde remplaçant juste le temps, la TKL a été proposée à plusieurs reprises pour extraire l'information utile d'une série d'images monochromes représentant les intensités pour des longueurs d'onde différentes. De telles images peuvent être issues de microscopie optique classique, confocale ou SNOM ([[Microscope optique en champ proche]])<ref>Analysis of optical near-field images by Karhunen—Loève transformation Daniel Charraut, Daniel Courjon, Claudine Bainier, and Laurent Moulinier, ''Applied Optics'', Vol. 35, Issue 20, {{p.|3853-3861}} (1996)</ref>.

=== Évolution de la topographie ===
De la même manière, la TKL permet de mettre en évidence des cinétiques différentes lors de l'analyse topographique dynamique, c'est-à-dire l'analyse de l'évolution du relief au cours du temps. Elle permet alors de déceler des phénomènes invisibles par simple observation visuelle, mais se distinguant par une cinétique légèrement différente (par exemple pollution d'une surface rugueuse par un dépôt)<ref>(en) [http://www.micronora.com/micronora_infos/hs_en.pdf Jean-Yves Catherin, Measure in 2D, visualise in 3D and understand in 4D dans ''Micronora Informations'' juin 2008, page 3]</ref>.

== Logiciels ==
=== R ===
Le package R FactoMineR inclut, en particulier, outre l'ACP, toutes les extensions décrites plus haut : AFC, ACM, AFDM et AFM). Ce logiciel est relié au livre {{Harvsp|texte=Husson, Lê & Pagès 2009|id=Husson et al}}.

Citons également les packages ade4<ref>{{Lien web|nom1=Penel,|prénom1=Daniel Chessel, Stephane Dray, Anne Beatrice Dufour , Jean Thioulouse,  Simon|titre=Accueil ade4|url=http://pbil.univ-lyon1.fr/ade4/|site=pbil.univ-lyon1.fr|consulté le=2017-11-25}}</ref> et adegraphics qui permettent d'analyser et visualiser (ACP inter, intra classe, etc.) des données issues de problématiques biologiques.

==== Exemple sous R ====
L'objectif est de réaliser une ACP sur des données physico-chimiques récoltées par 30 stations le long du Doubs<ref>{{Lien web|titre=R: Pair of Ecological Tables|url=https://pbil.univ-lyon1.fr/ade4/ade4-html/doubs.html|site=pbil.univ-lyon1.fr|consulté le=2017-11-25}}</ref>. 

Il y a 7 variables physico-chimiques d’intérêt, mesurées au niveau des 30 stations : le pH, la dureté de l’eau, le phosphate, le nitrate, l'ammoniaque, l'oxygène et la demande biologique en oxygène. Quelles données sont corrélées ? Peut-on représenter la pollution le long de la rivière ?

Les données étant exprimées dans des unités différentes, nous avons intérêt à réaliser une ACP normée (centrée et réduite) sur les 7 variables physico-chimiques.
<syntaxhighlight lang="r">
library(ade4)
library(adegraphics)
data(doubs)
ACP <- dudi.pca(doubs$env[,5:11], center = T, scale = T, scannf=FALSE, nf = 2)
#On récupère les 5 variables physico-chimiques, on centre les données, on les réduit, on défini le nombre d'axes pour l'ACP à 2)
ACP$eig #On récupère les valeurs propres
s.corcircle(ACP$co, xax = 1, yax = 2) #On affiche le cercle des corrélations sur le premier plan factoriel de l'ACP
s.label(ACP$li, xax = 1, yax = 2) #On affiche les stations dans le premier plan factoriel de l'ACP
s.value(doubs$xy,ACP$li[,1], symbol = "circle", col = c("blue","brown"))#Les stations sont symbolisées par des cercles pour la valeur du 1er axe de l'ACP (la pollution). Les stations peu polluées sont en bleu, les stations plus polluées sont en marron
</syntaxhighlight>L'interprétation de l'ACP sur le premier axe du <math>s.corcircle</math><ref>{{Lien web|titre=R: Plot of the factorial maps of a correlation circle|url=https://pbil.univ-lyon1.fr/ade4/ade4-html/s.corcircle.html|site=pbil.univ-lyon1.fr|consulté le=2017-11-25}}</ref> montre une corrélation négative entre le taux d'oxygène et la présence de polluants. Avec <math>s.value</math><ref>{{Lien web|titre=R: Representation of a value in a graph|url=https://pbil.univ-lyon1.fr/ade4/ade4-html/s.value.html|site=pbil.univ-lyon1.fr|consulté le=2017-11-25}}</ref>, on peut représenter les stations avec leurs valeurs sur cet axe de pollution, symbolisées par des cercles de différentes couleurs.
<gallery caption="Figures pour l'exemple d'ACP sous R" widths="300px" heights="300px">
Fichier:Stations le long du doubs.svg|Stations le long du doubs.
Fichier:Cercle des corélations d'une ACP normée.svg|Le cercle de corrélation (corcircle) tracé pour le premier plan factoriel de l'ACP sur les données du doubs. Le premier axe montre une corrélation négative entre l'oxygène et 3 polluants connus : nitrates (nit), phosphates (pho), ammoniaque (amm). On peut dire que cet axe est caractérisé par la pollution. Les eaux riches en oxygène n'ont pas de polluants et inversement. Le deuxième axe est un axe caractérisé essentiellement par le pH.
Fichier:Stations dans le 1er plan factoriel de l'ACP doubs.svg|Après interprétation du cercle des corrélations, on peut dire que les stations 25, 23, 24 sont les plus polluées. La station 15 a un pH élevé.
Fichier:Résultats du premier axe de l'ACP sur les stations.svg|Chaque station est remplacée par un cercle prenant la valeur du résultat de l'ACP sur le premier axe. Le cercle des corrélations a indiqué qu'il s'agit d'un axe de pollution. On remarque que la deuxième partie du fleuve est plus polluée que la première partie. La station 25 a la valeur de pollution maximale. Cela est cohérent avec les observations du 1er plan factoriel de l'ACP.
</gallery>

== Notes ==
<references />

== Voir aussi ==
* [[valeur propre|Valeurs propres]]
* [[Compression statistique]]
* Équilibre biais / variance
* [[Analyse de la variance]]
* [[Partitionnement de données]]
* [[Exploration de données]]
* [[Iconographie des corrélations]]
* [[Michel Loève]]
* [[Kari Karhunen]]
* {{Lien|Théorème de Karhunen-Loève|trad=Karhunen–Loève theorem}}
* [[Analyse discriminante linéaire]]

== Références ==

* [[Jean-Paul Benzécri]] ; ''Analyse des données. T2'' (leçons sur l'analyse factorielle et la reconnaissance des formes et travaux du Laboratoire de statistique de l'Université de Paris 6. T. 2 : l'analyse des correspondances), [[Éditions Dunod|Dunod]] Paris Bruxelles Montréal, 1973
* [[Jean-Paul Benzécri]] et Al. ''Pratique de l'analyse des données. T1'' (analyse des correspondances. Exposé élémentaire), Dunod Paris, 1984, 
* [[Jean-Paul Benzécri]] et Al. ''Pratique de l'analyse des données. T2'' (abrégé théorique. Études de cas modèle), Dunod Paris, 1984
* {{Ouvrage|auteur1=[[Brigitte Escofier-Cordier|Brigitte Escofier]]|auteur2=Jérôme Pagès|titre=Analyses factorielles simples et multiples ; objectifs, méthodes et interprétation|éditeur=Dunod, Paris|année=2008, {{4e}}{{éd.}} ({{1re}}{{éd.}} 1988)|pages totales=318|isbn=978-2-10-051932-3|id=Escofier&Pagès}}
* {{Ouvrage|auteur1=François Husson|auteur2=Sébastien Lê|auteur3=Jérôme Pagès|titre=Analyse des données avec R|éditeur=[[Presses Universitaires de Rennes]]|année=2009|pages totales=224|isbn=978-2-7535-0938-2|id=Husson et al}}
* Ludovic Lebart, Morineau Alain, Piron Marie; ''Statistique exploratoire multidimensionnelle'', Dunod Paris, 1995
* {{Ouvrage|auteur1=Jérôme Pagès|titre=Analyse factorielle multiple avec R|éditeur=EDP sciences, Paris|année=2013|pages totales=253 |isbn=978-2-7598-0963-9|id=Pagès 2013}}
* Mathieu Rouaud ; [http://www.incertitudes.fr/proba-stat-acp/livre.pdf Probabilités, statistiques et analyses multicritères] Un livre de 290 pages qui traite de l'ACP (les principes et de nombreux exemples concernant, entre autres, les isolants thermiques et les eaux minérales). Version numérique libre et gratuite.
* [[Michel Volle]], ''Analyse des données'', [[Economica]], {{4e|édition}}, 1997, {{ISBN|2-7178-3212-2}}

=== Liens externes ===
* [http://factominer.free.fr/ FactoMineR, une bibliothèque de fonctions R destinée à l'analyse des données]

{{Portail|probabilités et statistiques}}

[[Catégorie:Analyse des données]]
[[Catégorie:Exploration de données]]