Я выполнил частичную 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 объясняющих переменных, было бы довольно много времени вычислять каждую комбинацию параметров.