从 Ordistep 重现结果

Reproducing results from Ordistep

R 相对较新并且是第一次发帖,所以如果我的问题中缺少某些内容,我深表歉意。 我使用 Vegan 的 ordistep 函数使用默认的 "both" 方向方法进行变量选择。 我使用 Hellinger 转换的物种丰度数据框作为响应变量和 12 列自变量数据框。 如果我多次 运行 相同的代码(如下),我似乎会得到不同的包含变量。我假设有一个 set.seed 无法约束的伪数生成器。有没有办法重现一致的结果?

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

谢谢。

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

ordistep 的当前版本没有 pstep 参数。虽然我们现在有一个新的实现,但我们一直尊重 R set.seed().