t CUMMODURE 8U32 PROGRAMME LA TT.TCES ...
t
CUMMODURE 8U32
PROGRAMME LA TT.TCES BALANCEES
3-4-5
Lattice 543
c
Dans le cadre de dispositifs experimentaux de type ¡°lattice balan-
C&i3 3 x 3 ; 4 x 4 ; 5 x 5¡± d o n t le p l a n d¡¯expdrience sur le terrain est con-
forme aux mod&l@ en annsxes.
Le programme calculs, aprhs a v o i r i n t r o d u i t l e s donnees d a n s l¡¯ordre
(r¨¦p6titions p u i s t r a i t e m e n t s ) d¡¯un t a b l e a u d e r6sultats classiques :
- Les moyennes 3631186 et a just8es de chaque traitement ;
- Les divers termes de 1 ¡®analyse de variante de f¡¯eseai.
LB programme effectu8 auSSi un ChC¨¤Semsnt arithm8tiqus
des diverses
moyennes et compare les diverses moyennes 2 ¨¤ 2 par le test de Keulhs.
Annexe 1
@RH D¡¯EXPERIENCE LATTICE 3 x 3
(No cle traitement correspondant)
i
t
! -.
t
-
7
!
I
i
I
II
I
I
III
t
-
?u
I
I
!
I
r----- ?
I
i
t No
No ordre
ordre rc4p
dP I Na ordre
ordre r8p f
, No ordre rBp
r&p t1
! RBp.
RBKI.
No o r d r e
ordre rdp
rdp ,I
I
r-n
-i
!
¡®!
2i3 1
I ?nomq--l-t------- ?u ?----
,BI,,I,2i3f,f4t7~1f5i91
¡¯
¡¯
?I 8f6
;
I Lr- ?---------rr.rrr--L
10!2f4f5f6,~2f518f7f2i61
¡¯
!
¡¯
1
U! v3 !
L c,.-I-----------m
I
1513
.
17 je 19 33 f6 19 1.4 18 f3 ; 7; 513
¡¯.!
PLAN D¡¯EXPERIENCE LATTICE 4 x 4
(No dn traitement corr sspondan t)
RBP, ; No ordre
i
NQ
ordre
;
Mo ordre
2
PLAN D¡®EXPERIENCE LATTICE 5 x 5
(No de traitemant correspondant) tempe calcul 2mn 59
1
1
---I
¡®-
f
I
i
II
I
III
!
I
f
-!
f
-!
I R~P. l
No ordre r4p, I
No ordre r6p.
!
No ordre r8p.
1
!
-
-
!
I
I
-1
! RGp* !
No ordre r6p. !
No ordre r4p.
r
No
ordre r6p, i
3
Annexe IV : Tableau d'arrangement desblocs introduits en data (Lattice Y x 3)
QeP I
I
f
-T-=-I
h.Fa?J r
I 11
1 III [ IV 1
-
-
-
-
-
1
1
1
!
!--Y--!
?
J
¡¯ i 4 !
7 1 ¡®0 !
I -
e
-
m
-
!
!
f
! -f
! 2 J 1 J 5 J 6
J
¡®1
J
! -!-----
!
J
I
I
!
3 J
¡¯
!
6
I
9
!
¡®2
!
! -,- J -!-- I -
J
4 J 2 J 4 J 9,11;
-
w
-
m
-
-
J
I
!
I
I
I
J
5 J 2 J 5 J 7 1'2 I
J -
-
J
J -
-
C
I
-
!
I
I
J 6 ! 2 ! 6 ! 8
1 l0 1
! ---L3
I
! -
,
-
- l -1
! 7 J I5 J 4 J 8 !I*f
!------
l
I
1
J
1
J
o J 3 ! 5 ! 9 1 'O I
! -w
!
I -
-
-
-
J
f
J
f
9 I 3 I l5 I 7 1 ¡®I-J J
Annexe V : Tableau d'arrangement des blocs introduits en data (Lattice 4 x 4)
L
I
J-!--yry-!--?-!-!
J
¡¯
J
1
J
9
f ¡®3
! l7
!
-
- --a
J
J
-
-
!
!
J
J
!
1
2 !
¡¯
!
6
1 10 I 14 f 18 1
I -
-
!
J -!---
!
l
!
! 3 J ' 17 11'
! -15
J l9 !
-
- -
-
!
!
!
J
I -
-
!
J
l
4
t
¡¯
!
8 J
¡®2 J
¡®6 J
20
J
-
-
- -
-
-
!
I
I
!
!
f
J
!
5 ] 2 1 5 ] 10 ] 15 120 J
-
-
-
!
J
!
-
-
-
I
J
J
!
J
6 l 2 J 6 J 9
If
-If5 J
¡®9
J
-
- -
- -
-
!
I
!
1
Il
J
!
f7f2J7
!
l2
1
¡®3
!
¡®8
]
! -
,
- ! -
,
- !! -f----- J
1
8 J
2 J
8
]¡±
!
l4
J
.17
J
! -!-----
!
!
I
!
!
!
9 J 3 1 5 , '1 ! 16 , 18 ,
-
-
-
J
-
-
-
!
f
!'o , 3 I 6 ;12 : 15 f 17 ;
-
- -
- -
-
-
-
-
J
-"IL13-
1
,12 ! 3 f 8
f 10 ; 13 119 *1
!------
!
!
] 13 ] 4 J 5 il2 : 14 fl9
f
-
- -
v
s
-
.
-
f
J
1 14 I 4
f
6 f 11 ; 13 ;20
f
I-----'--
f
!
f 15 , 4 , 7 ; 10 f 16 f 17
;
--v
!
-.--
!
J
116 ; 4 , 8 , 9 115 f 78 ;
.
-
¡°-111
.-
..-¡°-..w.-I..me. ?.--..
-1-1-1-l-1--------..w ++l++ --..w-----*-.--I
.-
-]-]+j-p;
-
Im-------wm-u-w N N I
m--w -3 m 1---.w + I.---
4
¡°v-m
m
--
h,
..-
0)
-
.-
;o
m-
hl
e....
a
.-
h) -2.
.-
s:
.-
--
a
a
N N
. .
m
-
N
I-L- --a--- w w w l 1w---w- a
4
1
I-m--m.. a J 4Julc I 1--m--e -4 I I
0:
_l
I ruw
--m--c
w
_I
l-
---
m
--m-u.-
h) *
-
..L-
---------c
M-L
-.-
++++ --.---.-w-M..
-i-i-l-l- w.-m8.. m~41*p+ --.--.------ VI+I+ -----*.D----
---M---m-.-
f: e-.-m--
a
la
I hJ
I
I
---.- P e I-w-w a u3 1m--w J .a h)2 1o.-eIam J ht a w
m--I
---*a
r:
m
m-w.-.
c
a---
# m
5
Listinq .du proqramme "Lattice 543"
l@
rem
"LATTICES BALANCEES"
20
data 1, 4, 7, ID.......* )-e.---w- Tableau des blocs lattice 3 x 3
21
data 3, 6, 7, 11
>
30
data 1, 5, 9, 13*...***.
)
Tableau des blocs lattice I x 4
33
data 4, 8, 9, 15, le.... )-------
40
data 1, 6, 11, I~I.......
) ------Tableau dss blocs lattice 5 x 5
46
data.... 11, 18, 22, 29
>
100
dim. x (6,25), u3 (4,9), ~4 (5,16), u5 (6,25)
120
input "voulez vous conserver les donnees WI OU NON" ; a4
130
iF a$ = "non" then 200
140
print ¡°introduisez les donn¨¦es en DATA ¨¹ partir de 5OOD"
150
print "sn respectant l'ordre du tableau de resultats soit :"
160
p r i n t "d'abord toutes les repetitions du traitement i puis celles"
170
p r i n t "du traitement 2 etc.. ..¡±
175
print "POUR CONTINUER L'ANALYSE COMPOSER pun 260
100
stop
200
input "quel est le nombre de traitements a tester¡± ; k2
2 1 0
ko =o~g (K2) : r = ko + 1
220
for j = 1 to k2
230
for i = 1 to r
240
print "Resultat de la rr3pdtition No "i" pour le traitement No" j ; :inputx(i,j
2 5e
next i,j
2 ci!5
goto 360
260
input "quel est le nombre de traitements h tester" ; k2
265
ko = sqr (k2) : r z ko + 1
27 0
dim x (6,25), u3 (4,9), u4 (5,16), ~5 (6,25)
280
for j = 1 ta 9 : for i = 1 to 4 : read 03 (i, j) : next i,j
290
for j = 1 to 16 : fo:r i = 1 ta 5 : read u4 (i,j) : next i,j
300
for j = 1 to 25:for i, =lto6:read u5 (i,j> : next i,j
3 1 Cl
for j = 1 ta k2
32Cl
for i = 1 to r
3 3 ci
r e a d x (i,j)
34U
next i , j
34 5
goto 305
360
for j = 1 to 9 : for i = 1 to 4 : read ~3 (i,j> : next i,j
370
for j = 1 to 16 : for i = 1 to S : read u4 (i,j) : next i,j
380
for j = 1 to 25 : for i = 1 to 6 : read u5 (i, j) : next i, j
383
Print ch& (19) chr$ (19) chrp( (197).
6
385
rein ¡°choix du plan d'experience"
3 9 0
dim ul (r, k2)
4 00
on ko-2 goto 410, 420, 430
4 10
for j = 1 ta k2 : for i = 1 to r : ul (i,j> = u3 (i,j) : next i,j : goto 440
4 2 0
for j = 1 do k2 : for i = 1 ta r : ul (i,j> = u4 (i,j> : n@Xt i,j : goto 4dO
43c!
for j = 1 to k2 : for i = 1 to r : ul (i,j) = u5 (i,j) : next i,j
44 c!
rem "verification des donn¨¦es"
450
if k2 = 25 ttwn 550
461?
for j = 1 ta k2
46s
p r i n t j ;
47(3
for i = 1 to r
4 3 0
print tab (i*10) x (i,j) :
4 85
next i
49 u
Print
4 9 5
next j
s 00
print "v¨¦rifiez le tableau figurant ¨¤. l' ?cran et notez les erroure 6ventuelJes"
510
print "en rep¨¦rant les numeros de r¨¦p8tition et de traitement correspondants¡±
52C
input "POUR CONTINUER L'ANALYSE COMPOSER C" ; ad
5 3 0
if a6 = "c'" then 710
5 5 0
for j = 1 to 20
555
p r i n t j ;
S60
for i = 1 to r
5 7 0
print tab (i-l 1 10+2!) x (i,j) ;
58Q
next i
5 8 5
print
5873
next j
590
print ¡°v¨¦rifiez le tableau figurant ¨¤ l¡¯¨¦cran et notez les erreurs eventusl&s"
600
print ¡°en repdrant les num0ros de riSp¨¦tition et de traitement correspondants¡±
510
input "POUR CONTINUER LA VERIFICATION COMPOSER c" ; a+
6 3 0
if ad = r'~" then 640
fi40
for j = 21 to 25
6 4 5
print j ;
6 5 0
for i = 1 to r
660
print tab (i-l) * 10+2) x (i,j) ;
6 6 2
next i
6 6 5
print
C!I70
next j
6 7 5
print
680
print "v¨¦rifiez le tablfoau figurant ¨¤ l'¨¦cran et notez les errwrs ¨¦ventuelles"
690
print "sn rep¨¦rant les numiaro ? de repetition et de traitement correspondw?,i::
700
input "POUR CONTINUER L'AN4LYSE COMPOSER c" ; a$
'If! 5
if ag! = "c'. then 710
71C
input "combien de currections sont n¨¦cessaires" SS
715
if SS = 0 then 750
72C!
forua = 1 to ss
730
input "correction
: No de repetition, !\\jo dc traitement, N11.valeur""i,j,r¨ªi,j)
742
next ua
750
rem "calcul des divers totaux"
800
rein " l-total par r8p8titionN
PIO
for i 2 ? ta r
820
for j = 1 to k2
830
tr (i) = tr (i) + x (i,j>
846
next j,i
850
rem "2-total par traitement non ajuste"
i') t;s
dim ,t (k2)
860
for j = 1 to k2
i37U
for i = 1 to r
980
t Ij> = t (j> + x (i,j>
a90
next i,j
900
rem "3'total q¨¦n¨¦rals'
?10
for i = 1 to r
320
for j = 1 to k2
9 3 Il
g = 9 + x (i,j)
Y40
noxt j,i
9 50
rcIm " - 4 total par bloc incomplet"
C] ytj
dim tb (r*ko)
310
for u = 1 to r *ko
9 ? :.l
for i 2 1 to r
380
for j = 1 to k2
990
if ~1 (i,j) = u then tb (u) P tb fu) + x (i,j)
1000
next j, i, u
1010
retn "4- total des blocs incomplets pour chaque traitement¡±
1015
dim bt (k2)
1020
for j = 1 to k2
1030
for i = 1 to r
040
for u = 1 to r * ko
0513
if ul (i,j) i L' then bt (j) = bt (j) + tb (u)
060
next u, i, j
36G
rem "facteur d'ajustamsnt¡°
0
1365
dim IU (k2)
1370
for i = 1 to k2
1380
u)(i) z ko *t(i) - (ko + 1) * ht(i) t g
1390
next i
1400
rcm "variations diverses -1 variation totale"
1410
tc = g?2,&2/&420 far i = 1 t0 r
1430
for j = 1 to k2
1455
vt = vt t x (i,j>t 2
1460
n e x t j , i
14 7 0
vt = vt - tc
1480
rem "2'variatione dUBs aux rGp8titions"
1490
for i = 1 to r
1 SO0
vr = vr c t r
(i)fZ
1510
next i
1520
vr 5: wr/ko/ko - tc
1530
rem ¡°3-variations dues aux traitements"
1532
for i = 1 to kZ
153s
vc - v c t t(i)¡®2
1537
next i
1540
UC = vc/ (ko+l) - tc
1545
FB:I~ ¡°4-variations
dues aux blocs ajustfis"
1550
far i = 1 to k2
1560
vb I vb + w(i) 7 2
157[3
next i
15813
vb = vb/ko*3dko+l) = vb = vb/ko?3/(ko + 1)
1590
rem ¡°erreur intrablocs at facteur d'ajustement (uu
1600
ib = vt - vr - vc - vb
1610
mb = vb/fko+%l)
lG20
im = ib/(ko-1) * (k2-1))
1630
uu = (mb-im)/k2/mb : if uu(= 0 thon uu = 0
1635
dim mc (k2f, i (k2) : tt = 0
1540
rem ¡± totaux a jus t¨¦es par traitement¡±
1645
dim ta (k2), ma (k2)
1650
for i = 1 to k2
3 aio
ta(i) = t (i) + uu * tu(i) : ma (i) t ta (i>/r
1470
next i
1680
f u r i = 1 to k2
1690
tt = tt t ta (i)? 2
17flO
fwxt i
rem ';ci?lctil des wariancss et gain de ps¨¦cisionVf
wt = (tth - tcy(k2-1)
et = ib/(k2-l)/(r-2) : ex = et * (1 + kr-* uu)
~p = 100 ++ ((vb + ib)/(k~&Z-?f/ox) : $p = int (gp*lOO+ ,S)/lOO
rem "impression des rosultats"
print ch& (19) ch& c19j chr+ (747)
print ri
LATTICE "ko" X "ko : print
print tab (5) "moycnfws obs~rv¨¦9s'! tab (40) "moyennes corrig6Ses" : print
dirs mr (k2)
for i = 1 to k2 : mr (i) = t (i)/'c : ncxt i
for i = 1 to k2
print tab (3) "A" i" ' mr (i), teb (44) ma (i)
if i/20 = int (i/20) then 1840
naxt i
Print "notez les rdsultats figurant Gr l'¨¦cran f POUR CONTINUER"
input "LA NDTATIOM COMPOSER c SINON COMPOSER s" ; ad
if ai = "c" then 18313
rem "analyse de variancer'
print : print tab (2s) "ANALYSE DE VARIANCE" print
print "facteur" tab (25) ¡°variante" tab (5~0)~ D . D ?? L . "tab(65)" F ,
caIcu14" : print
19ocj
print iirBpEtitions" i;e"J (25) vr/ko tah (50)ko tab(65) vr/ko/et
? ¡®1¡® 1 fi
print "blocs adj." tm (25) vb/;kZ-1) tab(5C) k2-1 tab(65) vb/(k2-l)/at
1920
print ¡°traitements¡±
tab(2S)urt tab(50) k2-1 tab(65) wt/ex
193c
print ¡°erreur¡± tab(25) ex tab(50) (k2-1) * (r-2)
1 c.1 it, I-J
pri nt : print ¡°gain de pr¨¦cision s "gp - 100" $ 1, : mm = g/r/k2
¡®1 3 4 2
ax = sqr (ex) t p r i n t ¡°coc~ff. do v a r i a t i o n zl¡¯ int((ax/mm) * 104+~~5)/1~0¡± $ ¡±
1945
print "moyenne gBnBr2rle ='l mm : print
195@
p ?int tab(20) reclassement des moyennes ajust¨¦es"
1960
for i = '1 ta k2 : i(i) = i : me(i) = ma (i) : noxt i
1770
for i = 1 to k2-1
1980
for j = 1 to k2-'l
1930
if m3 (j) = no (j+l> thcn 2 ?40
2000
c = ma(j) : d = i(j)
2026
ma(j) = ma(j+l) : i(j) = i(j+l)
2 0 3 0
na(j+l) z c : i(j+l) = d
2 0 4 0
next j
2350
next i
2060
for i = 1 ta k2
10
Print !fAil i(i)" (il ;
next i
dim c3 (a), d3 (8), ct (15), d4 (15), cl (k2- 1), dl (k2-i), ~5 (Lti), d5(24;
on ko-2 9uto 2030, ZW5, 21133
3 c3 fi> = 3.CL........)
3
)tzr,bla de Keule pour 16
. ..r....d3 (8) = 6.22 )
d.d.1.
c4 (1) = bb.,b..*~**L.r.~...,*, ftable do K~ule pour 40 d.d.14
. . . . . . . . . . . . . . d4 (15) = 6,t)2 :)9oto 2115
55 (1) =.**.*o.* 1)tablc ds Kwls pour 120 d.d.1.
\\
..d5(24)= 6.02 '
on ko-2 yoto 212C, 2160, 2592
for i = 1 to k2-?
cl(i) = d(i) -+ sqr (ex/*)
dl(i) = d3(i) + sqr (ox/r)
next i
goto 2193
for i = 1 to k2-?
cl(i) = c4(i) + sqr (ex/r)
dl(i) = d4(i) * sqr (ex/r)
naxt i
yoto 2193
for i = 1 ta k2-1 : cl(i) = c5(i)*eqdox,/r) : dl(i) = d5(i) * sqr(ax/r);noxt i
print : input "waU1ez:-vous le test de Kouls ¨¤ la proba. C.05 OUI ou NGN a$
if a4 = "non" then end
n ZZ l:i=Otj=O
if j + n > k2 thon 2262
print i (j + 1) i
j z j t 1
if j>l thon print i (j) ;
i =i.tl
ifj+n=
k2 + 1 thon 2262
if ma (j+n) . ma(j)<= cl(i) then print II=" i(j+n) i : n = n + 1 : 9oto 2235
if ma (j+n) - me (j),) cl(i) then i = Cl : n = i : print : pt0 2230
print : input '*voulez voue le test de KEULS &I 10 proba. II.01 OUI ou NON";o$
if a$ = "non" then end
print t&b,(lO) "test da Keule PrOba o*[jl"
n = l:i=O:j=O
if j t n)k2 then end
229¡¯3
print i(j + 1) ;
2 2 0 0
j = j + 1
2301
if j >? then print i (j) ;
:,z -5 0 2
i =i+l
2305
ifj+n=
k2 + 1 then end
2 5 1 0
if mz (j t nj - ma(j){= d 1 (i) then print "=" i(j+n) ; : n z n + 1 : gotB.
2302
2320
if ma (j+n) - ma(j)>dl (i) ti&n i = 0 : n = 1 : print : goto 2330
12
- Formulation mathdmatique utilis¨¦e
. totaux :
-
-
r6sultat individuel : x )n ombra de traitements = k2
total lg¨¦n¨¦ral : g
ho
indice de la kattica =Vnbtrait>
total par bloc incomplet : tb
Il existe (r
x ka > valeurs tb
tata1 par rf5p¨¦tition
: tr
Il existe r valeur tr
total par trai teiilent
: t
If existe ko x ka
v a l e u r s t
. pour chaque traitement : total des blocs incomplets 03 il apparait bt
. facteur de correction pour chaque traitement : UI = kot - (ko + 7) bt + g
. kscteurs de variation :
. terme correctif ¨¦c = g
Y (r x kcr2)
. casr4 variation globale vt = c x2 - tc
. carr¨¦ variation r¨¦p¨¦tition vr = jtr2/ko2 - tc
. carr¨¦ variation traitement va = &t2/(ko + 1) - tc
. blocs ajustes carr6is vb = < w2/ko3(ko + 1)
. intrabloc somme des c'arrgs = ib = vt - vr - va - vb
. Variantes :
. blocs ajust¨¦s : vb/(ko2 - 1) Ir) a
. intrablocs vari¨¦s : ib/( 1(0*-l > -43 t
. facteur d'ajustement uu = a - et
liY2-y
si uu<O on consid¨¨re que uu = 0
correctif traitements
?
total ajuste traitement = t +wIo
. variante sur les traitements ajust¨¦s
variation traitements ajust¨¦s = ;5 t ajust4s totaux¡¯/r - tc
variante correspondante
: vz variation trait. ajust&/k2-l
. variante effective de l¡¯erreur : ex = et (1 + ka
:x uu)
Tableau d'analyse de variante
!
!
!
1
variation
!
Variante
!
D.D.L.
!
f
!
!
!
!
!
f
I
!
!
1 R ¨¦pdti tions
I
vr/ke
1
ka
! ur/ko/et
t
1
!
!
!
1
?
l
!
!
I
I Olocs ajust¨¦s
!
vb/(ko2 - 1)
I ko2 - 1
!vb/(ko2-l)/et !
!
!
!
!
!
!
!
!
!-
I
v a r i a t i o n t r a i t
! Traitements ajus@IajustBs/'ko2 - 1
I ko2 . 1
1
variance/ex
!
!
I
!
!
1
I
!
!
t
! Erreur ajust¨¦e
!
63 x
!(k 24)(k 4cll) !
1
!
I
I
I
!
13
- jode d'emploi du proqramme "Lattice 543"
load "lattice 543" return
loading >
rsady
) 3 mn 05 sec
run return
Voulez vous conserver les donn¨¦es OUI ou NON ? oui return
introduisez les donn¨¦es en DATA ¨¤ partir de 5000 en respectant l¡¯ordre ;
du tableau do r$sultats soit :
!
d¡¯abord tout¨¦s las rbp¨¦titions du
2 etc..,
traitement '1 puis celles du Eraitem ?nti
POUR CONTINUER L"ANALYSE COMPOSER ru-l 260
1
I
break in 180
1
f
ready
I
Tableau de r¨¦sultats
I
1
LATTICE
J
!
!
!
I
II
l
I
3 x 3
1
r
IX1
!
IV
f
!
-!
.
f
!
I
!
! Traitements , 1 I
?-!
2299
*
23,5
;
16,9
]
24,2
f
!
!
f
1'
2
34,4
;
3284
I
!
i
!
!
I - ?
?
-3
1
2697
!
3696
I
!
(1
I
!
!
! 3 !
2896
1
31,4
1
26,2
f
27,9
;
f
!
!
II
-!
2891 -!
28,6
!
I
!
!
!
! 4 !
[
I
25,1
1
3397
,
f
!
-!
i
I
!
II
I- ?
!
! 5!
3397
i
30,7
;
2693
!
3091
1
I
!-!
29,0
i
;
34,3
I
!
I
!
'9
I 6
32,0
I
1
.
31¡®0
f
I
!-
If
? - I
?
!
1
!!-
! 7 f
2591
,
30,2
1
24¡®9
f
29,4
;
I
!-!
If
-1
I
1
!
!
! 8,
31¡®6
f
33,l
i
1
-*
27,8
1
3119
I
!
!
-1
II
r
!
Pf
2996
1
33,5
;
29,4
f
29,7
f
\\
Selon Cochran et Cox.
1
"Experimental designs" Page 304
1
I
SO¨¹O data 22.9, 23.5, 16.9, 24.2, 34.4,......., 27.9 rtiturn
I
(3)
t
5001 data 28.1, 20.6, 25.1, 33.7,....,.., 25.1, 30.2 return
I
3hi.s
t
*un 260 return
!
(4)
Dutil est le nombre de traitement ¨¤ tester ? 9 return
I
(5)
1
22.9
23.5
16.9
24.2
I
2
34.4
32.4
26.7
36.6
f
3
28.6
31.4
26.2
27.9
I
4
28.1
28.6
25.1
33.7
I
5
33.7
30.7
26.3
30.1
t
6
29
32
34.3
31
1
7
25.1
30.2
24.9
29.4
I
a
3 1 . 6
3 3 .l
2 7 . 8
31.9
1
9
2 9 . 6
3 3 . 5
29.4
29.7
I
14
!
No R6f.
f
v¨¦rifiez le tableau figurant ¨¤ lIecran et notez 10s erreurs kventuclles
!
un rcperant les num¨¦ros de r&p¨¦tition et de traitements C:orrespondants
!
!
POIJR CONTINUER L'ANALYSE COMPOSER c ? c return
!
Cwnbicn de correction sont n¨¦cessaires ? 0 roturn
!
f
LATTICE 3 x 3
1
I
!
moyonncs
observ¨¦es
moyennes corrig¨¦es
I
A 1
21 .a75
21.875
I
A 2
!
32.525
32.525
!
Cl 3
23.525
28.525
!
P, 4
!
20.875
28.875
!
A 5
30.2
30.2
!
!
A 6
31.575
31.575
!
A 7
27.4
27.4
I
I
A 8
31.1
31.1
!
A 9
30.55
30.55
!
I
notez les r8sultats figurant ¨¤ 1"Bcran. POUR CONTINUER LA NOTATION
t
!
COMPOSER c SIPJON COMPOSER s 5' s rsturn
I
(9)
!
ANALYSE DE VARXANCE
!
I
r
facteur
variante
D.D.L.
F c a l c u l ¨¦
!
1
r-ipt! ti ti ons
34.4669266
3
5.37297662
f
Blocs adj.
3.40611107
8
,530971482
!
traitements
!
40.2911053
8
6.28089557
!
erreurs
6.4 148663
i6
I!
goin d e p r ¨¦ c i s i o n = - 1 5 . 6 3 $
I
cocff. d e v a r i a t i o n =
8 . 6 8 $
I
!
moy g¨¦n¨¦rale
= 29.1805556
!
Classement des moyennes a.iustdes
!!
Al< A? (A3 (A4<AS (A9( A04 MC: A2(
I
!
Voulez vous le test de Keuls ¨¤ la proba 0.05 OUI ou NON ? oui return
t
(9)
Voulez vous le test de Keul ¨¤ la proba 0.01 OLII ou NON ? non return
I
(10)
ready
i
Variantes dans le mode dtemploi du proqramme
1
I
ru17 return
f
(1)
voulez vous conserver les donnbes OUI ou NON ¡®1 non return
1
(2)
quo1 est le nombre de traitements ¨¤ tester ? 9 return
I
I
(31
K¨¦sultat de la rbp¨¦tition no 1 par le traitement no 1 ? 2 2 . 9 r e t u r n
I
(4)
,I
l
2
1 ? 23.5 return
!
;t
-3
n
-> ,..
15
1 NO H¨¦f
Resultat de la rep¨¦tition na 3 par le traitement no 9 ? 29.4 return
!
.
RQsultat de la r¨¦petition n" 4 par le traitement no 9 ? 29.7 raturn
1
I
Verifiez le tableau figurant ¨¤ 1¡¯6cran et notez les erreurs Qventuelles I
!
i.;n rep¨¦rant les numdros de rep¨¦tition et de traitements correspondants I
POUR CONTINUER L'ANALYSE COMPOSER c ? c
(5)
!Zambien de corrections sont n6cessaires ? 1 return
i
!
(6)
c 3rr ec ti on
: No d e r¨¦pAtition,
oJo d e t r a i t e m e n t , Nil. valeur ? 2,2,32,4 I
(7)
return !I
I
Les rQsultats apparaissent ensuite
B l¡¯ ?cran en tenant compte de la correction effectu¨¦e f
Remarques :
1) ¨¤ 1'5tape (6) de la variation, on s'¨¦tait aper ?u en v¨¦rifiant le ta- I
bleau de r¨¦sultats que l'on avait fait une erreur sur la 2Bme r¨¦pdti-.f
tion du traitement 2 c'est Pourquoi on a r¨¦pondu 1 B la question po&td
et ensuite on a introduit la correction ( r e p ) Nll.valeur(la virgule !
2 , 2 , 3 2 . 4
I
separe 2 donn¨¦es)
( t r a i t . )
!
!
!
2) Dans le cas d'une lattice 5 x 5, l'ecran ne pouvait pas contenir tout
le tableau de r¨¦sultats pour le v¨¦rifier, i l a p p a r a i t l e s 2 0 p r e m i e r s
traitements,
puis dans un deuxi¨¦me temps les 5 autres. Il en va de
I
m&me pour les moyennes, o¨´ il est alors n8cessaire de composer c pour1
voir apparaitre les 5 derni¨¨res moyennes.
3) Cans l¡¯exemple traite le coefft de correction
est negatif ce qui
I!
Explique ltQgalite des moyennes observdes et corrigees.
1
4) Si l'essai est non significatif on rgpondra non a la question No Y et:
l'appareil affichera ready.
5) En (7) l¡¯ordinateur calcule pendant 1 B 4 mn environ selon Ie type de! !
lattice (3 x 3 - 4 x 4 - 5 x 5).
I
(g) Le test de Keuls fait apparaitre certains classements inutiles,
I
il ne faudra pas en tenir compte dans l¡¯interpr¨¦tation.
Test de Keuls proba 0.05
1
7=3=4=5=9=8=6=2
CUMMODURE 8U32
PROGRAMME LA TT.TCES BALANCEES
3-4-5
Lattice 543
c
Dans le cadre de dispositifs experimentaux de type ¡°lattice balan-
C&i3 3 x 3 ; 4 x 4 ; 5 x 5¡± d o n t le p l a n d¡¯expdrience sur le terrain est con-
forme aux mod&l@ en annsxes.
Le programme calculs, aprhs a v o i r i n t r o d u i t l e s donnees d a n s l¡¯ordre
(r¨¦p6titions p u i s t r a i t e m e n t s ) d¡¯un t a b l e a u d e r6sultats classiques :
- Les moyennes 3631186 et a just8es de chaque traitement ;
- Les divers termes de 1 ¡®analyse de variante de f¡¯eseai.
LB programme effectu8 auSSi un ChC¨¤Semsnt arithm8tiqus
des diverses
moyennes et compare les diverses moyennes 2 ¨¤ 2 par le test de Keulhs.
Annexe 1
@RH D¡¯EXPERIENCE LATTICE 3 x 3
(No cle traitement correspondant)
i
t
! -.
t
-
7
!
I
i
I
II
I
I
III
t
-
?u
I
I
!
I
r----- ?
I
i
t No
No ordre
ordre rc4p
dP I Na ordre
ordre r8p f
, No ordre rBp
r&p t1
! RBp.
RBKI.
No o r d r e
ordre rdp
rdp ,I
I
r-n
-i
!
¡®!
2i3 1
I ?nomq--l-t------- ?u ?----
,BI,,I,2i3f,f4t7~1f5i91
¡¯
¡¯
?I 8f6
;
I Lr- ?---------rr.rrr--L
10!2f4f5f6,~2f518f7f2i61
¡¯
!
¡¯
1
U! v3 !
L c,.-I-----------m
I
1513
.
17 je 19 33 f6 19 1.4 18 f3 ; 7; 513
¡¯.!
PLAN D¡¯EXPERIENCE LATTICE 4 x 4
(No dn traitement corr sspondan t)
RBP, ; No ordre
i
NQ
ordre
;
Mo ordre
2
PLAN D¡®EXPERIENCE LATTICE 5 x 5
(No de traitemant correspondant) tempe calcul 2mn 59
1
1
---I
¡®-
f
I
i
II
I
III
!
I
f
-!
f
-!
I R~P. l
No ordre r4p, I
No ordre r6p.
!
No ordre r8p.
1
!
-
-
!
I
I
-1
! RGp* !
No ordre r6p. !
No ordre r4p.
r
No
ordre r6p, i
3
Annexe IV : Tableau d'arrangement desblocs introduits en data (Lattice Y x 3)
QeP I
I
f
-T-=-I
h.Fa?J r
I 11
1 III [ IV 1
-
-
-
-
-
1
1
1
!
!--Y--!
?
J
¡¯ i 4 !
7 1 ¡®0 !
I -
e
-
m
-
!
!
f
! -f
! 2 J 1 J 5 J 6
J
¡®1
J
! -!-----
!
J
I
I
!
3 J
¡¯
!
6
I
9
!
¡®2
!
! -,- J -!-- I -
J
4 J 2 J 4 J 9,11;
-
w
-
m
-
-
J
I
!
I
I
I
J
5 J 2 J 5 J 7 1'2 I
J -
-
J
J -
-
C
I
-
!
I
I
J 6 ! 2 ! 6 ! 8
1 l0 1
! ---L3
I
! -
,
-
- l -1
! 7 J I5 J 4 J 8 !I*f
!------
l
I
1
J
1
J
o J 3 ! 5 ! 9 1 'O I
! -w
!
I -
-
-
-
J
f
J
f
9 I 3 I l5 I 7 1 ¡®I-J J
Annexe V : Tableau d'arrangement des blocs introduits en data (Lattice 4 x 4)
L
I
J-!--yry-!--?-!-!
J
¡¯
J
1
J
9
f ¡®3
! l7
!
-
- --a
J
J
-
-
!
!
J
J
!
1
2 !
¡¯
!
6
1 10 I 14 f 18 1
I -
-
!
J -!---
!
l
!
! 3 J ' 17 11'
! -15
J l9 !
-
- -
-
!
!
!
J
I -
-
!
J
l
4
t
¡¯
!
8 J
¡®2 J
¡®6 J
20
J
-
-
- -
-
-
!
I
I
!
!
f
J
!
5 ] 2 1 5 ] 10 ] 15 120 J
-
-
-
!
J
!
-
-
-
I
J
J
!
J
6 l 2 J 6 J 9
If
-If5 J
¡®9
J
-
- -
- -
-
!
I
!
1
Il
J
!
f7f2J7
!
l2
1
¡®3
!
¡®8
]
! -
,
- ! -
,
- !! -f----- J
1
8 J
2 J
8
]¡±
!
l4
J
.17
J
! -!-----
!
!
I
!
!
!
9 J 3 1 5 , '1 ! 16 , 18 ,
-
-
-
J
-
-
-
!
f
!'o , 3 I 6 ;12 : 15 f 17 ;
-
- -
- -
-
-
-
-
J
-"IL13-
1
,12 ! 3 f 8
f 10 ; 13 119 *1
!------
!
!
] 13 ] 4 J 5 il2 : 14 fl9
f
-
- -
v
s
-
.
-
f
J
1 14 I 4
f
6 f 11 ; 13 ;20
f
I-----'--
f
!
f 15 , 4 , 7 ; 10 f 16 f 17
;
--v
!
-.--
!
J
116 ; 4 , 8 , 9 115 f 78 ;
.
-
¡°-111
.-
..-¡°-..w.-I..me. ?.--..
-1-1-1-l-1--------..w ++l++ --..w-----*-.--I
.-
-]-]+j-p;
-
Im-------wm-u-w N N I
m--w -3 m 1---.w + I.---
4
¡°v-m
m
--
h,
..-
0)
-
.-
;o
m-
hl
e....
a
.-
h) -2.
.-
s:
.-
--
a
a
N N
. .
m
-
N
I-L- --a--- w w w l 1w---w- a
4
1
I-m--m.. a J 4Julc I 1--m--e -4 I I
0:
_l
I ruw
--m--c
w
_I
l-
---
m
--m-u.-
h) *
-
..L-
---------c
M-L
-.-
++++ --.---.-w-M..
-i-i-l-l- w.-m8.. m~41*p+ --.--.------ VI+I+ -----*.D----
---M---m-.-
f: e-.-m--
a
la
I hJ
I
I
---.- P e I-w-w a u3 1m--w J .a h)2 1o.-eIam J ht a w
m--I
---*a
r:
m
m-w.-.
c
a---
# m
5
Listinq .du proqramme "Lattice 543"
l@
rem
"LATTICES BALANCEES"
20
data 1, 4, 7, ID.......* )-e.---w- Tableau des blocs lattice 3 x 3
21
data 3, 6, 7, 11
>
30
data 1, 5, 9, 13*...***.
)
Tableau des blocs lattice I x 4
33
data 4, 8, 9, 15, le.... )-------
40
data 1, 6, 11, I~I.......
) ------Tableau dss blocs lattice 5 x 5
46
data.... 11, 18, 22, 29
>
100
dim. x (6,25), u3 (4,9), ~4 (5,16), u5 (6,25)
120
input "voulez vous conserver les donnees WI OU NON" ; a4
130
iF a$ = "non" then 200
140
print ¡°introduisez les donn¨¦es en DATA ¨¹ partir de 5OOD"
150
print "sn respectant l'ordre du tableau de resultats soit :"
160
p r i n t "d'abord toutes les repetitions du traitement i puis celles"
170
p r i n t "du traitement 2 etc.. ..¡±
175
print "POUR CONTINUER L'ANALYSE COMPOSER pun 260
100
stop
200
input "quel est le nombre de traitements a tester¡± ; k2
2 1 0
ko =o~g (K2) : r = ko + 1
220
for j = 1 to k2
230
for i = 1 to r
240
print "Resultat de la rr3pdtition No "i" pour le traitement No" j ; :inputx(i,j
2 5e
next i,j
2 ci!5
goto 360
260
input "quel est le nombre de traitements h tester" ; k2
265
ko = sqr (k2) : r z ko + 1
27 0
dim x (6,25), u3 (4,9), u4 (5,16), ~5 (6,25)
280
for j = 1 ta 9 : for i = 1 to 4 : read 03 (i, j) : next i,j
290
for j = 1 to 16 : fo:r i = 1 ta 5 : read u4 (i,j) : next i,j
300
for j = 1 to 25:for i, =lto6:read u5 (i,j> : next i,j
3 1 Cl
for j = 1 ta k2
32Cl
for i = 1 to r
3 3 ci
r e a d x (i,j)
34U
next i , j
34 5
goto 305
360
for j = 1 to 9 : for i = 1 to 4 : read ~3 (i,j> : next i,j
370
for j = 1 to 16 : for i = 1 to S : read u4 (i,j) : next i,j
380
for j = 1 to 25 : for i = 1 to 6 : read u5 (i, j) : next i, j
383
Print ch& (19) chr$ (19) chrp( (197).
6
385
rein ¡°choix du plan d'experience"
3 9 0
dim ul (r, k2)
4 00
on ko-2 goto 410, 420, 430
4 10
for j = 1 ta k2 : for i = 1 to r : ul (i,j> = u3 (i,j) : next i,j : goto 440
4 2 0
for j = 1 do k2 : for i = 1 ta r : ul (i,j> = u4 (i,j> : n@Xt i,j : goto 4dO
43c!
for j = 1 to k2 : for i = 1 to r : ul (i,j) = u5 (i,j) : next i,j
44 c!
rem "verification des donn¨¦es"
450
if k2 = 25 ttwn 550
461?
for j = 1 ta k2
46s
p r i n t j ;
47(3
for i = 1 to r
4 3 0
print tab (i*10) x (i,j) :
4 85
next i
49 u
4 9 5
next j
s 00
print "v¨¦rifiez le tableau figurant ¨¤. l' ?cran et notez les erroure 6ventuelJes"
510
print "en rep¨¦rant les numeros de r¨¦p8tition et de traitement correspondants¡±
52C
input "POUR CONTINUER L'ANALYSE COMPOSER C" ; ad
5 3 0
if a6 = "c'" then 710
5 5 0
for j = 1 to 20
555
p r i n t j ;
S60
for i = 1 to r
5 7 0
print tab (i-l 1 10+2!) x (i,j) ;
58Q
next i
5 8 5
5873
next j
590
print ¡°v¨¦rifiez le tableau figurant ¨¤ l¡¯¨¦cran et notez les erreurs eventusl&s"
600
print ¡°en repdrant les num0ros de riSp¨¦tition et de traitement correspondants¡±
510
input "POUR CONTINUER LA VERIFICATION COMPOSER c" ; a+
6 3 0
if ad = r'~" then 640
fi40
for j = 21 to 25
6 4 5
print j ;
6 5 0
for i = 1 to r
660
print tab (i-l) * 10+2) x (i,j) ;
6 6 2
next i
6 6 5
C!I70
next j
6 7 5
680
print "v¨¦rifiez le tablfoau figurant ¨¤ l'¨¦cran et notez les errwrs ¨¦ventuelles"
690
print "sn rep¨¦rant les numiaro ? de repetition et de traitement correspondw?,i::
700
input "POUR CONTINUER L'AN4LYSE COMPOSER c" ; a$
'If! 5
if ag! = "c'. then 710
71C
input "combien de currections sont n¨¦cessaires" SS
715
if SS = 0 then 750
72C!
forua = 1 to ss
730
input "correction
: No de repetition, !\\jo dc traitement, N11.valeur""i,j,r¨ªi,j)
742
next ua
750
rem "calcul des divers totaux"
800
rein " l-total par r8p8titionN
PIO
for i 2 ? ta r
820
for j = 1 to k2
830
tr (i) = tr (i) + x (i,j>
846
next j,i
850
rem "2-total par traitement non ajuste"
i') t;s
dim ,t (k2)
860
for j = 1 to k2
i37U
for i = 1 to r
980
t Ij> = t (j> + x (i,j>
a90
next i,j
900
rem "3'total q¨¦n¨¦rals'
?10
for i = 1 to r
320
for j = 1 to k2
9 3 Il
g = 9 + x (i,j)
Y40
noxt j,i
9 50
rcIm " - 4 total par bloc incomplet"
C] ytj
dim tb (r*ko)
310
for u = 1 to r *ko
9 ? :.l
for i 2 1 to r
380
for j = 1 to k2
990
if ~1 (i,j) = u then tb (u) P tb fu) + x (i,j)
1000
next j, i, u
1010
retn "4- total des blocs incomplets pour chaque traitement¡±
1015
dim bt (k2)
1020
for j = 1 to k2
1030
for i = 1 to r
040
for u = 1 to r * ko
0513
if ul (i,j) i L' then bt (j) = bt (j) + tb (u)
060
next u, i, j
36G
rem "facteur d'ajustamsnt¡°
0
1365
dim IU (k2)
1370
for i = 1 to k2
1380
u)(i) z ko *t(i) - (ko + 1) * ht(i) t g
1390
next i
1400
rcm "variations diverses -1 variation totale"
1410
tc = g?2,&2/&420 far i = 1 t0 r
1430
for j = 1 to k2
1455
vt = vt t x (i,j>t 2
1460
n e x t j , i
14 7 0
vt = vt - tc
1480
rem "2'variatione dUBs aux rGp8titions"
1490
for i = 1 to r
1 SO0
vr = vr c t r
(i)fZ
1510
next i
1520
vr 5: wr/ko/ko - tc
1530
rem ¡°3-variations dues aux traitements"
1532
for i = 1 to kZ
153s
vc - v c t t(i)¡®2
1537
next i
1540
UC = vc/ (ko+l) - tc
1545
FB:I~ ¡°4-variations
dues aux blocs ajustfis"
1550
far i = 1 to k2
1560
vb I vb + w(i) 7 2
157[3
next i
15813
vb = vb/ko*3dko+l) = vb = vb/ko?3/(ko + 1)
1590
rem ¡°erreur intrablocs at facteur d'ajustement (uu
1600
ib = vt - vr - vc - vb
1610
mb = vb/fko+%l)
lG20
im = ib/(ko-1) * (k2-1))
1630
uu = (mb-im)/k2/mb : if uu(= 0 thon uu = 0
1635
dim mc (k2f, i (k2) : tt = 0
1540
rem ¡± totaux a jus t¨¦es par traitement¡±
1645
dim ta (k2), ma (k2)
1650
for i = 1 to k2
3 aio
ta(i) = t (i) + uu * tu(i) : ma (i) t ta (i>/r
1470
next i
1680
f u r i = 1 to k2
1690
tt = tt t ta (i)? 2
17flO
fwxt i
rem ';ci?lctil des wariancss et gain de ps¨¦cisionVf
wt = (tth - tcy(k2-1)
et = ib/(k2-l)/(r-2) : ex = et * (1 + kr-* uu)
~p = 100 ++ ((vb + ib)/(k~&Z-?f/ox) : $p = int (gp*lOO+ ,S)/lOO
rem "impression des rosultats"
print ch& (19) ch& c19j chr+ (747)
print ri
LATTICE "ko" X "ko : print
print tab (5) "moycnfws obs~rv¨¦9s'! tab (40) "moyennes corrig6Ses" : print
dirs mr (k2)
for i = 1 to k2 : mr (i) = t (i)/'c : ncxt i
for i = 1 to k2
print tab (3) "A" i" ' mr (i), teb (44) ma (i)
if i/20 = int (i/20) then 1840
naxt i
Print "notez les rdsultats figurant Gr l'¨¦cran f POUR CONTINUER"
input "LA NDTATIOM COMPOSER c SINON COMPOSER s" ; ad
if ai = "c" then 18313
rem "analyse de variancer'
print : print tab (2s) "ANALYSE DE VARIANCE" print
print "facteur" tab (25) ¡°variante" tab (5~0)~ D . D ?? L . "tab(65)" F ,
caIcu14" : print
19ocj
print iirBpEtitions" i;e"J (25) vr/ko tah (50)ko tab(65) vr/ko/et
? ¡®1¡® 1 fi
print "blocs adj." tm (25) vb/;kZ-1) tab(5C) k2-1 tab(65) vb/(k2-l)/at
1920
print ¡°traitements¡±
tab(2S)urt tab(50) k2-1 tab(65) wt/ex
193c
print ¡°erreur¡± tab(25) ex tab(50) (k2-1) * (r-2)
1 c.1 it, I-J
pri nt : print ¡°gain de pr¨¦cision s "gp - 100" $ 1, : mm = g/r/k2
¡®1 3 4 2
ax = sqr (ex) t p r i n t ¡°coc~ff. do v a r i a t i o n zl¡¯ int((ax/mm) * 104+~~5)/1~0¡± $ ¡±
1945
print "moyenne gBnBr2rle ='l mm : print
195@
p ?int tab(20) reclassement des moyennes ajust¨¦es"
1960
for i = '1 ta k2 : i(i) = i : me(i) = ma (i) : noxt i
1770
for i = 1 to k2-1
1980
for j = 1 to k2-'l
1930
if m3 (j) = no (j+l> thcn 2 ?40
2000
c = ma(j) : d = i(j)
2026
ma(j) = ma(j+l) : i(j) = i(j+l)
2 0 3 0
na(j+l) z c : i(j+l) = d
2 0 4 0
next j
2350
next i
2060
for i = 1 ta k2
10
Print !fAil i(i)" (il ;
next i
dim c3 (a), d3 (8), ct (15), d4 (15), cl (k2- 1), dl (k2-i), ~5 (Lti), d5(24;
on ko-2 9uto 2030, ZW5, 21133
3 c3 fi> = 3.CL........)
3
)tzr,bla de Keule pour 16
. ..r....d3 (8) = 6.22 )
d.d.1.
c4 (1) = bb.,b..*~**L.r.~...,*, ftable do K~ule pour 40 d.d.14
. . . . . . . . . . . . . . d4 (15) = 6,t)2 :)9oto 2115
55 (1) =.**.*o.* 1)tablc ds Kwls pour 120 d.d.1.
\\
..d5(24)= 6.02 '
on ko-2 yoto 212C, 2160, 2592
for i = 1 to k2-?
cl(i) = d(i) -+ sqr (ex/*)
dl(i) = d3(i) + sqr (ox/r)
next i
goto 2193
for i = 1 to k2-?
cl(i) = c4(i) + sqr (ex/r)
dl(i) = d4(i) * sqr (ex/r)
naxt i
yoto 2193
for i = 1 ta k2-1 : cl(i) = c5(i)*eqdox,/r) : dl(i) = d5(i) * sqr(ax/r);noxt i
print : input "waU1ez:-vous le test de Kouls ¨¤ la proba. C.05 OUI ou NGN a$
if a4 = "non" then end
n ZZ l:i=Otj=O
if j + n > k2 thon 2262
print i (j + 1) i
j z j t 1
if j>l thon print i (j) ;
i =i.tl
ifj+n=
k2 + 1 thon 2262
if ma (j+n) . ma(j)<= cl(i) then print II=" i(j+n) i : n = n + 1 : 9oto 2235
if ma (j+n) - me (j),) cl(i) then i = Cl : n = i : print : pt0 2230
print : input '*voulez voue le test de KEULS &I 10 proba. II.01 OUI ou NON";o$
if a$ = "non" then end
print t&b,(lO) "test da Keule PrOba o*[jl"
n = l:i=O:j=O
if j t n)k2 then end
229¡¯3
print i(j + 1) ;
2 2 0 0
j = j + 1
2301
if j >? then print i (j) ;
:,z -5 0 2
i =i+l
2305
ifj+n=
k2 + 1 then end
2 5 1 0
if mz (j t nj - ma(j){= d 1 (i) then print "=" i(j+n) ; : n z n + 1 : gotB.
2302
2320
if ma (j+n) - ma(j)>dl (i) ti&n i = 0 : n = 1 : print : goto 2330
12
- Formulation mathdmatique utilis¨¦e
. totaux :
-
-
r6sultat individuel : x )n ombra de traitements = k2
total lg¨¦n¨¦ral : g
ho
indice de la kattica =Vnbtrait>
total par bloc incomplet : tb
Il existe (r
x ka > valeurs tb
tata1 par rf5p¨¦tition
: tr
Il existe r valeur tr
total par trai teiilent
: t
If existe ko x ka
v a l e u r s t
. pour chaque traitement : total des blocs incomplets 03 il apparait bt
. facteur de correction pour chaque traitement : UI = kot - (ko + 7) bt + g
. kscteurs de variation :
. terme correctif ¨¦c = g
Y (r x kcr2)
. casr4 variation globale vt = c x2 - tc
. carr¨¦ variation r¨¦p¨¦tition vr = jtr2/ko2 - tc
. carr¨¦ variation traitement va = &t2/(ko + 1) - tc
. blocs ajustes carr6is vb = < w2/ko3(ko + 1)
. intrabloc somme des c'arrgs = ib = vt - vr - va - vb
. Variantes :
. blocs ajust¨¦s : vb/(ko2 - 1) Ir) a
. intrablocs vari¨¦s : ib/( 1(0*-l > -43 t
. facteur d'ajustement uu = a - et
liY2-y
si uu<O on consid¨¨re que uu = 0
correctif traitements
?
total ajuste traitement = t +wIo
. variante sur les traitements ajust¨¦s
variation traitements ajust¨¦s = ;5 t ajust4s totaux¡¯/r - tc
variante correspondante
: vz variation trait. ajust&/k2-l
. variante effective de l¡¯erreur : ex = et (1 + ka
:x uu)
Tableau d'analyse de variante
!
!
!
1
variation
!
Variante
!
D.D.L.
!
f
!
!
!
!
!
f
I
!
!
1 R ¨¦pdti tions
I
vr/ke
1
ka
! ur/ko/et
t
1
!
!
!
1
?
l
!
!
I
I Olocs ajust¨¦s
!
vb/(ko2 - 1)
I ko2 - 1
!vb/(ko2-l)/et !
!
!
!
!
!
!
!
!
!-
I
v a r i a t i o n t r a i t
! Traitements ajus@IajustBs/'ko2 - 1
I ko2 . 1
1
variance/ex
!
!
I
!
!
1
I
!
!
t
! Erreur ajust¨¦e
!
63 x
!(k 24)(k 4cll) !
1
!
I
I
I
!
13
- jode d'emploi du proqramme "Lattice 543"
load "lattice 543" return
loading >
rsady
) 3 mn 05 sec
run return
Voulez vous conserver les donn¨¦es OUI ou NON ? oui return
introduisez les donn¨¦es en DATA ¨¤ partir de 5000 en respectant l¡¯ordre ;
du tableau do r$sultats soit :
!
d¡¯abord tout¨¦s las rbp¨¦titions du
2 etc..,
traitement '1 puis celles du Eraitem ?nti
POUR CONTINUER L"ANALYSE COMPOSER ru-l 260
1
I
break in 180
1
f
ready
I
Tableau de r¨¦sultats
I
1
LATTICE
J
!
!
!
I
II
l
I
3 x 3
1
r
IX1
!
IV
f
!
-!
.
f
!
I
!
! Traitements , 1 I
?-!
2299
*
23,5
;
16,9
]
24,2
f
!
!
f
1'
2
34,4
;
3284
I
!
i
!
!
I - ?
?
-3
1
2697
!
3696
I
!
(1
I
!
!
! 3 !
2896
1
31,4
1
26,2
f
27,9
;
f
!
!
II
-!
2891 -!
28,6
!
I
!
!
!
! 4 !
[
I
25,1
1
3397
,
f
!
-!
i
I
!
II
I- ?
!
! 5!
3397
i
30,7
;
2693
!
3091
1
I
!-!
29,0
i
;
34,3
I
!
I
!
'9
I 6
32,0
I
1
.
31¡®0
f
I
!-
If
? - I
?
!
1
!!-
! 7 f
2591
,
30,2
1
24¡®9
f
29,4
;
I
!-!
If
-1
I
1
!
!
! 8,
31¡®6
f
33,l
i
1
-*
27,8
1
3119
I
!
!
-1
II
r
!
Pf
2996
1
33,5
;
29,4
f
29,7
f
\\
Selon Cochran et Cox.
1
"Experimental designs" Page 304
1
I
SO¨¹O data 22.9, 23.5, 16.9, 24.2, 34.4,......., 27.9 rtiturn
I
(3)
t
5001 data 28.1, 20.6, 25.1, 33.7,....,.., 25.1, 30.2 return
I
3hi.s
t
*un 260 return
!
(4)
Dutil est le nombre de traitement ¨¤ tester ? 9 return
I
(5)
1
22.9
23.5
16.9
24.2
I
2
34.4
32.4
26.7
36.6
f
3
28.6
31.4
26.2
27.9
I
4
28.1
28.6
25.1
33.7
I
5
33.7
30.7
26.3
30.1
t
6
29
32
34.3
31
1
7
25.1
30.2
24.9
29.4
I
a
3 1 . 6
3 3 .l
2 7 . 8
31.9
1
9
2 9 . 6
3 3 . 5
29.4
29.7
I
14
!
No R6f.
f
v¨¦rifiez le tableau figurant ¨¤ lIecran et notez 10s erreurs kventuclles
!
un rcperant les num¨¦ros de r&p¨¦tition et de traitements C:orrespondants
!
!
POIJR CONTINUER L'ANALYSE COMPOSER c ? c return
!
Cwnbicn de correction sont n¨¦cessaires ? 0 roturn
!
f
LATTICE 3 x 3
1
I
!
moyonncs
observ¨¦es
moyennes corrig¨¦es
I
A 1
21 .a75
21.875
I
A 2
!
32.525
32.525
!
Cl 3
23.525
28.525
!
P, 4
!
20.875
28.875
!
A 5
30.2
30.2
!
!
A 6
31.575
31.575
!
A 7
27.4
27.4
I
I
A 8
31.1
31.1
!
A 9
30.55
30.55
!
I
notez les r8sultats figurant ¨¤ 1"Bcran. POUR CONTINUER LA NOTATION
t
!
COMPOSER c SIPJON COMPOSER s 5' s rsturn
I
(9)
!
ANALYSE DE VARXANCE
!
I
r
facteur
variante
D.D.L.
F c a l c u l ¨¦
!
1
r-ipt! ti ti ons
34.4669266
3
5.37297662
f
Blocs adj.
3.40611107
8
,530971482
!
traitements
!
40.2911053
8
6.28089557
!
erreurs
6.4 148663
i6
I!
goin d e p r ¨¦ c i s i o n = - 1 5 . 6 3 $
I
cocff. d e v a r i a t i o n =
8 . 6 8 $
I
!
moy g¨¦n¨¦rale
= 29.1805556
!
Classement des moyennes a.iustdes
!!
Al< A? (A3 (A4<AS (A9( A04 MC: A2(
I
!
Voulez vous le test de Keuls ¨¤ la proba 0.05 OUI ou NON ? oui return
t
(9)
Voulez vous le test de Keul ¨¤ la proba 0.01 OLII ou NON ? non return
I
(10)
ready
i
Variantes dans le mode dtemploi du proqramme
1
I
ru17 return
f
(1)
voulez vous conserver les donnbes OUI ou NON ¡®1 non return
1
(2)
quo1 est le nombre de traitements ¨¤ tester ? 9 return
I
I
(31
K¨¦sultat de la rbp¨¦tition no 1 par le traitement no 1 ? 2 2 . 9 r e t u r n
I
(4)
,I
l
2
1 ? 23.5 return
!
;t
-3
n
-> ,..
15
1 NO H¨¦f
Resultat de la rep¨¦tition na 3 par le traitement no 9 ? 29.4 return
!
.
RQsultat de la r¨¦petition n" 4 par le traitement no 9 ? 29.7 raturn
1
I
Verifiez le tableau figurant ¨¤ 1¡¯6cran et notez les erreurs Qventuelles I
!
i.;n rep¨¦rant les numdros de rep¨¦tition et de traitements correspondants I
POUR CONTINUER L'ANALYSE COMPOSER c ? c
(5)
!Zambien de corrections sont n6cessaires ? 1 return
i
!
(6)
c 3rr ec ti on
: No d e r¨¦pAtition,
oJo d e t r a i t e m e n t , Nil. valeur ? 2,2,32,4 I
(7)
return !I
I
Les rQsultats apparaissent ensuite
B l¡¯ ?cran en tenant compte de la correction effectu¨¦e f
Remarques :
1) ¨¤ 1'5tape (6) de la variation, on s'¨¦tait aper ?u en v¨¦rifiant le ta- I
bleau de r¨¦sultats que l'on avait fait une erreur sur la 2Bme r¨¦pdti-.f
tion du traitement 2 c'est Pourquoi on a r¨¦pondu 1 B la question po&td
et ensuite on a introduit la correction ( r e p ) Nll.valeur(la virgule !
2 , 2 , 3 2 . 4
I
separe 2 donn¨¦es)
( t r a i t . )
!
!
!
2) Dans le cas d'une lattice 5 x 5, l'ecran ne pouvait pas contenir tout
le tableau de r¨¦sultats pour le v¨¦rifier, i l a p p a r a i t l e s 2 0 p r e m i e r s
traitements,
puis dans un deuxi¨¦me temps les 5 autres. Il en va de
I
m&me pour les moyennes, o¨´ il est alors n8cessaire de composer c pour1
voir apparaitre les 5 derni¨¨res moyennes.
3) Cans l¡¯exemple traite le coefft de correction
est negatif ce qui
I!
Explique ltQgalite des moyennes observdes et corrigees.
1
4) Si l'essai est non significatif on rgpondra non a la question No Y et:
l'appareil affichera ready.
5) En (7) l¡¯ordinateur calcule pendant 1 B 4 mn environ selon Ie type de! !
lattice (3 x 3 - 4 x 4 - 5 x 5).
I
(g) Le test de Keuls fait apparaitre certains classements inutiles,
I
il ne faudra pas en tenir compte dans l¡¯interpr¨¦tation.
Test de Keuls proba 0.05
1
7=3=4=5=9=8=6=2