Les contenus sont © P. Calmant et E. Depiereux - 2004; G. Vincke B. De Hertogh et E. Depiereux 2008.
Imprimé le
8/1/2025
Comme vu précédemment, le protocole d'une expérience peut être modélisé en une équation mathématique.
Ce chapitre nous permet de déterminer, par un test statistique de Fisher-Snedecor, si chacun des composants de variance de l'équation a une influence significative sur le paramètre étudié.
Méthodologie :
Énoncé de l'expérience:
Pour fixer les quotas laitiers de la Région wallonne, l'expérimentateur prospecte l'Ardenne, le Condroz et la Hesbaye. Dans chaque région, il recueille les statistiques de production de 4 vaches prises au hasard, en été et en hiver. Sur chacune de ces vaches, il réalise 3 traites pour mesurer la quantité de lait. |
été |
hiver |
||||||
traite 1 |
traite 2 |
traite 3 |
traite 1 |
traite 2 |
traite 3 |
||
Ardenne
|
vache 1 | qtité lait |
|||||
vache 2 | |||||||
vache 3 | |||||||
vache 4 | |||||||
Condroz |
vache 1 | ||||||
vache 2 | |||||||
vache 3 | |||||||
vache 4 | |||||||
Hesbaye |
vache 1 | ||||||
vache 2 | |||||||
vache 3 | |||||||
vache 4 |
Critères:
Région: ai où i = 3
- critère fixe avec 3 niveaux
- croisé au critère SaisonVache: B(i)j où j = 4
- critère aléatoire avec 4 niveaux
- hiérarchisé au critère Région
- croisé au critère SaisonSaison: ck où k = 2
- critère fixe à 2 niveaux
- croisé au critère Région
- croisé au critère VacheNombre de réplicats: l = 3
Équation:
x(ijk) l = mx + ai + B(i)j + ck + acik + Bc(i)jk + E(ijk)l
la production laitière serait donc influençable par:
ne pas oublier de tenir compte des erreurs de mesures E(ijk)l
(important quand la mesure est répétée)
Application de la règle 1 :
δ2 : pour les critères fixes Si l'indice :
|
i = 3 | j = 4 | k = 2 | l = 3 | ||
δ2a | ai | 0 | 4 | 2 | 3 |
S2B | B(i)j | 1 | 1 | 2 | 3 |
δ2c | ck | 3 | 4 | 0 | 3 |
δ2ac | acik | 0 | 4 | 0 | 3 |
S2Bc | Bc(i)jk | 1 | 1 | 0 | 3 |
S2 |
E(ijk)l | 1 | 1 | 1 | 1 |
Application de la règle 2 :
Considérer uniquement les sources de variabilité |
exemple:
|
i = 3
|
j = 4
|
k = 2
|
l = 3
|
||||||||
δ2a
|
ai
|
0
|
4
|
2
|
3
|
δ2a
|
S2B
|
δ2ac
|
S2Bc
|
S2
|
|
S2B
|
B(i)j
|
1
|
1
|
2
|
3
|
S2B
|
S2Bc
|
S2
|
|||
δ2c
|
ck
|
3
|
4
|
0
|
3
|
δ2c
|
δ2ac
|
S2Bc
|
S2
|
||
δ2ac
|
acik
|
0
|
4
|
0
|
3
|
δ2ac
|
S2Bc
|
S2
|
|||
S2Bc
|
Bc(i)jk
|
1
|
1
|
0
|
3
|
S2Bc
|
S2
|
||||
S2
|
E(ijk)l
|
1
|
1
|
1
|
1
|
S2
|
Espérances des carrés moyens de chaque composant:
Application de la règle 3 :
|
i = 3
|
j = 4
|
k = 2
|
l = 3
|
||||||||
δ2a
|
ai
|
0
|
4
|
2
|
3
|
24δ2a
|
6S2B
|
0δ2ac
|
0S2Bc
|
1S2
|
|
S2B
|
B(i)j
|
1
|
1
|
2
|
3
|
6S2B
|
0S2Bc
|
1S2
|
|||
δ2c
|
ck
|
3
|
4
|
0
|
3
|
36δ2c
|
0δ2ac
|
3S2Bc
|
1S2
|
||
δ2ac
|
acik
|
0
|
4
|
0
|
3
|
12δ2ac
|
3S2Bc
|
1S2
|
|||
S2Bc
|
Bc(i)jk
|
1
|
1
|
0
|
3
|
3S2Bc
|
1S2
|
||||
S2
|
E(ijk)l
|
1
|
1
|
1
|
1
|
1S2
|
Rappelons que l'intérêt de cette démarche, est de déterminer quels composants de l'équation initiale sont pertinents pour expliquer la production laitière.
Chaque composant devra donc faire l'objet d'un test d'hypothèses
et donc il faut pour chacun de ces composants calculer un F observé.
x(ijk)l = mx + ai + B(i)j + ck + acik + Bc(i)jk + E(ijk)l
Application de la règle 4 pour déterminer les degrés de liberté de chaque test d'hypothèses
Effectuer le produit de la valeur maximale de tous
les indices représentés en tête de
|
i=3
|
j=4
|
k=2
|
l=3
|
dl
|
|||||||||
test 1
|
δ2a
|
ai
|
0
|
4
|
2
|
3
|
24δ2a
|
6S2B
|
0δ2ac
|
0S2Bc
|
1S2
|
2
|
|
test 2
|
S2B
|
B(i)j
|
1
|
1
|
2
|
3
|
6S2B
|
0S2Bc
|
1S2
|
9
|
|||
test 3
|
δ2c
|
ck
|
3
|
4
|
0
|
3
|
36δ2c
|
0δ2ac
|
3S2Bc
|
1S2
|
1
|
||
test 4
|
δ2ac
|
acik
|
0
|
4
|
0
|
3
|
12δ2ac
|
3S2Bc
|
1S2
|
2
|
|||
test 5
|
S2Bc
|
Bc(i)jk
|
1
|
1
|
0
|
3
|
3S2Bc
|
1S2
|
9
|
||||
S2
|
E(ijk)l
|
1
|
1
|
1
|
1
|
1S2
|
48
|
H0: δ2a = 0
H1 : δ2a ≠ 0si δ2a = 0, ai = 6S2B + 1S2
donc, la source de variabilité ai a la même équation que la source de variabilité B(i)j
donc, F observé = E ( CMai ) / E ( CMB(i)j)
H0: S2B = 0
H1 : S2B ≠ 0
si S2B = 0, B(i)j = 1S2
donc, la source de variabilité B(i)j a la même équation que le Carré Moyen résiduel
donc, F observé = E ( CMB(i)j ) / E ( CMR )
H0: δ2c = 0
H1 : δ2c ≠ 0
si δ2c = 0, ck = 3S2Bc + 1S2
donc, la source de variabilité ck a la même équation que la source de variabilité Bc(i)jk
donc, F observé = E ( CMck ) / E ( CMBc(i)jk )
H0: δ2ac = 0
H1 : δ2ac ≠ 0
si δ2ac = 0, acik = 3S2Bc + 1S2
donc, la source de variabilité acik a la même équation que la source de variabilité Bc(i)jk
donc, F observé = E ( CMacik ) / E ( CMBc(i)jk )si le test est significatif, alors δ2ac ≠ 0
alors, il faut maintenir acik dans l'équation finale puisque l'interaction région/saison influence la production laitière de manière significative
si le test est non significatif, alors δ2ac = 0
alors, nous pouvons supprimer acik de l'équation finale puisque l'interaction région/saison n'influence pas la production laitière
H0: S2Bc = 0
H1 : S2Bc ≠ 0
si S2Bc = 0, Bc(i)jk =1S2
donc, la source de variabilité Bc(i)jk a la même équation que le Carré Moyen résiduel
donc, F observé = E ( CMBc(i)jk )/ E ( CMR )