Re: ACP ADE4-SAS

From: Cecile Rebourg (rebourg@moulon.inra.fr)
Date: Mon Jun 07 1999 - 11:42:43 MET DST


Bonjour,
Merci de votre reponse. Apparemment, mon message precedent a ete mal percu, et je
m'en excuse. Il ne s'agissait en aucun cas d'un commentaire affirmant que l'ACP sur
matrice de correlations dans ADE4 etait differente de celle sur SAS, mais je
cherchais une explication possible aux differences que j'observais avec mes donnees.

Par ailleurs, je m'etais bien sur assuree de comparer des sorties comparables avant
de poser cette question. Les sorties de SAS sont directement comparables avec celles
d'ADE4.
Pour l'exemple que vous citez, j'obtiens avec SAS :

Eigenvalues of the Correlation Matrix: Total = 6 Average = 1

                       1 2 3
4 5 6
Eigenvalue 4.6560 0.9152 0.2404 0.1027 0.0647 0.0210

Difference 3.7408 0.6748 0.1377 0.0381 0.0437
Proportion 0.7760 0.1525 0.0401 0.0171 0.0108
0.0035
Cumulative 0.7760 0.9285 0.9686 0.9857 0.9965 1.0000

Factor Pattern

           FACTOR1 FACTOR2 FACTOR3 FACTOR4 FACTOR5 FACTOR6

VAR1 0.95849 0.03249 -0.19684 -0.01603 -0.20306 -0.00157
VAR2 0.89438 0.40297 -0.01940 -0.15695 0.07798 -0.08130
VAR3 0.74128 0.63469 0.18139 0.10253 0.00183 0.06516
VAR4 0.92854 -0.24411 -0.23755 -0.03947 0.12017 0.07613
VAR5 0.92830 -0.28278 0.02157 0.22814 0.04220 -0.06333
VAR6 0.81485 -0.45759 0.33393 -0.11707 -0.03340 0.01720

Ces deux sorties sont identiques aux cnvp et cnco.
La seule sortie que je ne suis pas arrivee a rabouter a celle d'ADE4 est celle
appelee "Standardized Scoring Coefficients"., que voici :

 FACTOR1 FACTOR2 FACTOR3 FACTOR4 FACTOR5 FACTOR6

VAR1 0.20586 0.03550 -0.81869 -0.15609 -3.14068 -0.07503
VAR2 0.19209 0.44030 -0.08068 -1.52806 1.20603 -3.87867
VAR3 0.15921 0.69348 0.75445 0.99821 0.02826 3.10886
VAR4 0.19943 -0.26673 -0.98802 -0.38427 1.85857 3.63199
VAR5 0.19938 -0.30897 0.08970 2.22117 0.65272 3.02133
VAR6 0.17501 -0.49998 1.38890 -1.13981 -0.51660 0.82054

OBS FACTOR1 FACTOR2 FACTOR3 FACTOR4 FACTOR5 FACTOR6

1 -0.91775 0.32739 1.05936 -1.24437 -0.45528 0.74498
2 -0.35320 -0.13600 0.88154 0.65112 -0.86798 0.01198
3 -0.58515 -0.44432 0.93634 -0.58031 -0.66312 1.34690
4 -0.12423 -0.47551 0.37212 1.86705 0.10710 0.19177
5 0.08189 -0.65198 -0.12889 1.93912 -1.26282 0.88541
6 -0.22895 -0.21244 0.30158 1.26886 -0.23720 1.18506
7 1.79006 0.85484 -1.02934 -0.42458 -0.46368 0.97586
8 1.44901 0.78940 -0.02825 0.02318 -0.17757 -0.44063
9 1.56049 0.62343 1.25401 -0.59004 0.31099 -0.27614
10 0.68033 -0.79374 1.02216 1.22411 1.14438 -0.60219
11 0.26610 0.14650 0.69644 0.58379 0.57407 -0.68333
12 0.33542 -0.44472 0.21867 0.97449 0.41640 -0.47122
13 0.77406 1.39447 -1.95150 -0.48249 0.33068 -1.04127
14 0.64130 1.11469 0.29205 0.10124 2.09739 -0.36764
15 1.26806 -0.12456 -0.76968 -0.13699 -2.85300 -0.47852
16 0.41739 -0.90922 0.51882 -1.21548 0.58198 0.85729
17 -0.01123 -1.86283 -2.49284 -0.29293 1.64627 1.83956
18 0.55440 -2.44037 0.60077 -2.04837 -0.07341 -1.72793
19 -1.58737 -0.92538 -0.14155 -0.35993 -0.10422 -1.22936
20 -1.70479 -0.01526 -1.01419 0.43746 0.28154 -1.50836
21 -1.18810 0.42474 -1.59519 -0.00570 -0.94485 -0.80252
22 -1.03837 1.51322 0.20134 0.33654 0.83269 -0.43889
23 -0.87700 0.92241 -0.09492 -1.08337 0.32486 1.68403
24 -1.20237 1.32526 0.89115 -0.94241 -0.54524 0.34516

Concernant les variables, j'ai vu que cela correspondait aux valeurs du fichiers
cnco divisees par la valeur propre du facteur en question, mais en ce qui concerne
les observations, je n'ai pas trouve. Si vous pouviez m'eclairer sur ce point, ...

Quoiqu'il en soit, je pense avoir resolu mon probleme. Il semble en fait que les
differences que j'observais etaient dues aux donnees manquantes. Les observations
presentant une donnee manquante pour n'importe quel caractere sont systematiquement
eliminees sous SAS. Il semble que sous ADE4, elles soient remplacees (j'avais des
points comme donnee manquante) par un 0, avant que ce soit centre-reduit. Est-ce
cela?

Encore merci pour votre aide

Cecile Rebourg

Daniel Chessel a écrit :

>
> L'ACP sur matrice de corrélations dans ADE-4 est exactement la même que
> celle de S-PLUS
> Mais pas celle de SAS ? Je ne permettrais pas ce genre de commentaire
> Mais peut-être c'est plus difficile qu'il n'y parait de comparer deux programmes
> Cordialement
>
> > Bonjour,
> > A partir d'un meme fichier, une ACP sur matrice de correlation
> > realisee sous ADE4 ne me donne pas le meme resultat que la
> > meme ACP realisee sous SAS.
> > Avez-vous une explication a cette difference ?
> > Merci beaucoup
> >
> > Cecile Rebourg
> Daniel Chessel
> ---------------------------------
> Universite Lyon 1 - Bat 401C
> 69622 Villeurbanne CEDEX - France
> Tel : 04 72 44 82 77
> Fax : 04 72 43 11 41
> ---------------------------------



This archive was generated by hypermail 2b30 : Sat Feb 10 2001 - 10:35:59 MET