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