Апостериорный t-тест в R
Я сделал между субъектами повторное измерение ANOVA в R, дав мне те же результаты, что и SPSS.
Дизайн выглядит следующим образом: Позитивность (2 уровня) x Возникновение (2 уровня) x Группа (3 уровня)
В результате у меня трехстороннее взаимодействие. Теперь я хочу сделать все возможные дополнительные тесты, отдельно для трех групп. Это возможно с помощью следующего кода и использования пакетов lsmeans и afex:
anova_all <- aov_ez("code", "err", stim.lng, +
within=.("recency", "positivity"), between=.("group"))
ref <- lsmeans(anova_all,~recency*positivity|group)
comps <- contrast(ref,method="pairwise", adjust="none")
summary(comps)
Код выполняется, но я получаю другие результаты по сравнению с выводом SPSS (с точки зрения p-значений и t-значений). Это смущает меня, и я считаю, что я должен получить результаты, аналогичные SPSS, поскольку проводятся только t-тесты парных образцов.
Я предполагаю, что решающим моментом могут быть разные степени свободы в обоих программах. Например, когда я выполняю парный t-тест в конкретной группе из 9 субъектов, в SPSS степень свободы составляет 9, тогда как она находится между 89-102 в R (см. Прикрепленные изображения). Почему это так высоко в R? Есть ли еще какие-либо репликации результатов SPSS с R?
Вот мой дизайн и три различных результата из SPSS и R в конкретной группе пациентов для всех возможных парных специальных тестов:
Выход из R:
Выход из SPSS:
