私は多変量データを扱っており、土壌細菌の多様性 (メタバーコーディングによって取得) を環境要因に関連付けて、細菌の多様性の最も大きな変動を説明する要因を確認しようとしています。ほとんどの分析には phyloseq パッケージを使用しています。たとえば、さまざまな種類の序列をコーディングすることができました
ps.rda <- ordinate(ps, formula = ps ~ ., method = "RDA")
ただし、重要な環境変数を段階的に選択して RDA を実行したいと考えています。（パッケージordistep()内の機能のように）。veganphyloseq パッケージでこれをコーディングする方法を知っていますか? どんなアドバイスでも大歓迎です。