我正在处理多变量数据，试图将土壤细菌多样性（通过元条形码获得）与环境因素联系起来，看看哪些因素可以解释细菌多样性的最大变化。我正在使用 phyloseq 包进行大部分分析。我设法通过例如编写不同类型的协调
ps.rda <- ordinate(ps, formula = ps ~ ., method = "RDA")
但是，我想做一个 RDA，逐步选择重要的环境变量。（如包ordistep()中的功能vegan）。您知道在 phyloseq 包中对此进行编码的方法吗？任何意见是极大的赞赏。