544
Le logiciel R - Effectuer des analyses statistiques
> plot(delai∼traitement)
●
●
Placebo
T2
T3
T4
T5
3
4
5
6
7
8
9
10
traitement
delai
Figure 13.1 – Boˆ ıtes ` a moustaches des d´ elais de cicatrisation pour chaque traitement.
13.1.3 Table d’ANOVA et estimations des param` etres
Instruction R pour la table d’ANOVA : La fonction ` a utiliser est aov().
Comme pour le mod` ele de r´ egression, l’ANOVA fonctionne avec des formules R. Il faut donc sp´ ecifier le mod` ele ` a utiliser.
> mon.aov <- aov(delai∼traitement)
> summary(mon.aov)
Df Sum Sq Mean Sq F value Pr(>F)
traitement
4 36.47
9.117
3.896 0.0136 *
Residuals
25 58.50
2.340
--Signif. codes: 0 ‘ *** ’ 0.001 ‘ ** ’ 0.01 ‘ * ’ 0.05 ‘.’ 0.1 ‘ ’ 1
Attention
Si cela n’est pas le cas, pensez ` a d´ eclarer votre variable facteur (ici la
variable traitement) comme un objet R du type factor au moyen
de la fonction factor().
> class(traitement)
[1] "factor"
Comme l’ANOVA est en fait un mod` ele lin´ eaire, notons qu’il est aussi possible d’effectuer l’analyse de la variance du mod` ele lin´ eaire sous-jacent :
Le logiciel R - Effectuer des analyses statistiques
> plot(delai∼traitement)
●
●
Placebo
T2
T3
T4
T5
3
4
5
6
7
8
9
10
traitement
delai
Figure 13.1 – Boˆ ıtes ` a moustaches des d´ elais de cicatrisation pour chaque traitement.
13.1.3 Table d’ANOVA et estimations des param` etres
Instruction R pour la table d’ANOVA : La fonction ` a utiliser est aov().
Comme pour le mod` ele de r´ egression, l’ANOVA fonctionne avec des formules R. Il faut donc sp´ ecifier le mod` ele ` a utiliser.
> mon.aov <- aov(delai∼traitement)
> summary(mon.aov)
Df Sum Sq Mean Sq F value Pr(>F)
traitement
4 36.47
9.117
3.896 0.0136 *
Residuals
25 58.50
2.340
--Signif. codes: 0 ‘ *** ’ 0.001 ‘ ** ’ 0.01 ‘ * ’ 0.05 ‘.’ 0.1 ‘ ’ 1
Attention
Si cela n’est pas le cas, pensez ` a d´ eclarer votre variable facteur (ici la
variable traitement) comme un objet R du type factor au moyen
de la fonction factor().
> class(traitement)
[1] "factor"
Comme l’ANOVA est en fait un mod` ele lin´ eaire, notons qu’il est aussi possible d’effectuer l’analyse de la variance du mod` ele lin´ eaire sous-jacent :