Bonjour,
Voici donc quelques petites précisions :
>Quelle est la fiche de données ? Il y en a deux.
J'ai importé les données à partir de R pour être sûre d'avoir le même jeu
de données :
sous R, la fonction a utilisé est la suivante :
> Rtoade4(meau$mil)
File creation meau$mil.txt
File creation meau$mil_col_lab.txt
File creation meau$mil_row_lab.txt
Je transforme le fichier en binaire :
TextToBin : Text -> Binary
Text input file : C:/.../meau$mil.txt
Binary output file : meaumil
*------------------------------------------------------------*
| WinADE-4 * Metrowerks CodeWarrior C * CNRS-Lyon1 * JT & DC |
| TextToBin: Text->Binary 24/06/04 09/22 |
*------------------------------------------------------------*
-----------------------------------------------------
Option: Excel-Binary
Input text file: C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meau$mil.txt
-----------------------------------------------------
Row 1 Words 10 : 10×41×8.5×295×110×2.3×1.4×0.12×3.4×0.11
Row 2 Words 10 : 11×158×8.3×315×13×7.6×3.3×2.85×2.7×1.5
Row 3 Words 10 : 11×198×8.5×290×113×3.3×1.5×0.4×4×0.1
Row 4 Words 10 : 12×280×8.6×290×126×3.5×1.5×0.45×4×0.73
Row 5 Words 10 : 13×322×8.5×285×117×3.6×1.6×0.48×4.6×0.84
Row 6 Words 10 : 11×303×8.5×245×100×1.7×0.9×0.05×2.7×0.16
Row 7 Words 10 : 13×62×8.3×325×95×2.3×1.8×0.11×3×0.13
Row 8 Words 10 : 13×80×7.6×380×20×21×5.7×9.8×0.8×3.65
Row 9 Words 10 : 15×100×7.8×385×46×15×2.5×7.9×7.7×4.5
Row 10 Words 10 : 16×140×8×360×76×12×2.6×4.9×8.4×3.45
...
240 numbers found in file C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meau$mil.txt
Row number: 24 Column number: 10
-----------------------------------------------------
Output binary file C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil
-----------------------------------------------------
Je crée un fichier texte contenant 1 colonne et 4 lignes afin de définir
les blocs:
6
6
6
6
Je le transforme en binaire :
TextToBin : Text -> Binary
Text input file : C:/.../bloc.txt
Binary output file : bloc
*------------------------------------------------------------*
| WinADE-4 * Metrowerks CodeWarrior C * CNRS-Lyon1 * JT & DC |
| TextToBin: Text->Binary 24/06/04 09/13 |
*------------------------------------------------------------*
-----------------------------------------------------
Option: Excel-Binary
Input text file: C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\bloc.txt
-----------------------------------------------------
Row 1 Words 1 : 6
Row 2 Words 1 : 6
Row 3 Words 1 : 6
Row 4 Words 1 : 6
4 numbers found in file C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\bloc.txt
Row number: 4 Column number: 1
-----------------------------------------------------
Output binary file C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\bloc
-----------------------------------------------------
>quelles sont les indications utilisées comme indicateurs de blocs ?
KTabUtil : InitKTab
Matrix Input File : meaumil
Row indicator : bloc
Col indicator :
Option : output file name
*------------------------------------------------------------*
| WinADE-4 * Metrowerks CodeWarrior C * CNRS-Lyon1 * JT & DC |
| KTabUtil: InitKTab 24/06/04 09/23 |
*------------------------------------------------------------*
Option: InitKTab
Input file: C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil
Row indicator file: C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\bloc
Col indicator file:
Output ASCII file: C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil!.ktp
4:6/6/6/6/
1:10/
-----------------------------------------------------
Qualitative variable file: C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil!TLl
Number of rows: 24, variables: 3, categories: 34
Auxiliary ASCII output file
C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil!TLl.label: label file
----------------------------------------------
Qualitative variable file: C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil!TCc
Number of rows: 40, variables: 3, categories: 24
Auxiliary ASCII output file
C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil!TCc.label: label file
----------------------------------------------
Qualitative variable file: C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil!T4a
Number of rows: 16, variables: 2, categories: 8
Auxiliary ASCII output file
C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil!T4a.label: label file
----------------------------------------------
Normalisation du tableau :
KTabUtil : NormelizeKTab
---.ktp input file : meaumil!.ktp
matrix input file : meaumil
row weight (default : 1/nk) :
column weight (default : 1/pj) :
Row-Col option : 3
Option output file name :
*------------------------------------------------------------*
| WinADE-4 * Metrowerks CodeWarrior C * CNRS-Lyon1 * JT & DC |
| KTabUtil: NormelizeKtab 24/06/04 09/25 |
*------------------------------------------------------------*
Option NormelizeKtab
Input file C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil
Number of rows: 24, columns: 10
------------------------------------------
Each row has uniform weight (Sum inside block = 1)
File C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil.ktpl contains the row
weights
It has 24 rows and 1 column
Each column has uniform weight (Sum inside block = 1)
File C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil.ktpc contains the
column weights
It has 10 rows and 1 column
------------------------------------------
File C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil.ktta contains the
transformed table
block-centered and normalized by columns
It has 24 rows and 10 columns
---------------------------------------------------------
Bouroche, J.M. (1975) Analyse des donnŽes ternaires:
La double analyse en composantes principales.
Thse de 3. cycle, UniversitŽ de Paris VI. 1-57 + annexes.
---------------------------------------------------------
>Peut-on comparer les tableaux normalisés par blocs avant toute autre
>observation ?
>A-t'on les mêmes tableaux, les mêmes poids sur les deux marges, la même
>division ?
Et je pense que mon problème vient effectivement de là... la normalisation
des tableaux...
Le poids des colonnes :
-------------------------------------------------
Binary input file: C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil.ktpc - 10
rows, 1 cols.
1 | 0.1000
2 | 0.1000
3 | 0.1000
4 | 0.1000
5 | 0.1000
6 | 0.1000
7 | 0.1000
8 | 0.1000
9 | 0.1000
10 | 0.1000
Alors que sur R, j'obtiens :
> wit1$cw
[1] 1 1 1 1 1 1 1 1 1 1
Le poids des lignes :
-------------------------------------------------
Binary input file: C:\ade4\ade-4\ade4tools\dir_Try\Meaudr\meaumil.ktpl - 24
rows, 1 cols.
1 | 0.1667
2 | 0.1667
3 | 0.1667
4 | 0.1667
5 | 0.1667
6 | 0.1667
7 | 0.1667
8 | 0.1667
9 | 0.1667
10 | 0.1667
11 | 0.1667
12 | 0.1667
13 | 0.1667
14 | 0.1667
15 | 0.1667
16 | 0.1667
17 | 0.1667
18 | 0.1667
19 | 0.1667
20 | 0.1667
21 | 0.1667
22 | 0.1667
23 | 0.1667
24 | 0.1667
Sur R :
> wit1$lw
[1] 0.04166667 0.04166667 0.04166667 0.04166667 0.04166667 0.04166667
[7] 0.04166667 0.04166667 0.04166667 0.04166667 0.04166667 0.04166667
[13] 0.04166667 0.04166667 0.04166667 0.04166667 0.04166667 0.04166667
[19] 0.04166667 0.04166667 0.04166667 0.04166667 0.04166667 0.04166667
J'avais cru comprendre que l'on pouvait réaliser une normalisation de
Bouroche sur R à l'aide de l'option "total" du paramètre "scaling" de la
fonction "within.pca", mais apparemment les options de pondérations des
lignes et des colonnes ne sont pas les mêmes dans les 2 programmes...
Pour replacer ma question dans son contexte, je désire réaliser une analyse
Statico à l'aide de R.
Je voudrais coupler une série de K-tableaux faunistiques relevant d'AFC
avec une série de K-tableaux environnement relevant d'ACP. Se pose à ce
stade le problème de pondération des lignes qui doivent être les mêmes dans
les 2 séries de K-tableaux (Cf. mon mail du 27 mai 2004 :
http://pbil.univ-lyon1.fr/ADE-4/adelisthtmlannuel/03/0283.html).
Je voulais également obtenir les coordonnées des lignes-stations communes
aux 2 tableaux, ces sorties n'étant pas fournies pas R (Cf. mail du 1/06/04
: http://pbil.univ-lyon1.fr/ADE-4/adelisthtmlannuel/03/0286.html ).
A ce sujet Pierre Bady m'a fait passer une fonction "maison" (mail du 1er
juin : http://pbil.univ-lyon1.fr/ADE-4/adelisthtmlannuel/03/0287.html )
permettant de calculer ces coordonnées, mais comme il n'était pas très sûr
du calcul, j'ai essayer de comparer les 2 analyses réalisées avec ADE4 et
R, et j'ai commencé par comparer les 2 analyses STATIS...
Donc pour résumer,
- est-il possible de coupler un K-tableau relevant d'AFC et un
K-tableau relevant d'ACP ?
- comment calculer les coordonnées des lignes-stations communes
aux 2 séries de K-tableaux ? (le calcul de P.Bady est-il correct ?)
merci
Charline Laurent
This archive was generated by hypermail 2b30 : Tue Sep 07 2004 - 13:30:57 MEST