vegan::simper 中的多重比较选项

Question

我正在使用vegan 包中的adonis 函数来确定社区（PCB 同源物）中几个不同因素之间的差异差异。我还决定使用simper 函数来评估哪些社区成员对观察到的差异贡献最大。有没有办法在simper 函数中包含多个因素？我成功运行了如下所示的adonis 模型，但是simper 函数的相应代码不起作用。非常感谢。

# Adonis model
pcbtest3 <- adonis(pcbcong ~ FISH_CLASS+REACH+BASIN,
                   data=pcbcov,
                   method="bray",
                   permutations=999
                   )
# pcbcong=matrix of community dissimilarities
# FISH_CLASS,REACH, BASIN are factors
# Simper
simp <- with(pcbcov, simper(pcbcong,(FISH_CLASS, BASIN, REACH)))

Answer 1

simper 只能与一个分组变量一起使用。

一种解决方法是使用interaction() 从所有三个因素中创建一个新的分组变量。

对此的解释会很尴尬（取决于你有多少因子水平），我不知道这是否有意义 - 所以要小心。