Я начал использовать apaTables для работы.И мне было кое-что интересно.
APAtable вычисляет другие значения SS, MS и p, чем функция psych :: aov в R. Модель представляет собой анову 2x2.ген кодируется 0,1, пол 1,2.
- APAtable Output
Результаты ANOVA с использованием PSST0score в качестве зависимой переменной
Predictor SS df MS F p partial_eta2 CI_90_partial_eta2
(Intercept) 761.88 1 761.88 67.13 .000
gene 97.64 1 97.64 8.60 .005 .11 [.02, .23]
sex 97.35 1 97.35 8.58 .005 .11 [.02, .23]
gene x sex 77.46 1 77.46 6.83 .011 .09 [.01, .21]
Error 783.09 69 11.35
Вывод R формулы
Df Sum Sq Mean Sq F value Pr(>F)
gene 1 23.5 23.51 2.072 0.155
sex 1 20.3 20.26 1.786 0.186
gene:sex 1 77.5 77.46 6.825 0.011 *
Residuals 69 783.1 11.35
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Когда я сравнил вывод R с выходом Apatables, я получил разные значения p и сумму значений квадратов.Мне интересно, почему это так.какие-либо предложения?на какой вывод мне положиться?спасибо за ваш ответ.
ОБНОВЛЕНИЕ: после перезапуска модели с помощью функции lm я получаю те же p-значения, что и пакет ApaTable.ОБНОВЛЕНИЕ 2: я думаю, что apaTable не обрабатывает биномиальные переменные как факторы автоматически, тогда как lm / aov делает.Кто-нибудь может это подтвердить?