556
Le logiciel R - Effectuer des analyses statistiques
Anova Table (Type II tests)
Response: rdt
Sum Sq Df F value
Pr(>F)
region
327.19 3 112.1810 0.00000000000002955 ***
engrais
0.89 2
0.4571
0.6385
region:engrais 99.56 6 17.0667 0.00000013594213061 ***
Residuals
23.33 24
--Signif. codes: 0 ‘ *** ’ 0.001 ‘ ** ’ 0.01 ‘ * ’ 0.05 ‘.’ 0.1 ‘ ’ 1
Remarque
La formule region*engrais, utilis´ ee dans aov() et lm(), correspond
en fait `
a la formule region+engrais+region:engrais, c’est-` a-dire
le facteur r´ egion, le facteur engrais et l’interaction entre ces deux
facteurs.
La valeur-p associ´ ee au test de la pr´ esence d’interaction est significative.
Cela entraˆ ıne par exemple que l’effet de l’engrais sur le rendement peut
ˆ etre diff´ erent selon les r´ egions.
Attention
Nous consid´ erons qu’il n’y a aucun effet d’interaction si la valeur-p
associ´ ee est sup´ erieure ` a 5 %. Dans ce cas, nous effectuons une analyse
de la variance sans terme d’interaction, ce qui rend l’effet principal
plus facile ` a interpr´ eter. En pr´ esence d’interaction, il ne faut pas interpr´ eter les tests des effets principaux qui sont pr´ esent´ es dans les
sorties de la table d’analyse de la variance.
Test des effets conditionnels en pr´ esence d’interaction : On souhaite par
exemple savoir s’il existe un effet engrais dans la r´ egion 1. Pour cela, on
peut utiliser la fonction subset() qui permet de n’utiliser que les donn´ ees
relatives `
a une r´ egion particuli` ere.
> engrais.region1 <- summary(aov(rdt∼engrais,subset=
+
region=="Region 1"))
> engrais.region1
Df Sum Sq Mean Sq F value Pr(>F)
engrais
2
14
7.000
3 0.125
Residuals
6
14
2.333
Attention
Le test effectu´ e dans ce tableau d’analyse de la variance correspond
` a une ANOVA `
a 1 facteur (engrais) sur le rendement d’´ epeautre de
la r´ egion 1. Cependant, il ne prend pas en compte l’information des
donn´ ees pour les autres r´ egions qui permettent une meilleure estimation de la variance r´ esiduelle. Pour tester l’effet de l’engrais dans la
Le logiciel R - Effectuer des analyses statistiques
Anova Table (Type II tests)
Response: rdt
Sum Sq Df F value
Pr(>F)
region
327.19 3 112.1810 0.00000000000002955 ***
engrais
0.89 2
0.4571
0.6385
region:engrais 99.56 6 17.0667 0.00000013594213061 ***
Residuals
23.33 24
--Signif. codes: 0 ‘ *** ’ 0.001 ‘ ** ’ 0.01 ‘ * ’ 0.05 ‘.’ 0.1 ‘ ’ 1
Remarque
La formule region*engrais, utilis´ ee dans aov() et lm(), correspond
en fait `
a la formule region+engrais+region:engrais, c’est-` a-dire
le facteur r´ egion, le facteur engrais et l’interaction entre ces deux
facteurs.
La valeur-p associ´ ee au test de la pr´ esence d’interaction est significative.
Cela entraˆ ıne par exemple que l’effet de l’engrais sur le rendement peut
ˆ etre diff´ erent selon les r´ egions.
Attention
Nous consid´ erons qu’il n’y a aucun effet d’interaction si la valeur-p
associ´ ee est sup´ erieure ` a 5 %. Dans ce cas, nous effectuons une analyse
de la variance sans terme d’interaction, ce qui rend l’effet principal
plus facile ` a interpr´ eter. En pr´ esence d’interaction, il ne faut pas interpr´ eter les tests des effets principaux qui sont pr´ esent´ es dans les
sorties de la table d’analyse de la variance.
Test des effets conditionnels en pr´ esence d’interaction : On souhaite par
exemple savoir s’il existe un effet engrais dans la r´ egion 1. Pour cela, on
peut utiliser la fonction subset() qui permet de n’utiliser que les donn´ ees
relatives `
a une r´ egion particuli` ere.
> engrais.region1 <- summary(aov(rdt∼engrais,subset=
+
region=="Region 1"))
> engrais.region1
Df Sum Sq Mean Sq F value Pr(>F)
engrais
2
14
7.000
3 0.125
Residuals
6
14
2.333
Attention
Le test effectu´ e dans ce tableau d’analyse de la variance correspond
` a une ANOVA `
a 1 facteur (engrais) sur le rendement d’´ epeautre de
la r´ egion 1. Cependant, il ne prend pas en compte l’information des
donn´ ees pour les autres r´ egions qui permettent une meilleure estimation de la variance r´ esiduelle. Pour tester l’effet de l’engrais dans la