{"word_count": 558, "character_count": 4157}
Hier soir, on voyait comment faire une analyse des correspondances à partir du tableau de contingence. Mais on peut avoir le problème autrement, à partir des individus. Certes, ces derniers ne sont pas observés (vraiment), mais peu importe.&#62; data(HairEyeColor)&#62; N = HairEyeColor[,,"Male"] + HairEyeColor[,,"Female"]&#62; Hair = rep(rep(rownames(N),ncol(N)),+ as.vector(N))&#62; Eye =  rep(colnames(N),apply(N,2,sum))&#62; df = data.frame(Ind=1:sum(N),Hair,Eye)&#62; tail(df)    Ind  Hair   Eye587 587 Blond Green588 588 Blond Green589 589 Blond Green590 590 Blond Green591 591 Blond Green592 592 Blond GreenOn a ici une base classique, mais pas encore une matrice avec laquelle on pourra faire de l'analyse de données. Pour récupérer la matrice, on va disjoncter cette base,&#62; XH &#60;- model.matrix(Ind~0+Hair,data=df)&#62; XE &#60;- model.matrix(Ind~0+Eye,data=df)&#62; X &#60;- cbind(XE,XH)&#62; tail(X[,1:5])    EyeBlue EyeBrown EyeGreen EyeHazel HairBlack587       0        0        1        0         0588       0        0        1        0         0589       0        0        1        0         0590       0        0        1        0         0591       0        0        1        0         0592       0        0        1        0         0On a cette fois une vraie matrice de nombres, avec laquelle on peut faire des calculs. Parmi les représentation usuelles, on peut définir la matrice de Burt, qui est une forme de tableau de contingence en grande dimension (à partir des variables disjonctées).&#62; Burt &#60;- t(X) %*% X &#62; Burt[,1:5]     EyeBlue EyeBrown EyeGreen EyeHazel HairBlackEyeBlue  215        0        0        0        20EyeBrown   0      220        0        0        68EyeGreen   0        0       64        0         5EyeHazel   0        0        0       93        15HairBlack 20       68        5       15       108HairBlond 94        7       16       10         0HairBrown 84      119       29       54         0HairRed   17       26       14       14         0(etc) Mais on en reparlera.On avait noté que pour l'Analyse des Correspondances, la distance naturelle était celle du chi-deux. Prenons notre tableau disjonctif. L'écriture classique de la distance du chi-deux entre deux lignes,donne ici (comme les sommes des lignes donne  où  est le nombre de variables qualitatives - ici 2, mais on espère bien pouvoir intégrer plus de variables qualitatives - et où les nombres qui apparaissent dans la matrice sont les ),avecVoilà la distance entre deux individus. Si on voulait faire un ACP, on utilisait comme distance entre individusBref, notre analyse des correspondance peut être obtenue en faisant une ACP, avec les bons poids. Étant donné quepour voir l'analyse des correspondances comme une ACP (entre des vrais individus, contrairement à hier), on va utiliser les poids suivants, sur les colonnes (i.e. les variables indicatrices disjonctées), i.e. on veut utiliser des distances de la formeavec ici Bon, lançons nous... soit on utilise juste le terme de gauche et on demande une ACP sur données normalisées, soit on utilise ce poids, avec une ACP sur données brutes. Essayons le premier.Pour les poids, on utilise&#62; n=sum(N)&#62; (nk=apply(X,2,sum))  EyeBlue  EyeBrown  EyeGreen  EyeHazel HairBlack      215       220        64        93       108&#62; n*p[1] 1184On pose alors&#62; omega=(n-nk)/(n*p)et on fait ensuite une régression, dont on va chercher à représenter les variables (qui sont ici toutes nos modalités)&#62; library(FactoMineR)&#62; PCA(X,col.w = omega)On retrouve ici - à une symétrie près - la sortie que l'on a quand on fait une analyse des correspondances,&#62; library(FactoMineR)&#62; acm=CA(N)On a ici une approche très différente de celle d'hier:hier, on faisait une ACP, en projetant les individus, qui étaient les profils lignes de notre tableau de contingence (puis une analyse des profils colonnes, qu'on superposait ensuite)aujourd'hui, on fait une ACP, on projetant cette fois les variables, dans un tableau disjonctif complet.L'avantage de cette dernière technique est qu'elle peut se faire avec plus de deux variables qualitatives, contrairement à la première approche.