1

我正在使用包中的adonis2函数运行 perm-manova vegan。我的模型中有一个重要的交互,所以想对使用lsmeans包的交互进行 tukey 调整比较。

但是,虽然lsmeans接受许多模型类,但它不接受adonis对象的类 ( [1] "anova.cca" "anova" "data.frame")。

有没有办法可以将我的adonis对象强制为lm, aov,manova对象(或任何其他接受的类lsmeans),以便我可以使用这个函数?谢谢

library(vegan)
library(lsmeans)

data(dune)
data(dune.env)

man<-adonis2(dune ~ Management*A1, data = dune.env)

lsmeans(man,pairwise ~ Management:A1,adjust="tukey")

Error in ref.grid(object = list(Df = c(3, 1, 3, 12), SumOfSqs = c(1.46859175179317,  : 
  Can't handle an object of class  “anova.cca” 
 Use help("models", package = "lsmeans") for information on supported models.
4

0 回答 0