Дисперсионный анализ (ANOVA)
До сих пор вы создавали объекты с результатами вычислений, чтобы сохранять и извлекать из них элементы для удобного и детального представления результатов. Однако некоторые такие объекты — например, возвращаемые функцией aov() при проведении дисперсионного анализа (ANOVA) — можно использовать как входные данные для других функций: summary() и plot() для форматированного вывода и визуализации результатов. Кроме того, объект aov() можно передать в функцию TukeyHSD() для выполнения пост-хок попарного тестирования.
Набор данных abaloneKeep и пакет dplyr уже загружены.
Это упражнение является частью курса
R для пользователей SAS
Инструкции к упражнению
- Вычислите количество наблюдений
n, среднее значение, стандартное отклонение и дисперсию возраста (age) по группам пола (sex). - Выполните дисперсионный анализ переменной
ageпо группамsexи сохраните результат в переменнуюabaov. - Запустите функцию
summary()для объектаabaov. - Выполните попарные пост-хок тесты с помощью
TukeyHSD()для объектаabaov.
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# Compute n(), mean(), sd() and var() of age by sex
abaloneKeep %>% group_by(___) %>% select(___, ___) %>%
summarise(across(everything(), list(mean = ___,
___ = ~sd(.x),
var = ___)),
N = ___)
# Run aov() of age by sex, save as abaov
abaov <- aov(___ ~ ___, data = ___)
# Run summary() of abaov
___
# Perform TukeyHSD posthoc pairwise tests on abaov
TukeyHSD(___)