R包ap​​aTables方差分析输出与psych :: aov不同

时间:2018-11-29 11:50:21

标签: r export anova

我已经开始使用apaTables工作。我一直在想一些事情。
与R中的psych :: aov函数相比,APAtable计算的SS,MS和p值不同。 该模型是2x2方差分析。基因编码为0,1,性别为1,2。

  1. 可分配的输出

使用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
  1. 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会这样做。有人可以确认吗?

0 个答案:

没有答案