Я имею дело с несбалансированным дизайном/образцом и изначально узнал aov()
. Теперь я знаю, что для тестов ANOVA мне нужно использовать Type of Squares, который включает использование фикчирования с использованием lm()
вместо использования aov()
.
Проблема заключается в получении пост-hoc-тестов (в частности, Tukey HSD) с использованием lm()
. Все проведенные исследования говорят, что использование simint
в пакете multcomp
будет работать, но теперь, когда он обновил эту команду, похоже, не будет доступным. Похоже, он полагается также на то, чтобы пройти aov()
для вычисления.
По сути, все тесты Tukey HSD, которые я нашел для R, предполагают, что вы используете aov()
для сравнения, а не для lm()
. Чтобы получить сумму квадратов типа III, которые мне нужны для несбалансированного дизайна, я должен использовать:
mod<-lm(Snavg~StudentEthnicity*StudentGender)
Anova(mod, type="III")
Как использовать Tukey HSD-тест с модом, используя lm()
? Или, наоборот, вычислить мой ANOVA с помощью Type III и все еще иметь возможность запускать тест Tukey HSD?
Спасибо!