phyloseq пошаговый выбор вперед
Я работаю с многомерными данными, пытаясь связать разнообразие бактерий в почве (полученное путем метабаркодирования) с факторами окружающей среды, чтобы увидеть, какие факторы объясняют наибольшее разнообразие бактериального разнообразия. Я использую пакет phyloseq для большей части анализа. Мне удалось закодировать различные типы рукоположения, например,
ps.rda <- ordinate(ps, formula = ps ~ ., method = "RDA")
Однако я хотел бы сделать RDA с поэтапным выбором значимых переменных среды. (словно ordistep() функция в vegan пакет). Вы знаете, как кодировать это в пакете phyloseq? Любые советы высоко ценится.