Fonction aov() de R - ANOVA
Exemples du langage RLa fonction aov() de R est utilisée pour effectuer une analyse de variance (ANOVA) afin de comparer s'il existe des différences significatives entre les moyennes de trois groupes ou plus.
Lorsqu'il y a plus de deux groupes de données à comparer, utiliser une ANOVA est plus approprié que des tests t multiples.
La syntaxe de la fonction aov() est la suivante :
aov(formula, data)
Description du paramètre :
formulaFormule, au format variable réponse ~ variable de regroupement.
dataData frame, contenant les variables de la formule.
Exemple
# Effet de trois engrais sur le rendement des cultures
fertilizer <- factor(rep(c("A", "B", "C"), each = 5))
yield <- c(18, 20, 19, 22, 21, # Engrais A
25, 27, 24, 26, 28, # Engrais B
15, 17, 16, 14, 18) # Engrais C
# Créer un DataFrame
df <- data.frame(fertilizer, yield)
print("Aperçu des données :")
print(df)
# Effectuer l'ANOVA
result <- aov(yield ~ fertilizer, data = df)
print(summary(result))
fertilizer <- factor(rep(c("A", "B", "C"), each = 5))
yield <- c(18, 20, 19, 22, 21, # Engrais A
25, 27, 24, 26, 28, # Engrais B
15, 17, 16, 14, 18) # Engrais C
# Créer un DataFrame
df <- data.frame(fertilizer, yield)
print("Aperçu des données :")
print(df)
# Effectuer l'ANOVA
result <- aov(yield ~ fertilizer, data = df)
print(summary(result))
L'exécution du code ci-dessus produit la sortie suivante :
[1] "数据预览:"
fertilizer yield
1 A 18
2 A 20
3 A 19
4 A 22
5 A 21
6 B 25
...
Df Sum Sq Mean Sq F value Pr(>F)
fertilizer 2 190.0 95.00 28.61 2.73e-05 ***
Residuals 12 39.9 3.32
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
La valeur p est bien inférieure à 0.001, ce qui indique que les trois engrais ont une différence hautement significative dans leurs effets sur le rendement.
Autres extensions