如何用固定的嵌套因子解释permanova adonis的输出?

时间:2019-01-17 21:47:51

标签: vegan

当同时使用嵌套因子和交互时,我不太确定如何在adonis2输出(纯素食包)中解释“:”。

例如,没有嵌套因素:

        adonis_data <- adonis2(Y ~ 
                       A*B*C , data=factors, 
                       permutations = 9999, method = "euclidean")

        adonis_data

因子之间的相互作用在结果输出中用“:”表示,例如:

   Terms added sequentially (first to last)
        A       ***
        B       ***
        C       ***
        A:C     ***
        B:C     ***

但是,在集成嵌套效果(使用“:”或“ /”)时:

        adonis_data <- adonis2(Y ~ 
                       A:B*C , data=factors, 
                       permutations = 9999, method = "euclidean")

        adonis_data

我难以解释此输出,因为我不知道这里显示的嵌套因子(带或不带“:”)

   Terms added sequentially (first to last)
        A      ***
        C      ***
        A:B    ***
        A:C    ***
        A:C:B  ***

所以在这里,我不太确定A:B是否显示了嵌套因子的作用而不是交互作用...

我对此事的最后想法是,如果嵌套因子被写为输出中的相互作用(例如A:B),则考虑到这些术语是按顺序添加的,我们可以通过结果的顺序来区分这两种情况。

所以,如果这是正确的,那意味着嵌套的B应该在A和C之间,而A:B在这里是一种交互作用?

在这里得出A对Y的B效应根据A值而不同是否正确?

有人能启发我吗? 非常感谢, 劳拉

1 个答案:

答案 0 :(得分:0)

这是 R 公式的常见问题,对adonis纯素食主义者而言,并不是什么特别的问题。无论如何,在默认的adonis2(目前仅为adonis)中,测试是顺序的,这实际上意味着之前的术语会影响结果,但是以后(随后的)条款无效。术语是什么,取决于您定义模型的方式,就像您在问题中所展示的那样。