我正在使用 vegan 包中的 adonis 函数来确定不同环境因素对不同地区森林植物群落组成的影响。我想首先使用 adonis 来消除区域效应，这是为了适合像这样的模型
adonis_region<- adonis(community ~ region,permutations=999,method="bray")
其中社区是属于不同区域的森林斑块中物种的存在/不存在数据矩阵。
然后我想使用这个模型的残差作为其他一些分析的响应变量。我知道我可以使用地层将区域建模为块变量，但这对我不感兴趣，因为之后我想执行另一种分析，我希望消除区域效应。
我的问题是我无法弄清楚如何从 adonis 模型中获取残差值。
任何可以解决此问题的其他类型的多变量分析的提示都非常受欢迎。