RÉPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

UNIVERITE FERHAT ABBAS SETIF

FACULTE DES SCIENCES DE LA NATURE ET DE LA VIE

DEPARTEMENT DES SCIENCES AGRONOMIQUES

02

COURS DE BIOSTATISTIQUE

Mr. Rouabhi. A

Année universitaire 2014 - 2015

 Le plan factoriel
Plan d’expérience factoriel
• une expérience menée sur le rendement d’une plante, englobant deux facteurs, mesurés à deux niveaux,
par exemple deux types engrais à deux dosages différents : est une expérience factorielle 2 x 2 ou 2 2.
• Les traitements sont constitués des quatre combinaisons possibles de chacun des deux facteurs, aux
deux niveaux considérés.

Facteur A (type Facteur B (dose

d’engrais) d’engrais)

Azote Dose 1

Dose 2
Phosphore

Numéro du traitement Combinaison des traitements

Engrais (facteur A) Doses (facteur B)
1 Azote Élevée
2 Azote Faible
3 phosphore Élevée
2 phosphore Faible
Pour une expérience factorielle comprenant:
5 niveaux du facteur A et 3 niveaux du facteur B, le nombre total de traitements sera 5 x 3
NB: chaque traitement doit avoir des répétitions.

Anova à 2 facteurs
*C'est l'extension d'une ANOVA à un (1) facteur de classification. On regarde l'effet de plusieurs facteurs
*Ex: une ANOVA à deux facteurs de classification permet de déterminer l'effet de l'addition d’un engrais
(facteur A) et à différentes doses (facteur B) sur le rendement.
*Contrairement à l'ANOVA à un facteur, ici on doit tester plusieurs H0
• Pour deux facteurs A et B à 2 niveaux chacun (a1 et a2, b1 et b2, respectivement), donc
• une hypothèse nulle concernant :
• 1. l'effet du facteur A (H0 : ma1 = ma2 ; H1 : ma1 ≠ ma2)
• 2. l'effet du facteur B (H0 : mb1 = mb2 ; H1 : mb1 ≠ mb2)
• 3. l'effet d'interaction AxB (H0 : pas d’interaction ; H1 : interaction)
L'effet d'interaction traduit simplement le fait que l'effet d'un facteur dépend de l'autre facteur.
C’est-à-dire que l'effet d'un facteur n'est pas le même, il change selon les niveaux de l'autre facteur.

Modèle sans interaction

Dose (Facteur B)
Engrais (Facteur traitement Faible (b1) Élevée (b2)
A)
Azote (a1) 20 30
Phosphore (a2) 40 52

Modèle avec interaction

(Facteur B)
Engrais (Facteur traitement Faible (b1) Élevée (b2)
A)
Azote (a1) 20 40
Phosphore (a2) 50 12

Plus d’exemples
 Variable de réponse: rythme cardiaque au repos
Facteur A : dose de médicament
Facteur B : catégorie d’âge

Interaction

Anova à 2 facteur REP

Nous allons étudier les différentes étapes de la procédure d’analyse de la variance d’une expérience à 2 facteurs
sur le rendement d’une céréale (q/ha), avec deux sorte d’engrais (Facteur A) et trois niveaux de doses (faible,
moyenne, élevée) (facteur B), définis dans un PER, à trois répétitions.

Type d’engrais (facteur A)

Azote (a1) Phosphore (a2) Total rép
a1b1 a1b2 a1b3 a2b1 a2b2 a2b3
Dose Rép1 46.5 49.5 127.7 49.3 65.5 67.9 406.4
d’engrais
Rép2 55.9 59.5 134.1 53.2 65.0 112.7 480.4
(facteur B)
Rép3 78.7 78.7 137.1 65.3 74.0 129.0 562.8

Total traitement (Ti) 181.1 187.7 398.9 167.8 204.5 309.6 G=1449.6

Avant de dresser le tableau d’Anova il faut donner les paramètres:

« r » =3 (le nombre de répétitions)
« a » =2 (le nombre de niveaux du (facteur A) azote et phosphore),
« b » = 3 le nombre de niveaux du facteur B (doses). (faible, moyenne, élevée)

Source de variation Degré de liberté (dll) Somme des carrés (SC) Carré Moyen CM=SC/ddl F calculé

Total rab – 1 SCT

Traitement ab – 1 SCF CMF CMF/CMR

 A a-1 SCA CMA CMA/CMR

B b-1 SCB CMB CMB/CMR

AB (a-1)(b-1) SCAB CMAB CMAB/CMR

Erreur (rab - ab) SCR CMR

Calculer les totaux des traitements (Tij), les totaux des doses (Rk), et le total général (G)
TC= G²/rab=(1449.6)²/3*2*3=116741.1

SCT = -TC = [(46.5)²+ (55.9)²+…..(129)²]-116741.1= 17479.1

SCF = -TC = [[(181.1)²+(187.7)²+….(309.6)²]/3]-116741.1= 14251.8

Tableau de l’Analyse de variance préliminaire

Source de variation Degré de liberté Somme des carrés Carré moyen F calculé F tabulaire

Total 17 SCT=17479.1 5%

Traitement 5 SCF=14251.8 2850.3 10.5** 3.3

Erreur 12 SCR=3227.3 268.9

Le traitement affiche un effet significatif; s’agit –il du facteur A ou B ou des deux.

Consolider un nouveau tableau à 2 entrées

Facteur A
Azote Phosphore Total (Bj)
Facteur B Dose faible 181.1 167.8 348.9
Dose Moyenne 187.7 204.5 392.2
Dose élevée 398.9 309.6 708.5
Total (Ai) 767.7 681.9 G=1449.6
Le traitement affiche un effet significatif; s’agit –il du facteur A ou B ou des deux
On va calculer les trois composantes factorielles de la somme des carrés des traitements:

SCA = -TC = [[(767.7)²+(681.9)²]/3*3]-116741.1= 408.9

SCB = -TC = [[(348.9)²+(392.2)²+(708.5)²]/3*2]-116741.1= 12846.2
SCAB = SCF-SCA-SCB = 14251.8- 408.9-12846.2= 996.6

Source de Degré de liberté (dll) Somme des carrés (SC) Carré Moyen CM=SC/ddl F calculé
variation

Total 17 (rab – 1) SCT=17479.1

Traitement 5 (ab – 1) SCF=14251.8 CMF=2850.3 CMF/CMR=10.5*

A 1 (a-1) SCA=408.9 CMA=408.9 CMA/CMR=1.5ns

B 2 (b-1) SCB=12846.2 CMB=6423.1 CMB/CMR=23.8*

AB 2 (a-1) (b-1) SCAB=996.6 CMAB=498.3 CMAB/CMR=1.8ns

Erreur 12 (rab-ab) SCR=3227.2 CMR=268.9

Révision du dispositif factoriel

Plan factoriel
• Dispositif à un seul facteur
• Dispositif à un deux facteurs (A et B)
• Facteur A : plusieurs niveaux
• Facteur b : plusieurs niveaux
*Facteur A (fertilisant à 2 niveaux)
*Facteur b: (dose à 3 niveaux)

Répétition (1er
Dose élevée
groupe)

répétitions
Azote Dose modérée
(2ème groupe)

Répétitions
Dose faible
(3ème groupe)
Fact AxB
Répétitions
Dose élevée
(4ème groupe

Répétitions
phosphore Dose modérée
(5ème groupe)

Répétitions
Dose faible
(6ème groupe)
Disposition sur champs
Azote dose Phosphore dose Azote dose Azote dose Phosphore dose Phosphore dose
modérée modérée élevée faible élevée faible
Azote dose Phosphore dose Phosphore dose Azote dose Phosphore dose Azote dose
modérée modérée élevée élevée faible faible
Phosphore dose Azote dose Phosphore dose Phosphore dose Azote dose Azote dose
modérée modérée faible élevée élevée faible

A1b3(127.7) A2b3(67.9) A1b1(46.5) A1b2(49.5) A2b1(49.3) A2b2(65.5)

A1b3(134.1) A2b3(112.7) A2b1(53.2) A1b1(55.9) A2b2(65) A1b2(59.5)
A2b3(129) A1b3(137.1) A2b2(74) A2b1(65.3) A1b1(78.7) A1b2(78.7)

Type d’engrais (facteur A)

Azote (a1) Phosphore (a2) Total
rép
a1b1 a1b2 a1b3 a2b1 a2b2 a2b3
Dose Rép1 46.5 49.5 127.7 49.3 65.5 67.9 406.4
d’engrais
Rép2 55.9 59.5 134.1 53.2 65.0 112.7 480.4
(facteur B)
Rép3 78.7 78.7 137.1 65.3 74.0 129.0 562.8

Total traitement (Ti) 181.1 187.7 398.9 167.8 204.5 309.6 G=1449.6

Avant de dresser le tableau d’Anova il faut donner les paramètres:

« r » =3 (le nombre de répétitions)
« a » =2 (le nombre de niveaux du (facteur A) azote et phosphore),
« b » = 3 le nombre de niveaux du facteur B (doses). (faible, moyenne, élevée)

Source de variation Degré de liberté (dll) Somme des carrés (SC) Carré Moyen CM=SC/ddl F calculé

Total rab - 1 SCT

Traitement ab - 1 SCF CMF CMF/CMR

A a-1 SCA CMA CMA/CMR

B b-1 SCB CMB CMB/CMR

AB (a-1)(b-1) SCAB CMAB CMAB/CMR

Erreur (rab - ab) SCR CMR

Calculer les totaux des traitements (Tij), les totaux des doses (Rk), et le total général (G)
TC= G²/rab= (1449.6)²/3*2*3=116741.1
SCT = -TC = [(46.5)²+(55.9)²+…..(129)²]-116741.1= 17479.1

SCF = -TC= [[(181.1)²+(187.7)²+….(309.6)²]/3]-116741.1= 14251.8

Tableau de l’Analyse de variance préliminaire

Source de variation Degré de liberté Somme des carrés Carré moyen F calculé F tabulaire

Total 17 SCT=17479.1 5%

Traitement 5 SCF=14251.8 2850.3 10.5** 3.3

Erreur 12 SCR=3227.3 268.9

Le traitement affiche un effet significatif; s’agit –il du facteur A ou B ou des deux
Consolider un nouveau tableau à 2 entrées :

Facteur A
Azote Phosphore Total (Bj)
Facteur B Dose faible 181.1 167.8 348.9
Dose Moyenne 187.7 204.5 392.2
Dose élevée 398.9 309.6 708.5
Total (Ai) 767.7 681.9 G=1449.6
Le traitement affiche un effet significatif; s’agit –il du facteur A ou B ou des deux
On va calculer les trois composantes factorielles de la somme des carrés des traitements:

SCA = -TC= [[(767.7)²+(681.9)²]/3*3]-116741.1= 408.9

SCB = -TC= [[(348.9)²+(392.2)²+(708.5)²]/3*2]-116741.1= 12846.2

SCAB = SCF-SCA-SCB = 14251.8- 408.9-12846.2= 996.6

Source de variation Degré de liberté (dll) Somme des carrés (SC) Carré Moyen CM=SC/ddl F calculé

Total 17 (rab – 1) SCT=17479.1

Traitement 5 (ab – 1) SCF=14251.8 CMF=2850.3 CMF/CMR=10.5*

A 1 (a-1) SCA=408.9 CMA=408.9 CMA/CMR=1.5ns

 B 2 (b-1) SCB=12846.2 CMB=6423.1 CMB/CMR=23.8*

AB 2 (a-1)(b-1) SCAB=996.6 CMAB=498.3 CMAB/CMR=1.8ns

Erreur 12 (rab-ab) SCR=3227.2 CMR=268.9

Exercice 3
on veut tester l’effet de l’altitude (fact-A) et l’effet variétal (Fact-B) sur la productivité du pommier. À cet effet,
nous avons testé à deux niveaux d’altitude (haute A1 et moyenne altitude A2) trois variétés de pommier (V1,
V2 et V3).
Les rendements moyens obtenus sont données en q/ha.
A1V1 (25.9) A2V2 (26.8) A2V1 (30.7) A1V1 (27.8) A1V2 (29.3) A1V3 (31.8)

A2V1 (31.2) A2V3 (19.1) A1V3 (35.4) A2V3 (22.8) A2V2 (27.2) A2V1 (33.8)

A1V2 (26.3) A1V1 (24.2) A2V2 (23.8) A2V3 (21.1) A1V3 (42.6) A1V2 (30.6)

Tests des effets inter-sujets

Variable dépendante:VAR00003

Source Somme des carrés de

type III ddl Moyenne des carrés D Sig.

Altitude 77,709 1 77,709 9,714 ,009

Variété 9,458 2 4,729 ,591 ,569

Alt* Var 351,898 2 175,949 21,994 ,000

Erreur 96,000 12 8,000

Total 535,064 17

a. R deux = ,821 (R deux ajusté = ,746)

Dispositif factoriel en blocs

Nous allons décrire les différentes étapes de la procédure d’analyse de la variance d’une expérience à 2 facteurs
sur le rendement, avec deux types d’engrais (Facteur A) et trois niveaux de doses (faible, moyenne, élevée)
(facteur B), définis dans un PEBAC, à trois répétitions.
 Dose d’engrais Type d’engrais (facteur A)
Azote (a1) Phosphore (a2)
Faible (b1) a1b1 a2b1
Moyenne (b2) a1b2 a2b2
Élevée (b3) a1b3 a2b3

A1B1(46,5) A1B3(127,7) A1B2(49,5) A2B1(49,3) A2B3(67,9) A2B2(65,5)

A1B3(134,1) A2B1 (53,2) A1B1(55,9) A2B3(112,7) A1B2(59,5) A2B2 (65)
A2B3 (129) A2B2 (74) A1B3(137,1) A1B2 (78,7) A1B1(78,7) A2B1(65,3)

Combinaison des traitements Rendement obtenu après récolte Total traitement (Tij)
Bloc1 (rép1) Bloc2 (rép2) Bloc2 (rép2)
a1b1 46.5 55.9 78.7 181.1
a1b2 49.5 59.5 78.7 187.7
a1b3 127.7 134.1 137.1 398.9
a2b1 49.3 53.2 65.3 167.8
a2b2 65.5 65.0 74.0 204.5
a2b3 67.9 112.7 129.0 309.6
total répétition (Rk) 406.4 480.4 562.8 G=1449.6

Avant de dresser le tableau d’Anova il faut donner les paramètres:

« r » =3 (le nombre de répétitions (blocs))
« a » =2 (le nombre de niveaux du (facteur A) azote et phosphore),
« b » = 3 le nombre de niveaux du facteur B (doses). (faible, moyenne, élevée).

Source de variation Degré de liberté (dll) Somme des carrés (SC) Carré Moyen CM=SC/ddl F calculé

Total rab - 1 SCT

Bloc r–1 SCB CMB CMB/CMR

Traitement ab - 1 SCF CMF CMF/CMR

A a-1 SCA CMA CMA/CMR

B b-1 SCB CMB CMB/CMR

AB (a-1)(b-1) SCAB CMAB CMAB/CMR

Erreur (r - 1)(ab - 1) SCR CMR

Calculer les totaux des traitements (Tij), les totaux des répétitions (Rk), et le total général (G)
TC= G²/rab=(1449.6)²/3*2*3=116741.1

SCT = -TC= [(46.5)²+(55.9)²+…..(129)²]-116741.1= 17479.1

SCB = -TC = [[(406.4)²+(480.4)²+(562.8)²]/(2)*(3)]-116741.1= 2040.3

SCF = -TC= [[(181.1)²+(187.7)²+….(309.6)²]/3]-116741.1= 14251.8

SCR =SCT-SCB-SCF = 17479.1-2040.3-14251.8 = 1186.8

Tableau de l’Analyse de variance préliminaire

Source de variation Degré de liberté Somme des carrés Carré moyen F calculé F tabulaire

Total 17 17479.1 5%

Bloc 2 2040.3 1020.1 8.5* 4.1

Traitement 5 14251.8 2850.3 24.0** 3.3

Erreur 10 1186.8 118.6

*Construire le tableau à double entrée des totaux (facteur A) X (facteur B), avec le calcul des totaux du (facteur
A) et les totaux du (facteur B).

Dose (Facteur B)
Engrais traitement Faible (b1) Moyenne (b2) Élevée (b3) Total (Ai)
(Facteur A)
Azote (a1) 181.1 187.7 398.9 767.7
Phosphore (a2) 167.8 204.5 309.6 681.9
Total (Bi) 348.9 392.2 708.5 G=1449.6
Calculer les trois composantes factorielles de la somme des carrés des traitements:

SCA = -TC= [[(767.7)²+(681.9)²]/3*3]-116741.1= 408.9

SCB = -TC= [[(348.9)²+(392.2)²+(708.5)²]/3*2]-116741.1= 12846.2

SCAB = SCF-SCA-SCB = 14251.8- 408.9-12846.2= 996.6
Source de variation Degré de liberté (dll) Somme des carrés(SC) Carré Moyen CM=SC/ddl F calculé

Total 17 (rab – 1) SCT=17479.1

Bloc 2 (r – 1) SCB=2040.3 CMB=1020.1 CMB/CMR=8.6*

Traitement 5 (ab – 1) SCF=14251.8 CMF=2850.3 CMF/CMR=24.0*

A 1 (a-1) SCA=408.9 CMA=408.9 CMA/CMR=3.4

B 2 (b-1) SCB=12846.2 CMB=6423.1 CMB/CMR=54.1*

AB 2 (a-1)(b-1) SCAB=996.6 CMAB=498.3 CMAB/CMR=4.2*

Erreur 10 (r - 1)(ab 1) SCR=118.6 CMR=118.6

A1B1 A1B3 A2B3 A1B2 A1B3 A1B2

A1B2 A2B1 A1B1 A2B1 A2B2 A1B1
A2B2 A2B3 A1B3 A2B2 A2B3 A2B3
A1B3 A1B2 A2B1 A1B1 A2B1 A2B2

A1B1 A2B2 A2B3 A1B2 A1B3 A2B1

A1B2 A2B1 A1B1 A1B3 A2B2 A2B3
A2B2 A2B3 A1B3 A1B2 A2B1 A1B1
A1B3 A1B2 A2B1 A1B1 A2B3 A2B2

Fact A: Variété
Fact B: Fertilisant

A1B1 A1B2 A1B3 A2B2 A2B3 A2B1

A1B2 A1B1 A1B1 A2B3 A2B2 A2B3
 A1B2 A1B3 A1B3 A2B2 A2B1 A2B1
A1B3 A1B2 A1B1 A2B1 A2B3 A2B2

Var A1 Var A2
A: Variété
B: Fertilisant
C: Dose d’irrigation

A1B1C1 A1B2C1 A1B3C1 A2B2C2 A2B3C2 A2B1C2

A1B2C1 A1B3C1 A1B1C1 A2B3C2 A2B1C2 A2B2C2
A1B1C2 A1B2C2 A1B3C2 A2B2C1 A2B1C1 A2B1C1
A1B3C2 A1B2C2 A1B1C2 A2B3C1 A2B3C1 A2B2C1

A1B1C1 A1B2C1 A1B3C1 A2B2C2 A2B3C2 A2B1C2

A1B2C1 A1B3C1 A1B1C1 A2B3C2 A2B1C2 A2B2C2

Var A1 A1B1C2 A1B2C2 A1B3C2 A2B2C1 A2B1C1 A2B1C1 Var A2

A1B3C2 A1B2C2 A1B1C2 A2B3C1 A2B3C1 A2B2C1

A1B1C1 A1B2C1 A1B3C1 A2B2C2 A2B3C2 A2B1C2

IRR:C1 A1B2C1 A1B3C1 A1B1C1 A2B3C2 A2B1C2 A2B2C2
A1B1C2 A1B2C2 A1B3C2 A2B2C1 A2B1C1 A2B1C1
IRR:C2
A1B3C2 A1B2C2 A1B1C2 A2B3C1 A2B3C1 A2B2C1

Traitement (A B C D)

B C D A
C D A B
A B C D
D A B C