Download Programmes Fortran pour distances de rangs, constellations et

j
i
1
GRAMMES F
POURDISTANCES
DE RANGS
CONSTELLATIONS
ET CORRÉLATIBN
- .... .
.-
~
__
. .
.
._
OFFICEDE LA RECHERCHESCIENTIFIQUE
-
.~
.
~
.
.-
Anne-Marie AUBRY
INITIATIONS
-
DOCL~MENTATIONSTECIINIWIES
No 30
O . R. S. T. O . M .
P A R I S
1976
PROGRAMMESFORTRAN
POURDISTANCES
DE RANGS
CONSTELLATIONS ET CORRELATION
Anne-Marie AUBRY *
Analyste-Programmeur
*
Banque de données pédologiques de I'O.R.S.T.O.M., 70-74, route d Aulnay, F 93140 Bondy
SOMMAIRE
Page
RESUME
- ABSTRACT -
7
RESUMEN
9
INTRODUCTiON
RAPPEL DES METHODES
11
PROGRAMMES D K A T ET H M V T POUR DISTANCES ET CONSTELLATIONS-
13
MODE D'EMPLOI PROGRAMMES
DES
i- H
AT M
K D
VT
CARTES CONTROLE
D'UN FORME
25
25
FICHIER D'ENTRCE STANDARD
27
EXEMPLE
PROGRAMME
AT
CS
POUR CORRELATIONS
M 0 3 E D'EMPLOI D U PROGRAMME
T
CS A
D
35
43
CARTES C ~ N T R O L E
43
EXEMPLE D~ENTREE
43
EXEMPLE DE SORTIE
44
PROGRAMME
E
M
D'EMPLOI
MODE
POUR MEDIANES
DU PROGRAMME
E
M
CARTES CONTROLE
45
DT
49
49
EXEMPLE
51
EXEMPLE DE SORTIE
BIBLIOGRAPHIE
55
7
II s’agit de fournir aux utilisateurs le texte et les moyens d’emploi de programmes, indépendamment d‘une application donnée. Ces programmes sont écritsen Fortran, langage de programmation
universel e t évolué. On en donne le mode d‘emploi complet, y compris :
1. les cartes contrôle système à l‘intention des utilisateurs du terminal de I’O.R.S.T.0.M.;
2. la façon dont doivent être présentées les données, quels que soient l’ordinateur e t
l’application.
Cette notevient en complément des notesantérieures qui ont décrit des cas concrets
elle devrait permettre aux lecteurs intéressésdes’enservirdans
d’utilisation decesprogrammes;
leur propre discipline.
ABSTRACT
A small statistical package including rank distances, a clustering procedure, rank correlations
and confidence intervals for the median is presented.
The computer programs and their control cards are given for a wide variety of applications.
The language used is FORTRAN I V and the processor is a UNIVAC 1108, 192K. Al1 F-type formats
are accepted. Missing data are allowed. For users working on the ORDO terminal of O.R.S.T.O.M.,
system control cards are added.
This report follows some papers on specific applications of the same programs.
RESUMEN
Se presenta un conjunto deprogramasde
calculo dedistancias,constelaciones,correlaciones
la estadistica matematica no paramétrica para
e intervalos de confianza demedianas.
el
Por su gran nlimero de aplicaciones sonprogramasdeinteresgeneral.Estanescritosen
FORTRAN I V para un UNIVAC 1108, 192K; y aceptantablas incompletas de datosescritosen
forma F.
Se dB la forma de aplicacidn completa e inclus0 las cartas “control sistema” para el utilizador
que trabaja sobre el terminal ORDO del O.R.S.T.O.M.
Este trabajo es el complemento deostrosanteriores,ya
aplicacidn especifica.
publicados, que se refieren a una
9
INTRODUCTION
Les quatreprogrammes
D K A T , H M V T, C S A T e t M E D T sontl'expression en
langage Fortran d'algorithmes de calcul qui ont été décritspar VAN DEN DRIESSCHE (1974).
Ce sont des méthodes multivariables non-paramétriques,c'est-à-dire
permettant d'étudier des
fichiers dedonnées indépendamment de la fonction de répartition desvariables, qui est souvent
uneinconnue.
D K A T calcule les distances entre m unités prises 2 à 2 dans un espace à v dimensions
(v étant le nombre de variables).
H M V T fait le regroupement en constellations deces unités (((clusteringn).
C S A T fait le test de corrélation de Spearman entre les v variables prises 2 B 2 du fichier
des m unités.
ME
DT
donne,parvariable,
la médiane e t les limites de confiance de la médiane.
Ces programmespeuventintéresser
I'hydrobiologie, entre autres disciplines.
la pédologie, la géographie, la foresterie, la géologie,
Ils sontutilisables sur le terminal OrdoprocesseurdesServices
ScientifiquesCentraux de
1'O.R.S.T.O.M. à Bondy. Ce terminal est relié à un ordinateur Univac 1108 qui exécute les travaux.
Bien entendu, le mode d'emploi et les cartes contrôle système restent valables pour travailler dans
le site de l'ordinateur central de la société S.T.A.D. (11. rue de la Vistule, 75013 Parisi.
Les valeurs limites des paramètres sontimposées par la t a i l l e de mémoire centraledisponible :
1061< mots de 36 bits. D K A T occupe 72K. H M V T 89K, C S A T 60K e t M E D T 57K. On
pourrait très facilement réduire ces occupations mémoire pour des applications moins volumineuses
(v. infra), e t qu'il y aurait lieu de faire surdes ordinateursplus petits que
quecellesautorisées
ce 1108. D'autre part, l'utilisation sur un IBM desséries 360 ou 370 ne poserait aucun problème,
car il suffirait de ne modifier qu'untrès petit nombre d'instructions.Enfin, la présentation des
textes Fortran permetaux utilisateurs extérieurs à 1'O.R.S.T.O.M. de les faire perforer, car le
nombre d'instructions parprogrammenedépassepas
500.
11
RAPPEL.DES MÉTHODES
Les programmes permettent de traiter des fichiers de données numériques positives ne satisfaisant aucune condition préalable. Les variables peuvent ëtre quantitatives pures
et/ou ordonnées
et/ou binaires. On accepte que des données manquent, en les repérant par le code -1.
D K A T calcule des ((distances)) - plus exactement des ((relations de dissemblance)) -entre
qui les caractérisent). Leprocessusde
calcul
proposépar VAN DENDRIESSCHE et GARCIA GOMEZ (1972) est le suivant :
m points d'un espace à v dimensions(desvariables
a) Remplacementde la matrice X desdonnées - brutes ou calculées - xki (1 G Ic < m.
1 < i < v) par une matrice R de rangs rki de mêmes dimensions. Le rangement est fait pour chaque
variable i indépendamment :
- remplacement desdonnées
effectives par les rangs 1 à m i (mi < m, mi = m
manquantes de i, mi est donc le nombre de points où a é t é relevée i);
- remplacement desrangs
correspondant auxdonnées
- données
identiquespar la moyenne de ces
rangs;
- les donneesmanquantesgardent
le code-1.
b ) Pour chaque i,on calcule un facteur de pondération fonction dunombre e t de la répartition
des données :
~i
[
1
ei
C
= mi3 -mi--.
(t,:
-tqi) /12
qi= 1
où : mi est l'effectif réelde i (nombre de points où i a été relevée), mi
ei est le nombre de lots qi dedonnées identiques (1 ei Q mi);
tqi est le nombre dedonnées du lot qi (mi 2 tqi
1).
>
<
< m;
Remarque : valeurs limites de li
-
qi, tqi = 1 : toutes les données sontdistinctes.
mi3 -mi
Alors ei = mi e t
Il
=-
- 3 qi,
12
tqi = mi : toutes les donnéessontidentiques,
sur I'échantillon des points étudiés.
Alors ei = 1 et li = O.
la variable i est ((constante))
12
c) L'expression de la distance entre 2 points d'indice h e t k (h # k, 1 Q h S m, 1 Q k G m)
adaptée de KENDALL e t STUART (1966)s'écrit alors :
dhk est un nombre réel de l'intervalle [0,2]
Si, pour un i donné, li = O ou l'un au moins des rhi ou rki = - 1, le ième terme est abandonné
e t la constante v est diminée de 1 dans l'expression (1). A la limite, si v devient êgal à O, la distance
est incalculable et on la note II.
H M V T regroupe en constellations m points sur la basedes C: distances les séparant,
La méthode imaginée par VAN DEN DRIESSCHE (1964)est décrite en détail (1974)et également
traitée entièrement sur un exemple (1965)
où les distances employéessont les D2de MAHALANOBIS.
L e processus est valable quel que soit le type des distances, mais on a adjoint dans ce texte
un enchaînement automatique sur les distantes de rangs D K A T.
II est à noter, pour mémoire, que H M V T peut aussi calculer entre m points les distances
Ag de HIERNAUX (1965)en vue de leur regroupement réalisé dans le mëme travail. Divers travaux
a, 1968 b).
ont fait appel à cet algorithme (1966, 1968
Les deux règlesabrégées,
peuvents'écrire :
à pas, les constellations,
à partir desquelles sont réalisées,pas
- Règ/e
7. Distance moyenne entre 2 constellations supérieureaux deux distances moyennes
intra-constellation.
- Règle 2.
Distanceentre points non regroupés e t toute constellation supérieure à toutes
les distances moyennes' intra-constellation.
C S A T calcule le test de corrélation derangsde
prises 2 à 2 caractérisant les m points de l'êtude.
SPEARMAN (1904)entre les v variables
a) On part de la matrice des rangs obtenue comme pour le calcul des distances (variable par
variable indépendamment).
b ) Puis, pour tout couple de variables d'indice i e t j (1 < i G v, 1 < j S v , i # j), on calcule les
facteurs de pondération Ji e t li :
où : mii est le nombrede points où i e t j ont été relevéesenmëmetemps;
ei, tqi pour i.
qi, ei, tqi sont les homologuesdeqi,
c) Enfin, le cœfficient de corrélation entre i et j est donnépar
la formule :
rki # - 1
rkj # - 1
Rii est un nombre réel appartenant à l'intervalle [ -1, 11. On compare ce nombre à la valeur
critique, facteur de l'effectif mii, lue dans la table extraite deBEYER (1966)e t incluse dans l e
programme, e t ondonne le résultatdu test au risque 5 76 ou 1 %.
On a adjoint le programme M E D T qui détermine au risque 5 % les limitesinférieure
(notée 1 ) et supérieure (S)de l'intervalle de confiance desmédianesdechaquevariable.
13
PROGRAMMES D K A T ET H M VT
POURDISTANCES ET CONSTELLATIONS
D K A T : DISTANCESDERANGS
IMPLICIT INTEGEcl(V)
IMPLICIT IPJTEGER( 5 )
DIMENSION ECH(13) (VAR(13) qV(13) rFICH(H) rVD(400)
DIMENSION V E Q 1 ( 1 5 ) ~ V E Q 2 ( 1 5 ) r V E Q 3 ( 1 5 )
EQUIVALEMCE (VEQl(1)*l'UT) ( V E Q l ( 2 ) rVDEM) 9 (VEQl(3)rVAK(1) 1
EQUIVALENCE (VEQ2(1) * B I D I ) 9 (VEO2(7)r V O E M ) 9 (VE!3?(3) 9\/4@(1))
EQUIVALENCE (VEQ3(l)rRIDl),(VE03(2)rRIn2)1(V~~3(3)rVAR(l))
COM~ON/C2/DI~100~100~rC(100~
./C3/Y(40O)/Cl/IWEAD~IWRIT~I~2~KFOR
./C4/JE*I*
JPEtIV
COMWON X(3000)~XRIS(42000)
DATA CP/'.'/
CALL PRTM(hbeO*O)
DEFINE FILE lO(S00~15tUtID)
XBIS(l)=O.
ID=1
READ(5r1003)IREAD9IWRIT1IW2
1003
1001
20
35
37
21
22
23
24
FORMAT (312)
READ(IREAD~20~END=lOOO)FICH
FORMAT ( R X q R A l )
READ(IREAD*35)TRAV*KFOR
FOKMAT(8XtAlrI2)
IF IKFOR.EQ.O)KFOR=I
800
IF(KFOR.CT.lO)GOTO
READ(IREAD*lOO)JE*IVM
IF(JE.LT.2.OR.JE.GT.lOO)GOTO
600
I F ~ I V ~ ~ L T ~ l ~ O R ~ I V M ~ ~ T601
~4OO~GDTO
HIOl=O
HID2=O
VDEM=l
VUT=O
DO 37 I=l*IVM
VD ( 1 1 =VDEM
DO 21 J=l*JE
READ(IREAD922)ECH
dRITE(1O'J)ECH
CONTINUE
FORtJAT(13A4)
DO 23 I = l * I V M
KEAD(IRE4Dv22)VAR
JEPI=JE+I
WRITE(lO'JSPI)VEB1
CONTINUE
FOKPAT(2FZ.Ot13A4)
READ(IREAD*25)IV
14
25
(14)
IF(IV.LT.1.OW.IV.GT.IVM)GOTO
IF ( I V . E O a 1VY)GDTO 500
FOYMPT
602
I=O
IVl=IVM+l
D O 34 L=lrJV
HEAD(lREAD922)V
24
I=I+1
I F ( I .EQ.IVI)GOTO
ZY
2 99
1002
JEPI=JE+I
QEAD(lO*JSPI)VEQ3
DO E8tLY=1r 1.3
IF(VAH(K).NE.V(K))GOTO
289
CONTINUE
GOTO 34
VDEM=O
NRITE(lO*JFfI)VEQE
VD (1 ) =VDEY
GOTO 26
29
FORMAT(~XIF~.~V~~A~)
CONTINUE
34
IF(I.EQ.IVV)GOTO
500
IPl=I+l
VDEM=D
DO 38 L=IPlrIVM
VO (L)=VDEM
JEPI=JE+L
READ(lO*JEPI)VEQ3
VDEF?=O
WRITE(lO'JFPI)VEQ2
READ ( 1 O * JEPI ) V E Q L
38
CONTINUE
500 JPE=O
CALL L€CT(IVYrJE~XrVO~$lOOO)
502
ID=JE+l
I=O
O 0 44 M=lrIVM
JEPI=JE+M
REAO(lO*JEPI)VEQ2
IF(VC)EM.EQ.O)GOTO
44
I=I+l
IF(I.6T.IV)GOTO 45
CALL RGEQLI (XJ)
IF (XJI4r6094
VUT=l
60
WRITE(lO*JEPI)VEQl
Y(I)=XJ
4
CONTINUE
A4
CONTINUE
45
JEMl=JE-l
DO 7 K=l*J€Ml
KPl=K+l
DO 6 L=KPlrJE
D=O
fil
30 1
62
401
63
JV=O
DO 4 I=lrIV
IF(Y(I))61,5*61
Kl=K-1
TK=X(Kl*IV+I)
IF(TK+l.)6?*5rh2
Ll=L-l
TL=X(LI*IV+I)
IF (TL+1s ) 6 3 9 5 . 6 3
DKL=TK-TL
O=O+OKL*DKL/Y(I)
JV=JV+l
15
5
CONTINUE
IF(JV.EQeO)GOTO
DI (K*L)=D/JV
8
GOTO h
n I (K*L)=11.
6 DI (LqK)=DI (KIL)
7 CONTINUE
FORMAT (1H1)
104
WRITE(IWRITI~OS)TRAV,FICH
105
FORMAT(lH1*3X*'.TRAVAIL '*Al,' SUR LE FICHIER DU @,8A1,63X*11*
./4%~'.'*9~*10~'1'~~10~'2'~*10~'3'~,10~*4'~,1~~*5'),~0('~')~~
~~lO~'A'~*~~~'~*~*'0'/4XI,.(r10(r1234567~9O~~/lX~lO4(~
DO 12 K=lrJE
DO 1 1 L=l*JE
IF(L.NF.K)COTO
10
C (L)=cp
GO10 11
CALL NOTA (K,L)
10
CONTINUE
11
8
WRITE(IWRITI~O~)K,(C(L)~L=~~JE)
12
13
106
100
102
17
CONTINUE
IF(K.EQ.100)GOTO
17
K=K+1
DO 13 J=K*100
~RITE(IWRIT~lO6)J
FORMAT(lX*I3,'.')
FOHMAT(I3rI4)
FORMAT(lX~I3~'.'*100Al)
CONTINUE
WRITE(IWRIT,109)TRAV,FICH
1O 9
108
15
111
112
36
113
114
19
115
18
110
116
16
117
FORMAT(1Hlr' ENTREES LIGNES ET COLONNESDE LA
.UTILISEES POUR LE TRAVAIL '*Al*' EFFECTUE SUR
~8Al/lX*llFJ('.'))
WRITE(IWRITtlO8)
/ * ' ENTREESLIGNES 1 )
FORMAT(
DO 15 J=lrJE
READ(1O'J)ECH
WRITE(IWRITrlll)J,ECH
FORYAT(lX,T3*'.'*13A4)
WRITE(IWRIT~ll2)
FORMAT( 1Hlr' ENTREES COLONNES
*)
DO 18 I=lrIVM
JEPT=JE+I
READ(1O'JEPI)VEQl
IF(VDEM.EQ.O)GOTO
19
IF(VUT.EQ.1)GOTO
36
WRITE(IWRIT*113)VAR
GOTO 18
WRITE(IWRIT,1141VAR
FORMAT(4X,'UTILISATION DE'r3Xv13A4)
FORWT(4X~'ABANDON DE'97X113A4)
GOTO 18
WRITE(IWRIT,115)VAR
FORMAT(21X913A4)
CONTINUE
WRITE(IWRITv104)
W R I T E ( I W 2 ~ 1 1 0 ) T R A V ~ F I C H ( lqFICH(2)
)
rFICH(4)
.FICH(8)
FORMAT(
'DKEND
TRAVAIL 'rAlr' SUR FICHIER
WRITE(IW2rl16) JE
FORV4T (
13)
DO 16 J=2*JE
JMl=J-1
WRITE(IW2,117) (DI(J*K).K=lrJYl)
FOWAT(
8F9.4)
GOTO 1001
MATRICE DES DONNEES
LE FICHIER Du t 9
rFICH(5)
',6Al)
rFICH(71,
16
1002
120
600
700
601
701
602
702
800
801
1000
2000
WRITE(IWRITI~~O)VIFICH
PA5DANS
F O R M A T ( ~ H ~ Y ~ ~ ANE
~ , SE
/ TROUVE
@
GOTO 1 0 0 0
WRITE(IWRIT,700)
FORMAT(lH1,'JE
HORS DES L I M I T E S ' )
GOTO 1 0 0 0
WRITE(IWRIT,701)
FORMAT(lH1,'IVM
HORS D E SL I M I T E S ' )
GOTO 1 0 0 0
WRITE(IWRITI~O~)
HORS D E SL I M I T E S ' )
FORMAT(lH1,'IV
GOTO 1 0 0 0
WRITE(I~RIT~BO1)
F O R M A T ( l H l r ' P L U . 5 DE 1 0 CARTESFORMAT')
STOP
END
SUBROUTINELECT(IVM,JEM*UIVD,B)
COMMON/C1/IREAD~IdRIT,IPUIKFOR
I M P L I C I TI N T E G E R ( V )
DIMENSION U(1) 'VD(1)
D I M E N S I O N 7(400)r F M T ( 1 3 0 )
KFOR=KFOR*13
READ(IREAD*2000) (FMT(I),I=lrKFOR)
FORMAT (13A6)
L=O
DO SOO J = l . J E M
READ(IREADtFMT,END=502) (Z(I)rI=l~IVM)
DO 5 0 1 I = l r I V M
501
I F ( V D ( 1 ) mNE.1)GOTO
L=L+1
501
500
502
1
U(L)=Z(I)
CONTINUE
CONTINUE
RETURN
WRITE(IWRIT*l)
F O R M A T ( @ I L MANQUE
RETURN 5
END
DES CARTESDONNEESO
SURROUTINERGEQU(XJ)
COMMON/C4/JE,I,JRE,IV
COMMOPJ X ( 3 0 0 0 ) r X B I S ( 4 2 0 0 0 )
D I M E N S I O NA ( 1 0 0 ) q R ( 1 0 0 )
XRIS(l)=O.
T=O.
DO 1 J = l r J E
H(J)=O.
DO 2 J = l r J E
1
Jl=J-1
2
3
A(J)=X(Jl*IV+I)
MP= J E
DO 13 J=l * J E
IF(R(J))3*3*13
IF(R(J)+1.)4r13*4
4
s=o.
6
E=O e
30
31
Y=A ( J I
IF(Y+1.)30912,30
DO 9 K = l , J E
IF(A(K)+1.)31,9131
If ( A ( K ) - Y ) 7 r 8 r 9
LE F I C H I E R ' r 8 A l )
17
7
8
9
32
10
11
12
13
17
18
22
11
12
1
2
3
4
5
R
IF ( D ) 2 r f J 9 7
I=35
GOTO 6
7
1= 3 4
6
IP=36-1
C (K)=CO ( IP)
RETURN
END
10
18
H M V T : CONSTELLATION§
(La premièreversion de ceprogrammea &té écrite par M.-C. MASBOU du C.N.R.S., en 1965)
. _.
.
-. . --..
.
. - - ~ .-. - . . ---_____ .-
__
COMMON /C1/IM(100~100)rN(100)
*/C~/IMIN*IINXIIREADIIWRIT
!
650
30
*/C2/IE(lOO)~IJ(1001
COMMON C(101r400)~0(100~10013)
DIMENSION 2(6)rFMT(130)
DATA AL/'
'/
DIMENSION NZEHO(17)
DIMENSION Sll(100)
DO 6 5 0 I=lr100
IE(I)=I
DO 30 1 ~ 1 9 1 7
NZERO (11 = O
READ(5r300)IREAD,IWRIT,IW2
FORMAT (312)
8
READ(IREAD~105,END=200) (2!(1)91=1e6)
105 FORMAT (6A6)
300
READ(IREAD~100)NG~NXtKFORIKIMP
100
FORMAT(I39312)
IF(KFOReGT.1O)GOTO 5 0 0
DO A A I=l,NG
N(I)=O
DO 8 7 J=l,NG '
87
IM(I,J)=O
DO 88J=lr400
C(I,J)=O
8 8 CONTINUE
K=l
NGl=NG*l
IF (NX.EQ.1) GOTO 4
DO11=2rNG
IN=I-l
READ(IREAD9101)
1
[D(IqJ,K),J=l,IN)
101 FORMAT(BF9.4)
NG2=NG-1
DO21=1,NG2
D(I,I,l)=O
IN=I*l
DO 2 JzINING
2 D(I,JIK)=D(J,I*K)
GOTO 17
4
READ(IREAD9103)NV
1 0 3 FORMAT (13)
DO 400 L=1,130
400 FMT (LI=BL
IF(KFOR.EO.O)KFOR=l
KFOR=KFOR*13
READ(IPEAD,600)(FMT(I),I=1,KFOR)
600
I
FORMAT(13Ah)
DO 5 1=19NG
READ ( IREADqFMT) ( C ( 19 J) 9 J=l ,NV)
5 CONTINUE
READ(IREADIFMT) (C(NGl,J),J=lrNV)
NIZNG-1
D061=1 ,NI
D(IIII~)=O
IL=I+1
DO 6JzILvNG
19
NVV=O
s=o
7 P = l rNV
IF(C(IrM).LT.O..OR.C(J,M).LT.O.)
GOTO 7
S= ( (C( 1 rM) -C (J9M)) / C (NGlr M ) 1 * * 2 + S
NVV=!J\(V+l
7 CONTINUE
IF(NVV)9,10,9
10
D(IrJ*l)=ll.
GOTO 6
9
CONTINUE
D~I~Jr1~~~S~lOOOO~/NVV
6 C~(J*Irl)=D(I~J*l)
17
IF(KIMP.EQ.O)GOTO
177
CALL PANG(MGr2)
177
WRITE(IWRITIZ~)NF
26
FORPAT(lHl*'NOMBRE DE GROUPES AVANT COMPACTAGE'rI4)
690
DO 700 I = l * Y G
Sll (I)=O
DO 700 J=l,NG
C!c!
IF~~~I1J~1~~11~~700~7011700
701
700
711
710
Y01
720
730
900
910
920
21
20
25
24
22
23
500
501
ZOO
SIl(I)=Sll(I)+l
CONTINUE
SllM=O
O0 710 I=l*NG
IF(S11 (I)-SllM)710r710r711
SllFO.=S11(1)
151 1M=I
CONTINUE
IF(S11M)901r900r901
DO 7 2 0 I=lrNG
D(IrISllMrI)=D(NGrIrl)
DO 730 I=lrNG
D~IS11~rIrI~~D~NG~Irl~
I=IE(IUG)
1E (NG)=IE (IS11M)
IE(ISllM)=I
NGZNG-1
60TO 6 9 0
WRITE(IWQIT,~~O)NG
FORtMAT ( 1 H 9 ' NOUVEAU NOMBRE DE GROUPES
CALL STAR ( 1 , N C r Z )
IF(IMIN.EO.O)GOTO
8
O0 20 K = l r IMIlJ
MK=N ( K )
IF(MK.GT.17)GOTO
22
DO 9 2 0 J=l,%K
L=IM (Kcr J)
IJ(J)=IE(L)
WRITE(IW2,Zl) (IJ(J)vJ=l,MK)
9
rI4)
FORM4T(17I3)
CONTINUE
IF(I~IIN.EO.I)GOTO 24
INIkPl=IMIN+l
0 0 2 5 K=IMINPltI
L=IMIK*1)
Jl=IE (L)
WHITE(IW2r?.l)Jl
UI?ITE(IW~~~~)NZERO
GOTO 8
WHITE(6123)K
ÇORMAT(1Hlr'PLUS DE 17 GROUPES DANS LA CnNSTELLATION'~I4)
WHITE(I~!QITI~O~)
FORMAT(' PLUS DE 1 0 CARTES FORMAT')
STOP
END
20
21
NI=NG-1
DMIN=i7(lr2,K)
IMIN=l
JI41 NE,?
DO 1 I=l*NI
IL=I+1
DO 1 J=IL.NG
IF (D(I*J*Y)-DMAX)4r4*2
2 DMAX=D(I*J*K)
IMAXZI
4
3
1
7
C
200
JMAX=J
GO TO 1
IF(D(I*J,K)-DMIN)3.1tl
DMIN=D(IIJ,K)
IMIN=I
JMII'l=J
CONTINUE.
IF((DMeX-OMIN).LT.DMIN)GO
TO 6
1=1
CIUlAX=DYIN
CM(I)=DMIN
SVIN=PYIN
IM(I?l)=IMIN
IM(II~)=JYIN
N(I)=2
NC=2
CALCUL DE D2 EXTHP CONSTELLATION *GRAND QUE CYPX
IF(I.EQ.1)GOTO
100
C(ItJJ)=C(I1JJ)+D(KI,KLIto
CONTINUE
C(I,JJ)=C(I*JJ)/(ILOIN)
C(JJ~I)=C(IIJJI
1 7 CONTINUE
DO 2 8 IJ=lrJI
JA=IJ+l
DO 2R JJ=JA*I
27
IF(C(IJIJJ).LE.CMAX.OR.C(IJIJ~JJ).LE.~~IN)
28 CONTINUE
100 6=100000
BE=" ( 1 1
C Y ( 1) =DMIN
RECI?EWHE D UN ACCROISSFYENT MINIMUM
C
DO 8 W l r N C
DO 4 9 L=lrI
NL=N (LI
DO 4 9 J=lrNL
IF(M.EQmIM(L*J) )GOTO 8
4 9 CONTINUE
NJ=N ( 1 )
A=O
J=lrElJ
NK=IM(I*J)
A=D(MqNKrK) + A
48 CONTINUE
IF(A.GT.R)COTO
DO 48
8
GO TO 130
22
IM(I,NJ+l)=M
tl=A
8 CONT 1NUE
CDISTANCEYOYENNE
INTRA CONSTELLATION
IF(R-100000)132~129r132
132
NO=1
O 0 47 L=2rNJ
47 ND=ND+L
SMIN=SMIN+R
DMIN=SMIN/ND
NJ=NJ+1
C
COt4PARAISON ENTRE 0 2 INTRAEXTRACONSTELL4TION ET INTPe CONSTELL4TIO
N (1 1 =NJ
DO 13 ."(=lrNG
DO 14 J=l r I
NI=N(J)
DO 14 L=lrNI
IF(IM(JrL)-M)14r13rl4
14 CONTINUE
DO 23 L=l,NJ
IK=IM ( 1 rL)
13 CONTINUE
IF(CM(I).GT.CMAX)CMAX=CM(I)
GOTO 200
130 DO131 J=lrJI
C(IqJ)=Y(J)
131 C(J*I)=C(I*Jl
128 N(I)=N(I)-l
IF(N(I).LE.l)GOTO
1000
129 DYIN=lOOOOO
NCAS=l
IL=1
M2=1
20 DO 81 M=IL*NG
DO 90 L=l*I
IN=N(L)
DO 90 NU=lrIN
IF(M.EQ.IM(LrNU))GOTO
81
Y0 CONTINUE
GOTO(a?r43,44)rNCAS
42 NCAS=2
M l =M
GOTO 81
.
43 NCAS=3
M2=M
44 I F (DM1N.LT.D (Ml r M , K ) ) G O 1 0 81
DMIN=D ( M l r M r K )
hl
.
IMIN=Ml.
JMIN=M
CONTIFjUE
IF~PYIN-10~000~133r1O~br133
133
IF(Ml.EO.M2)GOTO
IL=@
39
23
klRITE(1wR1T1109) (Z(Y)rM=lr6)
KMIN=IMIN/15+0~001
IRMIN=IMIN-15*KMIN+O.001
IF(K,MIN.EO.O)GOTO 300
DO 301 KMZlrKMIN
11=1+15*(tiY-1)+0.001
I~AX=lS*KM+O.OGl
JMAX=IMAX
CALL SPl(Ilr1MAXrJMAX)
CALL SP3 ( 11 r JMAX)
301
CONTINUE
300
IF(IHt-~IN.EQ.O)GOTO 2 0 0 0
I1=15*K~IN+l+0.001
IMAX=IMIN
CALL SPl(I1rIMAXrIMAX)
CALL SP3(11rIMAX)
GOTO 2000
6
WRITE(IdRITr108)
108 FORPAT(~OXIJ~HILY AN PAS DE CONSTELLATIONS
2 0 0 0 RETURN
END
103
10
104
50
)
24
*
O
O
IF(NX.EO.1)GO
DO 20 K = l r 5
FMTl (K)=FTI (K
FMT2 ( K I =FT2 (K
GOTO D O
20
10
40
30
RO
bh
11O
CONTINUE
~HITE(IWRIT~110)
FORVAT (1H1)
RFTUKN
END
75
MODE D'EMPLOI DES PROGRAMMES D K A T 5 .HM V T
CARTES CONTROLE
1
- CARTESCONTROLESYSTEMEETPROGRAMMEPOURL'OBTENTIONDESDISTANCES
DE RANGS,PUISREGROUPEMENTENCONSTELLATIONSDE
m UNITES.
26
-
2
CARTESCONTROLESYSTEME ET PROGRAMMEPOURL'OBTENTIONDESDISTANCES
DE RANGS ENTRE m UNITES.
"RUNIG
" ASG
ORDKAT~TP----~ORSTOM
A. 9 R C r 395i)L
*T
"MOVE
A.
r2
"FTNDrA 4.DKAT
"COPINIA A - O K A T r
"FREE A .
" A S G q T 10. r F 2
" A S G r T R.rFZ
"XOT
.DKAT
050608
Fichierd'entreestandard
"DATAIL
"END
(rn, v, v')
2 < rn < 1O0
1 G v'
< v < 400
8.
"F 1 N
3
-
CARTESCONTROLESYSTEME ET PROGRAMMEPOUR LE REGROUPEMENT EN
CONSTELLATIONS DE m UNITES,APARTIR DE DISTANCESQUELCONQUESENTRE
CES m UNITES.
"RUNrGORYMVTrTP---rORSTOM
"ASGIT
A.rbCq3952L
"NOVE e . * 2
" F I N O r A A.HMVT
"COPIN,A A - H M V T r
"FWE
A.
" A S G r T 4. rF2
llXOT .HMVT
050604
FICHIER 1 2 - 0 6 - 4 3
"'1"'
rnGID0 en 13
Moitiéinférieure de la matrice (rn x rn) dedistancesentre
prises 2 à 2,au format (8F9.4).
Lm
"DATAvL
"EWD
"F
1N
4.
les rn unités,
' 27
FORME D U N FICHIER DENTREE STANDARD : EXEMPLE
NUMERO
DE CARTE
El
Date &identification du fichier de v variables
relevéessur m unités,dans les colonnes 9 à 16.
a
kfor en colonnes 10-11 = nombre decartes format
de lecture desdonnées (l<kfor<lO), format
(12),
kfor = 1 pardéfaut. 1 = numéro du travail sur
le fichier précité <9, en colonne 9.
v : nombre total des variables du fichier, (141,
en colonnes 4-7.
m : nombre d'unités (ou groupes), (13),
en colonnes1-3.
El
El
m : cartes libellés desunités, 48 caractères
maximum, en colonnes 1-48.
Jm-t3)
a
v : cartes libellés desvariables,
maximum, en colonnes 1-48.
)m+v+4]
4
48 caractères
v' : nombre devariables à traiter parmi les v. (14),
en colonnes 1-4. Si v'<v, cettecarte doit êtresuivie
desv'noms
(ou libellés) des variables à retenir,
dans l'ordre.
kfor cartek) format de lecture d'unenregistrement
de v données,en colonnes 1-98.
F
Donnéesde m enregistrements de v variables à
lire au format précité. Si certainesdonnées
manquent, elles sont à remplacerpardes -1.
28
EXEMPLE DE SORTIE
L‘utilisation conjointe des programmes D K A T et H M V T donne les sorties d’imprimante que voici.
@COPIN
A.9
6 ABS
@Jl)ELETE*C Hm
FURPUH 0 0 % 6 - l l / l R - 1 2 : 4 8
R
IS N O 1 CATALOGUED OR ASSIGNE0
F A CS T A T U SG
: 00010000000
W X Q T .DI<AT
y:
n
a.
...
.
.
29
t t t
O-N
-
0
OO
vv
3
30
31
IdDATAqL 8.
DATA T 7 RL70-S 11/18-12:48:52
1.
DKEND
TWAV41L
1 SUR F I C H I E H 120643
2.
3.
4.
5.
6.
7.
8.
END
'
7
-1630
-5043
a1903
a0597
.3610
-6354
.5062
el732
,2342
*O520
-2177
-2117
.PO69
-6939
.5264
e3395 .4019
a3261
-7422
,0726
DATA.
Moitié inférieure de la matrice des distances, non plusreprésentées par leur classe d'appartenance,
mais par un nombre décimai comportant 4 chiffres après la virgule.
Les distances incalculables sont notées 11.
A titre d'exemple, lerésultat.O726
de la matrice desdistances).
est la distance entrel'unité7(ligne
I W A T A t N 8.99.
DATA T 7 RL70-5 11/18-12:48:53
ENI) DATA.
IMAGE
COUNT:
A
@ERS 8.
FllRPUP 00t6-11/1~-12:49
7) et l'unité6(colonne
6
d
d
"""0WW
......
(Y
0.
3
O
r:
VI
c
ü
I
C
c
6
32
4
33
l
I
I
I
I
l
I
I
I
l
I
l
l
I
l
I
l
I
l
I
l
I
I
I
I
I
I
I
I
I
I
l
I
I
I
I
I
l
l
I
I
I
I
I
I
I
I
I
I
I
I
I
I
l
I
l
l
I
I
I
l-
t
1
I
I
I
I
l
I
I
I
cl
w
U
34
LSOATAIL 8.
DATA T7 WL70-5~11/18-12:48:58
1.
1 7 6
2.
2 4 5 3
3.
0
0
0
0
0
0
0
0
0
0
0
0
END DATA.
Q
LFREE 8.
reEPS 0.
FURPUR 0026-11/30-12:48
I F 1N
0
0
0
0
0
35
PROGRAMME C S AT POUR CORRÉLATIONS
36
1
20
200
50 1
aooo
502
2001
20 1
90
722
81
a2
83
61
62
KM=O
QEAP(IREAD*l)JEqIZ
FORMAT (139 14)
IF(JE.GT.lP7.OR.JE.LT.l)GOTCJ
5002
IF(IZ.i_T.?.OR.I/.GT.345)GOTO
5001
RFAD(IREADt2) ((Z(KrI)rK=lq9).1=lrIZ)
FORMAT(lX19A6)
PEAD(IQEADq200)IV
FOQMAT (14)
IF(IV.LT.?.OR.IV.GT.~&~)~~T~
7000
IF(IV.YE.IZ)GOTO
ROOO
DO 5 0 1 J = l r I V
O0 5 0 1 1 ~ 1 - 9
Y(ItJ)=Z(I*J)
GOTO 5 0 2
READ(IREADt2) ((Y(ItJ)rI=lr9)rJ=lrIV)
READ(IREADt201)(FMT(I)rI=l*~FOR)
FORMAT(lXr13A6)
FORMAT ( 1 3 A h1
DO 90 J = l r J E
READ(1READtFMT) ( X ( I q J ) t I = l r I Z )
KM= 1
JQO=l
1=0
JP=O
J Q = J O*0J P
WRITE(IWRITt73)A4r451JQIA4*A5rJQ
I=I+l
WRITE(IWRITt511 ( Z ( K t l ) r K = l r 9 ) t ( Z ( K ~ T ) r K = l r 9 )
IF(I.EQ.IZ)GOTO
a3
JP= 1/44
IF ( ( I - J P o 4 4 ) eEQ.O)GOTO A l
GOTO 82
I=O
JP=O
JQ=JQ*l
WRITE(IWRITt60)A4rA5~Al~JQ~A4rA5rA1,JQ
WRITE(IWRITt50)
I=I*l
WRITE(IWRITt51) (Y (KII) qK=l*9) t ( Y ( K q I ) rY=lrcj)
IF(I.EQ.IV)GOTO
63
JP=I /44
IF((I-JP*44).EQ.O)GOTObl
63
GOTO 6 2
M=O
1M=O
N=O
IK=O
J=O
IVMl=IV-l
DO 1 0 0 I = l r I V M l
IPl=I+l
DO 1 0 0 K = I P l r I V
5 04
503
SOS
5000
J=J*l
IF(IV.NE.IZ)GOTO
503
DO 5 0 4 I J = I t J E
A(IJ)=X(IrIJ)
B(IJ)=X(KIIJ)
GOTO 505
C A L LE X T R ( I t K t A 9 H )
32
IF(IB.EO.0)GOTO
C A L LE F C O M ( A q 8 r L )
IF(L.GT.5.AND.L.LT.l27)60T0
1TAR=-6
M=M+ 1
WRITE (3'Jt500) ITABtL
103
37
GOTO 1 0 0
CALL S R A N K ( A I H I R I L ~ R S ~ $ ~ O O O )
CALL T E S T ( P S I R C S ~ R C ~ ~ L ~ I T A R )
WRITE(3'J~SOO)ITAB,L
FQ4MAT(IZ!rI3)
CONTINUE
IF (M.EQ.O)GOTO
13
IV2=IV*(IV-l)
IV2=(IV2+1)/2
J=O
JR=O
JP=O
JQ=
1JQ+
103
500
100
41
WRITE(IWRIT,60)A4,AS,AlrJOIA4rASI91,JO
42
13
43
WHITE(IWRIT,52)
J=J+ 1
HEAD(3'J*500)IT4@rL
IF(ITAR.NE.-6)GOTO
42
CALL EOUIV(JII,K)
W R I T E ~ I W R I T ~ ~ ~ ~ ~ Y ~ I J , I ) ~ I J ~ ~ I ~ ~ , ( Y ~ I J ,
S(Y(IJIK),IJ=~,~
),L,L
JR=JR+I
IF(JR.EQ.M
)GOTO 13
JP=JR/ 11
IF((JR-JP*ll)eEQ.O)COTO
41
GOTQ 42
IF(N.EQ.O)GOTO
30
J=O
JH=O
JP=O
JQ=JQ+1
dRITE(IWRIT,60)A4,A5,A1,JWIAL,AS,AlrJQ
44
30
49
NRITE(IWRITq53)
J=J+1
READ(J'JISOO)ITARIL
IF(ITPB.NE.O)GOTO
44
CALL EQUIV(J~IIK)
WRITE(IWRIT,56) ~ Y ~ I J ~ I ~ ~ I J ~ l r 9 ~ ~ ~ Y ~ I J ~ I ~ . I J = l r 9 ~ r ( Y
B(Y(IJ,K) tIJ=199 1 ,L*L
JR=JR+l
IF (JR.EQ.N )GOTO 30
JP=JR/ll
IF((JR-JP*ll).EQ.O)GOTO
43
GOTO 44
IF(IK.EQ.O)GOTO
31
J=O
JR=O
JP=O
JO=JQ+l
WRITE(IW41T,60)A4,A51A11JQIA4,AS,AlrJO
46
YRITE(IWRIT~54)
J=J+1
REA@(3'J~5001ITAB,L
IF(ITAR.LE.O)GOTO
46
CALL EQUIV(J,I,Y) .
dRITE(IWRIT,57) ~ Y ~ I J I I ~ ~ I J ~ ~ ~ ~ ~ I ~ Y ~ I J I I ) ,
A(Y(IJqK),IJ=lr9 )IL,ITAR~LIITAH
JR=JR+l
31
IF(JR.EQ.IK)GOTO
JP=JR/11
IF((JR-JP*ll).EQ.O)GOTO
49
GQTQ 46
IF(JM.EQ.O)GOTO .9999
'
31
J=O
.
.
'
,.
JR=O
JP=O
38
47
JO=JO+l
WWITE(IWRiTr60)A4rA5rAl~JOrA4.ASrAl~JQ
hn
50
51
52
53
NPITE(IWRIT~SS)
J=J+1
HFAD(3*J*SnO)ITAStL
IF(ITAS.NE.-5.4ND.ITAH.NE.-L)GOTO
48
IT4R=-ITAR
C A L LE O U I V ( J * I , K )
~ R I T € ( I W Q I T I S ~~) Y ~ I J ~ I ~ ~ I J ~ 1 ~ 9 ~ , ~ Y ~ I J ~ I ~ r I J ~ l ~ 9 ~ r (
I ( Y ( I J , K ) , I J = 1 * 9) r L t I T A H , L t I ' T A 8
JR= J R + 1
IF(JR.EQ.IM)GOTO
9999
JP=JR/11
IF((JR-JP~ll).EO.O)GOTO
47
bOTO 48
F O Y M A T ( l X ~ * C O R R E L A T I O NDE
SPEAkMAN
ENTRE"'r32Xr
u
*COWRELATIOY DE S P F A
ER
N kTARNE " * / / )
F O H M A T ( l X 9* A h r R X 9 9 A 4 )
FORMAT(~XI'CORRELATIONDESPEARMANINCALCULABLEPOUR"',20Xr
0
' C 9 H R E L A T I O N DE SPEARMAN I f i C A L C U LPAO
BU
LE
R"*//)
F O R M A T ( l X * ' A R S E N C E DE CORRELATION DE SPEARMAN POUR L E SV A R I A I 3 L E S " Q
Q*HX*
'AHSENCE D t CORaELATTGN DE SPFAKMAN POUR L EVSA R I A H L t S " *
"//)
54
55
F O H M A T ( l X * ' C O R R E L 4 T I O N DE SPFARMAN P O S I T I V E ENTRELESVAR1AbLES"'r
"9X *
'CORRELATION DE SPEAHMAN P O S I T I VEEN T RL EVSA R I A B L E S l l t /
"1)
F O R M A T ( 1 X v ' C O R H E L A T I O N DESPEAHMAhNEGATIVE
ENTFlE L F SV A R I A B L E S " ' *
"9X *
*CORRELATION DE SPEAHMAN
NEGATIVE
ENTRE
LES
VARIAHLES"'/
"/)
56
57
60
2
3
7
73
7000
hono
5001
6001
5002
6002
5500
550 1
32
F O R M A T ( ~ X I ~ A ~ ~ R X , ~ A S / ~ /Xl X ~* * E
~ FAF' E
r ~IC3 T
, 5I 0F X , ' E F
Q F E C T I F '*13/)
FORMATllX*9A6,9X*9A6/1~*9AS*~X*~A6
/ l X * ' E F F E C T I F* r I 3 , 1 0 X , * R I
*SQOE'.I3r'
POUR C E N T ' , 2 1 X * ' E F F E C T I F* t I 3 r l O X r * R I S Q U E ' r I 3 9 ,
POUR CE
"NT'/)
F O R M A T ( ~ H ~ ~ ~ ~ X I * F I C H I E R' e Z A 4 9 T' R A V A I L* 9 A l * 2 X , * P A G E' r 1 3 1 2 6 X * * F I
/ C H I E R * r 2 A 4 . * T R A V A I *L* A l , Z X , ' P A G E
'913//)
FORMAT (9Ah)
FORMAT(~X,?A~/HXIA~.I~)
FORMAT ( 8 4 )
~ O R P A T ( l H l v 2 Z X ~ ' F I C H I E K' r ? A 4 , 1 2 X * ' P P G F* , 1 3 * 2 6 X * ' F I C H I E P
**2A4r
"~ZXQ'PAGF '*13//)
dRITE(IWRIT~6000)
F O R b ' A T ( 1 H l r ' L E NOMBRE DE V A R I A B L E S A CORQELE2N*'ESTPASDANSLES
2-346')
. L I Y I T E SP E R M I S E S "
GOTO 32
WRITE(IWHIT~6001)
DANS Lb.S
F O R M & T ( l H l , * L E NOMRRE DE V A R I A B L E S OLJ F I C H I E RN * * E S TP A S
. L I M I T F S PEPE(1SES" 2 - 3 4 6 ' )
GOTO 32
WRITF(Ik/RIT,6002)
FORMAT(lHl.'LE
NOMBRE D " U N 1 T E SN * ' F S TP A SB A N SL E SL I M I T E SP E H Y I S
.ES'' 1-1 27
)
bQTO 32
WRITE(IhRITt5501)
F O R M P T ( l H l * ' I L Y A PLUS DE 1 0 CAGTES F O R M A T 0
STOP
END
-'
1=L
K=J+LP-(L-l)*IV
39
GOTO 2
CONTINUE
RETURN
END
SURROUTINEEFCOM(4r4,L)
DIMENSION A ( l ) r R ( l )
COMMON /C1/JEM
INTEGERQ2 L
1
2
L=O
DO 1 J = l * J E M
IF(A(J)+1)3t1,3
IF(P(J)+1.)2tlr2
L=L+1
A(L)=A(J)
R(L)=R(J)
CONTINUE
HETURN
END
SUBPOUTINE L O C ( I * I L )
COMMObJ/C6/Y(9r346)rZ(9,346)/C7/IZ1IVtIPtIWRIT
INTFGER Y t Z
IH=1
DO 1 K = l r I Z
DO Z J=199
IF(Z(J,K).NE.Y(JtI))GOTO
1
CONTINUE
1L=K
GOTO 3
CONTINUE
2
1
WRITE(I~RITtS)(Y(J,I)1J=1,9)r(Y(JII),J=1t9)
IR=O
FORMAT(lHlr9A6*8X~9A6
5
/'
Q
N ' @ A PAS € THEE C O N N U E ' t 4 2 % ,
* b i t ' ARECONNUE*
P A S ETE
RETURN
END
3
SUBWOUTINE E X T K ( I * K , A * B )
D I M E N S I O NP ( l ) , B ( l I
COMMON
/C7/IZtIV,IR /Cl/JE
COYWON X ( 3 4 6 r 1 2 7 )
C A L LL O C ( I 9 I L )
C A L LL O C ( K 9 K L )
IF(IR.EQ.O)GOTO
2
DO 1 J = l r J E
A(J)=X(ILtJ)
1
2
END
1
S U H Q O U T I N ET E S T ( R S , R 5 t W l r L , I T A B )
D I M E N S I O N RE; ( 11 qR1 ( 1 )
COMMON / C 2 / N , I M c , I K ~ J
I F (ARS ( R S ) .GT.RS ( L ) IGOTO 1
ITAR=O
&=N+ 1
GOTO 1 0
I F (ABS(PS).GT.H1 ( L ) )GOTO 3
IF(RS.GT.O)GOTO
2
.
ITAR=-S
,
IM=IM+l
. .
' GOTO 10
.
.
1TAQ=+5
IK=IK+l
GOTO 1 0
.
.
I F (RS.GT.O)GOTO
4
'
2
:
)
40
ITAa=-1
IM=IM+1
4
1O
GOTO 1 0
ITA\H=+l
IK=IK+l
RETllRN
EN0
-(Voir note infra paginale)
5
40
50
55
57
60
2
3
70
1
KETUHN 6
END
SURROUTIME T I E ( R * N * K T v T )
C
DIMENSION R (1)
T=O.
Y=O.
5
10
20
30
40
50
60
70
75
00
90
X=l.E3R
IND=O
00 30 I = l r N
IF(R(I)-Y)30,30*10
IF(R(I)-X)?Ot30*30
X=H(I)
IND=IND+l
CONTINUE
IF(IND)90*90,40
Y=x
CT=O.
DO 6 0 I = l r N
IF(R(I)-X)60*5Gr60
CT=CT+I
CONT 1 NUE
I F ( C T ) 7095.70
IF (KT-1) 75*80*7S
T=T+CT* (CT-1.) 1 2 .
GOTO 5
T = T + (CT*CT*CT-C7)/12.
GOTO 5
HETURN
..
EMI
* D’aprèsles(wbroutinesa SRANK. TIE, RANK dans IBM.1966. System1360ScientificSubroutinePackage
(360 A - CM - 03X) Programmer’sManual. IBM ApplicationProgramn’H20-0205-0.WhitePlains,
157 p.
41
SUAROUTINE RANK (AIRIN)
DIMFNSION A(l),R(l)
00 10 I=lrN
10
R(I)=0.0
DO 100 I = l * N
.?O
IF(R(I))%O*20,100
SVALL=O .O
EQUAL=O .O
-
-
X=4(I)
DO 50 J=lrbl
30
40
50
60
70
80
90
1 O0
IF (A(J)-X)30r40r50
SMALL=SMALL*l.
GOTO 50
EQUAL=EQUAL+I.
R(J)=-l.
CONTINIJE
IF(EOUAL-l.)60r6Or70
R(I)=SMALL+l.
5 0 T O 100
P=SMALL+EQUAL"(E~UAL*l.)/(E~lJAL+EQUL)
D O 90 JZlrN
IF (R(J)+1.)90*80*90
R(J)=P
CONTINUE
CONTINUE
RETURN
END
43
MODED'EMPLOI
'
CARTES CONTROLE
'
'
DU PROGRAMME C S A T
'
. .
Les cartes contrôle que voici permettent de calculer et de tester les coefficients de corrélation de
SPEARMAN (1904) entre v' variables prises 2 B 2.
"RUN,J/T
OKCSATtTP___tORSTOM._
"AS6 9 T A . t R C t 3953L
"MAOVE A. 9 2
"FINDtA A.CSAT
"COPIN.4
A-CSATt
"FREE A .
"XOT
.CS4T
0506
1
Fichier d'entréestandard moins les noms d'unités
EXEMPLE D'ENTREE
6Gm
< 127
1 < v'
< v G 346
44
EXEMPLE DE SORTIE
F I C H I E R 12-06-43
PAGE
1
1
PAGE
2
TPAVAIL 1
PAGE
3
CAEACTFRE 1
Liste des v variables présentes dans le fichier.
CAPACTFRE 3
F I C H I E R 12-06-43
TRAVAIL
CORRELATIONDESPEARMANENTRE:
CAHACTERE 1
\
F I C H I E R 12-06-43
ARSENCEDECORRELATIOY
DE SPFARt4P.N POUR LES VAPIAQLES:
C4RACTERE: 1
CAPACTFRE 7
'
CP'ACTERF:
ChPACTFRE 3
EFFECTIF
7
CAPACTERE 2
ChPACTERE 3
EFFECTIF
7
Résultats des .
:
C tests de corrélation.
Dans l'ordre,
quatre
rubriques
peuvent
apparaïtre
:
.- corrélations incalculables (effectif trop faible pour le couple
concerné);
- absences de corrélation, l'effectif est indiqué;'
- corrélations positives, l'effectif e t le risque (5 % ou 1 %) sontindiqués;
- corrélations négatives, l'effectif e t le risque (5 76 ou 1 X)sontindiqués.
La sortie est dédoublée verticalement.
45
PROGRAMME M E D T POUR MÉDIANES
D I M E N S I O N F I C H ( 2 ) 9 E C H ( 8 7 1 2 7 ) q V A K ( 9 ) r V ( 9 ) . X ( l ? 7 ) r X f l ( l Z 7 ) , I l t 1 ) ( 1 7 7 )*
"FYT(130) rVD(300) ,2(30O) ,Vi7(9.300)
COMMON K F O R q Y ( 3 8 1 0 0 )
I M P L I C I TI N T E G E H
1V)
D E F I N EF I L E2 ( 3 0 0 , 9 9 U , I D )
C A L LP R T M ( h h r O t 0 )
ID=O
hl
soo
100
101
102
50
205
103
51
52
104
READ(Sr61)N,NW
FQRMAT(212)
KEAD(N~100~END=1000)FICH
FOHMAT(HX,?Ah)
HEAD(N,lOl)TRAV,KFOR
FORMAT(8X,A1912)
IF(KFOR.GT.1O)GOTO5500
IF(KFOR.EQ.O)KFOR=l
KFOR=KFOPQI 3
READ(NvlO2)JE,IVM
FORMAT(I3rI4)
I F ( J E " I V M ~ L T ~ 2 ~ 0 R ~ J E Q I V M . G T ~ 3 R l O O3~0~0O
0TO
IF(JE.LT.2.OR.JE.GT.127)GOTO
4000
IF(IVM.LT.l.OR.IVM.GT.2000)GOTO
5000
DO 5 0 I = l r J E
P E A D ( N t 2 0 5 ) ( E C H ( J , I ) 9J=196)
FORMAT(8A61
FORPAT ( 9 A 6 )
0 0 51 I=l
,IVM
READ(N9103)VAR
WHITE(2'I)VAR
0 0 52 I = l , I V M
VD(I)=l.
READ(Nr104)IV
IF(IV.GT.IVM)GOTO
6000
FORMAT(I4)
I F ( I V . E O . I V M ) G O T O5 3 9
I=O
D O 54 L=19IV
READ ( K r 1 0 3 ) V
546
56
57
IF(I.FO.IVV)GOTO
7000
0 0 546 J = l r 9
VR(J,L)=V(J)
I=I+1
RE4l)(2*I)VAR
DO 57 K ~ 1 . 9
IF(VAR(K).NE.V(K))GOTO
CONTINUE
GOTO 5 4
58
46
58
54
59
539
5hO
53
105
501
502
204
Q
DO 10000
40
2
1
3
4
5
6
9
10
11
12
7
'FICHIER 'r2P6r'TRAVAIL
IVAQ=lrIV
'*Al
//lX*4Phr14X19fi6///)
47
YE=O
MEl=O
ME2=0
MED 1= O
MEDZ=O
MED=O
IQ=(~+O.OOl)/Z+0.001
IF(M-IQ*2)14,13,14
13
ME1=MQ+(M+~.001)/2+0.001
MEZ=ME 1 +1
VEDI=II\ID "€1)
t4€DZ=IND(YE2)
IN=MEl-IR
IS=ME2+IR
CJOTO 15
14
15
ME=~~+(M+1.001)/2+0.001
MED=IND (ME)
1N=ME- 1 R
IS=KE+IR
INF=IND(IN)
ISUP=IYD(IS)
INMlzIN-1
WHITE(~*~~~)FICHITRAVIFICH,TRAVI(VR(JIIVA~),J=~*~)*(V~(J*IVAR)*J=
*99)
16 I=19INMl
INDI=IND(X)
wRITE(6920O)X(I) 9 (ECH(J*Ih'DI) tJ=ltt))
9X(I) 9 (tCH(J9INDI) rJ=1+8)
WRITE(h*201)X(IIJ)9 (ECH(J9INF) ,J=l,L3) ,X(IN) 7 (ECH(J,INF),J=1,8)
DO
16
I F ( Y E ) 18917918
17
18
19
MEr"il=MEl-l
CiOTO 19
I.FEMl=ME-l
INPl=IN+l
DO 20 I=INPlrMEMl
I!JDI=IND ( 1 1
20
21
22
23
24
25
INDI=IND(I)
WRITE(6rEOO)X(I) ,(ECH(J,INDI) 9J=198),X(I) ,(LC4(J,INDI) ,J=lr8)
NRITE(NW1999) ( V R ( J , I W A R ) , J = 1 ~ 9 ) , M , X o , X ( I N ) ~ X ( I S ) , F I C ~
999
F@RMAT(~A~,I~YZF~O.~,~A~)
lonoo
CONTINUE
500
5500WRITE(6*5501)
5501' FORMPT(~HZI'PLUS DE 10 CARTES FORMAT
ANNONCEESV,)
GOTO
GOTO 1 0 0 0
200
202
FORMAT(3XqF10.4r lX98A6~9X9F10.4~1X*RA6)
FORVAT(1Xv'M'rF10.4,1X98A6,7X9'M8*F10.4,1&,8A6)
48
20 1
203
3000
3001
4000
4001
5000
5001
6000
6001
7000
7001
1000
49
MODE D'EMPLOI DU PROGRAMME M E D T
CARTES CONTROLE
'
Les cartes de contrôle présentées ci-après commandent, pour chacune des v' variables, le calcul de la
médiane de I'échantillon des données e t celui des limites (inférieure e t supérieure) de confiance de
lamédiane,aurisque
5%. Le calcul n'est fait que pour un effectif reel supérieur à 10.
EXEMPLE D'ENTREE
50
UNITE
UNITE
UNITE
UNTTE
UNTTE
UNITE
UNTTE
UNITE
UNITE
UNITE
UNTTE
UNITE
IJNTTE
UNITE
UNITE
UNTTE
UNITE
UNTTE
UNITE
UtvITE
UNITE
UNITE
11
12
13
14
15
16
17
la
19
20
21
22
23
24
25
26
27
28
29
30
31
32
C/N
PH
TP
3
(3F4.0)
12.34.9e30
13.34.5a23
15.0
4.5
.21
13.1 4.9 a 5 2
16.33.9
-41
12.24.4-1.
12.24.6
.65
9.3 6.1 - 7 8
10.0 5.4 .41
11.1 6.2 - 5 4
9.2 7 . 4 -1.
9.4 7.5 -1.
9.9 5.P -1.
10.7 5.8 -1.
14.16.3
.40
9.3 7 . 6- 7 5
14.84.9-30
1 3 . 74 . 5- 2 3
e21
13.44.7
1 3 . 2 4.8 .§2
15.83.9
-41
12.24.7-1.
13.4
4.6
-65
10.55.9.78
10.0 5.3 - 4 1
.§4
10.6
6.0
10.27.2
-1.
10.1 7 . 2 -1.
12.3 5.9 -1.
11.1 5.6 -1.
1 2 . 46 . 4- 4 0
1 0 . 4 7.6 -75
51
.
.
EXEMPLE DE SORTIE
9.2000
9.3000
tJidITE 11
UNTTE 8
9.3000
IJKTTE 16
IJNJTE 1 2
UNITE 13
1 o . n o n o UNITE
9
1 0 . 0 0 0 0 IJNTTE 2 5
1 0 . 1 0 0 0 IJNTTE
1 0 .?O00 U N I T E , 2 7
1 0 . 4 0 0 0 UNTTE 3 2
1 0 . 5 0 0 0 UNITE 24
1 0 . 6 0 0 0 , UNITE 25
10.7000 UNTTE .14
1 1 . 1 0 0 0 IJI'IITE 1 0
1 1 . 1 0 0 0 UIJITTE '30
1 ? . P 0 0 0 UNTTE 6
1 2 . 2 0 0 0 UNITE 7
12.?onr) ~ J N J T E2 2
1 2 . 3 0 0 0 W T T E 29
1 2 . 4 0 0 0 UIt!ITE 31
12.ei)OO UMITE 1
13.1000 UhITE 4
1 3 . 7 0 n o IJNTTE 2 0
1 3 . 3 0 0 0 UNITE 2
1 3 . 4 0 0 0 UNITE: 19
1 3 . 4 0 0 0 UNITE 2 3
1 3 . 7 0 0 0 UNITE
IH
1 4 . 1 0 0 0 UNTTE 15
1 4 . 9 0 0 0 IJlJTTE 17
IL;.OO00 WJITE 3
1 s . p o n n UNXTE 21
l6..3000 UNITE 5
'
4.4000
9.4000
a
1
S
52
EXEMPLE DE SORTIE (suite)
3.5000 (JtVTTE
3.9000UNITE
4 . 4 0 0 0 UNI r E
4.c000 CJYTTE
4.5onn
4.5000
&.hono
1
Y
M
S
4.6000
4.7000
4,7000
h.8000
4.9000
4.0000
4.5000
5
21
6
2
IINITE 3
UNTTE 18
L J N T T E7
UNTTE 2 3
U"1TTE 1')
IJNTTE 22
UNITE ;In
UNITE 1
UNTTE L
UNTTE 17
5,3000 UNITE 25
SeLt000 U N I T E 9
5.c000 t J N I T E 7 0
S.PO00 UNTTE 1 3
5 r 9 0 0 0 UNITE 1*
C;,aono W I T T E 24
5 , Q O G O U f i J t t E 29
h.onoo
6.1000
U W T T E 26
IJtJTTE
$ * 2 0 0 0 br!ITE
6.7000 UNITE
h.LO00 IINTTF:
7.2000 UI\!XTE
7.2000 UNTTE
7.L.000 Ur4XTF.
7.5000 UNITE
7.6000UNITE
7.6000 UNITE
A
10
15
31
27
28
11
12
14
32
53
EXEMPLE DE SORTIE (suite e t fin)
TRAVAIL 1
F I C k l I E P 13-06-43
TP
..
-1.0000
-1.0000
-1,0000
-1.0000
-1.oorlo
-1.0000
-1.0000
-1.0000
-1.onoo
- 1 .O000
.?IO0
.?~oo
.2300
.73(JO
.3060
.3llflO
1
.4ono
.4000
.4100
.41no
M
.41rlo
M
-4100
.5200
S
.520O
.5400
,5400
.65CiO
.65Ii0
.7500
-7500
.7HCO
.7n00
IJMITE
IJNIITE
UI!TTE
UIJTTE
h
11
12
13
UNTTE. 1 4
I.INTTE 2 2
IJNTTE %7
iJt!TTE 23
U N T T E 29
IINTTF 30
UNTTE 3
IJNTTE 14
IJNTTE 2
UNITE 18
Ub!ITE 1
UNITE 17
UNJTE 1s
UkIITE 31
UNITE 5
IJNTTE 9
UNITE 21
IJNITE 25
1IYlTE 4
l!NTTE 2 0
IJNITE 1 0
UI\ITTE 76
IINITE 7
UNITE 23
1INTTE 1b
LJNTTE 3 2
UIdTTE II
IJI~!IT€,2 4
Une synthèse des résultats apparaît in fine ; elle reprend, en regard du libellé de la variable, l'effectif
réel, la limite inférieure de confiance de la médiane, la limite supérieure de confiance e t I'identification du fichier.
55
BlBLlOGRAPHiE
Le numéro placé devant chaque référence bibliographique indique :
(1) l'utilisation de Palgorithme de distance de D K A T, conjointement
avec un processus de constellation différent de H M V T, que nous avons
programmé mais non inclus dans ce texte;
(2) l'utilisation de D K A T;
(3) I'utilisation de H M V
(4) Sutilisation de D K A T, H M V T et C S A T;
(5) Sutilisation de H M V T et C S A T;
(6) les sources bibliographiques.
AUBRY, A.M., VAN DEN DRIESSCHE, R., BAUZON, D., PERRAUD, A., DOMMERGUES,
tests of rank correlation in the study of biological.and chernical data from tropical forest soils (Ivory Coast).
Bull. Ecol. Res. Comm., Stockholm, vol. XVII, 433-442.
Y.- 1973 - Measures of rank distances followed by repeated clustering and
BAUZON, D., VAN DEN DRIESSCHE, R., DOMMERGUES, Y. - 1968a - Caractérisation respirornétrique e t enzymatique des horizons de surface des sols forestiers. Sc. So/, vol.VIII, 2, 55-78.
BAUZON, D., VAN DEN DRIESSCHE, R., DOMMERGUES, Y . 1969 - L'effet litière. I - Influence in situ des litières forestières sur quelques caractéristiques biologiques des sols. Oecol,
Plant., vol. IV, 99-122.
~
BAUZON, D., AUBRY, A.M., VAN DEN DRIESSCHE, R., DOMMERGUES, Y. - 1975 - Contribution à la connaissance de la biologie des sols de la station PB1de Lamto, Câte d'ivoire. Rev.
Ecol. Biol. Sols, sous presse.
BEYER, W.H. - 1966 - Handbook of tables for probability and statistics. Chemical Rubber,
Cleveland, 517 p.
COINTEPAS, J.P. - 1974 - Application d'une analyse rnultivariable à quelques problèmes de
pédogenèseen Tunisie. Cah. ORSTOM, sér. Pédol., vol. XII, 2, 145-163.
GIRARD, M.C. 1968b Approche statistique de la notion desérie,sur un exemple pris dans
la bordure septentrionale de la Beauce. Thèse, Faculté dessciencesdeParis, 2 vol., 417 p.,
rnultigr.
~
~
GIRARD, M.C.-- 1969 - Statistiques e t pédologie détaillée. Introduction de la mesuredes distances Ag de Hiernaux. Sc. Sol, 1, 37-62.
GIRARD, M.C. - 1973 - Les images ballons e t les études du milieu naturel. Coursde/'ÉCole
d'Été du CNESà Tarbes, C.N.E.S.,Paris.
GIRARD, M.C. 1974 Traitement informatique des unités de paysage détectees par photographies-ballon. Edition spéciale du Service de stat/st/quede la Pologne, Varsovie.
~
~
GIRARD, C.M., GIRARD, M.C. 1973 - Interprétation des photographies-ballon. Comparaison
avec les petites régions agricoles. Soc. fr. Photogramm., 52, pp. 23-36, St-Mandé.
~
56
GIRARD GANNEAU, C.M.. GIRARD, M.C. 1974
- Traitement de l'information fournie par
des photographies prises en ballons stratospheriques. Applications en Agroliomie et Amknage-
ment du territoire. Coll. rntcrtrat. cclnformac/que et Enr~ironnen~enra, Arion
GIRARD GANNEAU, C.M., GIRARD, M.C. - 1974 - Photographs trom balloons : tlleir use in
agronomy and management of environment. 9t international Symposium on Remote sensing of
environment. April. A m . Arhor. Michigan.
GIRARD, b1.C.. GIRARD-GANNEAU, C.M. - 1974 Les photographies prises a haute altitude.
1,
Leur interprétation. Un exemple : les photographies-ballons. Rev, pl~oto-~nrerpidratton,
Technlp, Paris.
GIRARD, C.M.. GIRARD, M.C.
phère Paris, Masson, 186 p.
.
1975 -Application de fa délédétection ri l'étude de la bios-
HIERNAUX, J. . 1965 - Une nouvelle mesurede distance anthropologlque entre populations
utilisant simultan6ment les fréquences geniquesdes pourcentages de traits descriptifs e t des
moyennes metriques. C. R. Acad. Sc., Paris, 26 Opt., 1748-1750.
KENDALL, M.C., STUART, A. - 1966 - The aclvanced theory of stahtrcs. vol. 3, Desrgn and
analysrs, and tirne s e r m . Griffin, London, 552 p.
SPEARMAN. C. - 1904 - The proof and measurement of xsociation betweerl two things. Am.
J. PSvch., 15, 72-101.
VAN DEN DRIESSCHE, R. 1964 - Mesurede la dissociation analytique des sol!; par les distancesg6néralisées Dz de Malahanobis. ORSTOM, Bondy, 14 p. multigr.
VAN DEN DRIESSCHE, R. - 1965 - La recherche des constellationS de groupes à partir des
distances génerdis8es D' de Malahanobis. B/om -Prax., vol. VI, 1 , 30-4/.
VAN DEN DRIESSCHE, R. - 1966
&r. Pédol., vol. IV, 4, 91-96,
Un p r o b l h e de classification numgrique. Cah. ORSTOM,
VAN DENDRIESSCHE, R. - 1971 . A soi1 information system. In : First SernniarandGeographical Information SVstem, Vancouver, Simon Fraser University, decembet'1971, duplicated,
4 P.
VAN DEN DRIESSCHE, R. - 1972 Current vvork on computerized soi1 ddta processing in
ORSTOM. FA0;UNESCO Ad Hoc
contultation on computerized soi1 data interpretation for
development Purposes.Rome, 19-21 avr. 1972, 1 p.
~
VAN DEN DRIESSCHE, R., GARCIA GOMEZ,A.
profils. Rev. €col. Biol. Soi. vol. lx, 4, 617-628.
-
1972
-
Distances non-parametriques entre
VAN DEN DRIESSCHE, R. - 1973 - Un systhmeanlangage naturel pour banque dedonnées
pécfologiques. In : Semaine d'etude Sol et Fertilisation, sept. 1973, Faculte desSciences Agronomiques de I ' i t a t à Gembloux.
VAN DEN DRIESSCHE, R., GARCIA GOMEZ, A. - 1973Comparaison multivariablenonparametriqua des profils sans nomenclature des horizons. Cah. ORSTOIlil, sdr. Pédol., vol. XI,
3i4, 257-2G4.
VAN DEpJ DRIESSCHE, R. - 1974 - La Ijdnyue de données pédologiques de I'ORSTOR/I. Cah.
ORSTOM, sér. PPdoI., vol. XII, 1. 125-132.
O.R.S.T.O.M.
Direction générale :
- 75008 PARIS
24, rueBayard
Services Scientifiques Centraux :
Service Central de Documentation :
70-74, route d'Aulnay
-
93140 BONDY
Imprimerie S.S.C. BONDY
Editeur -
- O.R.S.T.O.M.
DBpbt l6gal : 3e trim. 1976
- I.S.B.N.
2-7099-0420-9
-