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