У меня смешанный дизайн с 1 между фактором (Группа: ASD, CTR) и тремя внутри факторов (Время: до, после; Доверительный управляющий: хорошо, плохо; Шаг: 1,2,3,4,5).
Я выполнил смешанную модель anova (тип III) с aov_car:
anova <- aov_car(rating ~ time*trustee*step*group + Error(id/time*trustee*step), data = datatemp, anova_table = list(es = "ges"))
Я получил значительное трехстороннее взаимодействие: Время * Попечитель * Шаг
Я хочу разложитьвзаимодействие, делая запланированные контрасты для времени (меня просто интересуют контрасты между до и после для разных опекунов и этапов).
Это код, который я использовал:
trustee_time_step <- emmeans(anova, ~ time | trustee | step)
pairs(trustee_time_step)
результаты для сравнения со временем
По умолчанию, pair () должен выполнить коррекцию Tuckey.Однако коррекция контрастов не производится.
Коррекция Тукки появляется только в том случае, если я использую шаг для разделения контрастов:
step_time_trustee <- emmeans(anova, ~ step | trustee | time)
pairs(step_time_trustee)
результаты для сравнения с шагом
Я полагаю, это потому, что у Времени есть только два уровня (каждый контраст просто пост-против-до), в то время как у шага 5 уровней (см. Метод Тьюки для сравнения семейства из 5 человек).оценки »в конце рисунка).
Так вот мой вопрос: почему pair () не применяет коррекцию в первом случае?Даже если пост против предварительно - только один контраст, это сделано 10 раз (один для каждой комбинации опекуна и шага), и, таким образом, это должно быть исправлено по моему мнению.Что я должен делать?Должен ли я исправить их вручную?И если да, то как?
Я очень озадачен этим вопросом и не хочу сообщать о нескорректированных результатах, если это не является "статистически хорошим".
Надеюсь, мой вопросдостаточно ясно.Спасибо за помощь!