Функция aov() в R - ANOVA
Примеры на языке RФункция aov() в R используется для выполнения дисперсионного анализа (ANOVA), чтобы сравнить, есть ли значительные различия между средними значениями трех или более групп.
Когда нужно сравнить более двух групп данных, использование ANOVA более уместно, чем множественные t-тесты.
Синтаксис функции aov() выглядит следующим образом:
aov(formula, data)
Описание параметра:
formulaФормула в формате переменная отклика ~ группирующая переменная.
dataФрейм данных, содержащий переменные из формулы.
Пример
# Влияние трех удобрений на урожайность сельскохозяйственных культур
fertilizer <- factor(rep(c("A", "B", "C"), each = 5))
yield <- c(18, 20, 19, 22, 21, # Удобрение A
25, 27, 24, 26, 28, # Удобрение B
15, 17, 16, 14, 18) # Удобрение C
# Создать датафрейм
df <- data.frame(fertilizer, yield)
print("Предварительный просмотр данных:")
print(df)
# Выполнить 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, # Удобрение A
25, 27, 24, 26, 28, # Удобрение B
15, 17, 16, 14, 18) # Удобрение C
# Создать датафрейм
df <- data.frame(fertilizer, yield)
print("Предварительный просмотр данных:")
print(df)
# Выполнить ANOVA
result <- aov(yield ~ fertilizer, data = df)
print(summary(result))
Выполнение приведённого выше кода даёт следующий вывод:
[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
p-значение намного меньше 0.001, что указывает на то, что три удобрения имеют высокозначимые различия по влиянию на урожайность.
Другие расширения