我已经开始使用apaTables工作。我一直在想一些事情。
与R中的psych :: aov函数相比,APAtable计算的SS,MS和p值不同。
该模型是2x2方差分析。基因编码为0,1,性别为1,2。
使用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函数重新运行模型后,我得到了与ApaTable包相同的p值。 UPDATE2:我猜apaTable不会自动将二项式变量视为因子,而lms / aovs会这样做。有人可以确认吗?