TukeyHSD после внутривенных факторов ANOVA

Question

Я запускаю следующие команды

a1 <- aov(Ratio~(Condition*Event) + Error(Participant/(Condition*Event)), Data)
summary(a1)

, которые выводят правильные результаты (таблицы ANOVA и т. Д.).

При запуске:

TukeyHSD (a1, "Condition")

Я получаю: Ошибка в UseMethod ("TukeyHSD"): нет применимого метода для "TukeyHSD"

Почему работает ANOVA, но не TukeyHSD?Обе переменные фактора являются строками (Условие имеет 3 уровня, а Событие имеет 4 уровня).

РЕДАКТИРОВАТЬ: Когда я повторяю aov без термина Error, это работает, однако Тьюки не показывает существенной разницы между любой из пар (ANOVA был значимым для уверенности).Значит ли это, что Тьюки корректирует множественные сравнения?

Answer 1

aov() с Error() возвращает объект с классом "aovlist" "listof". У TukeyHSD, по-видимому, нет метода ни для одного из этих классов.

Answer 2

В качестве альтернативы традиционному повторному измерению ANOVA для внутрисубъектного дизайна, вы можете рассмотреть возможность использования линейного подхода со смешанными эффектами. Он все чаще используется в научном сообществе и позволяет избежать некоторых ловушек ANOVA, допуская при этом более сложные структуры ошибок. С переменной переменной ответа достаточно пакета nlme , но вы также можете использовать lme4 , что дополнительно позволяет справиться с категориальными переменными ответа. Для множественных сравнений (включая специальные тесты Тьюки) пакет multcomp (см. Функцию glht()) можно использовать с моделями со смешанными эффектами, оснащенными nlme::lme, как описано здесь: Повторные измерения ANOVA с использованием R .

Одно краткое замечание о вашем дизайне: если ваша (зависимая) переменная ответа, Ratio, является пропорцией или ограниченным значением, вы можете подумать об использовании другой функции связи.

Answer 3

library(laercio)
some.aov <- aov(Ratio~(Condition*Event) + Error(Participant/(Condition*Event)), Data)
anova(some.aov)
LTukey(some.aov,"Condition")