我正在使用包中的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.