r - 从 Ordistep 复制结果

Question

R 相对较新，并且是第一次发布，所以如果我的问题中缺少某些内容，我深表歉意。我正在使用 Vegan 的ordistep函数使用默认的“双向”方向方法进行变量选择。我使用 Hellinger 转换的 Species Abundance 数据框作为响应变量和 12 列自变量数据框。如果我多次运行相同的代码（如下），我似乎会得到不同的包含变量。我假设涉及 set.seed 无法约束的伪数生成器。有没有办法重现一致的结果？

set.seed(2000)
step.both <- ordistep(pre_met.rda, pstep=1000)

谢谢你。

Answer 1

ordistep使用标准 R RNG 和 Honors set.seed()。我无法重现您的问题，但是当我在调用之前将种子设置为相同的值时，我得到了完全相同的结果。

当前版本ordistep没有pstep参数。尽管我们现在有了新的实现，但我们一直尊重 R set.seed()。