Выбор модели для частичной RDA
Я выполнил частичный RDA с помощью vegan пакет. Я знаю, что могу выбрать самую экономную модель с ordiR2step функция. OrdiR2step использует скорректированный R2 как критерий соответствия. Однако, очевидно, эта функция не работает для частичных RDA. Об этом говорится в документальном фильме о веганах, также я получаю сообщение об ошибке при попытке следующего кода:
> step_both_R2 <- ordiR2step(rda(species_data ~ treatment +
> Condition(stream), data=chemical_data), scope = formula(rda_1),
> direction="both", pstep = 1000 )
Кто-нибудь знает какой-либо другой способ найти наиболее экономную модель для частичного RDA?
Поскольку в моей модели около 10 объясняющих переменных, было бы довольно много времени вычислять каждую комбинацию параметров.