我在vegan
软件包的帮助下执行了部分RDA。我知道我可以使用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个解释变量,因此计算每个参数组合都将非常耗时。