Ordistep的复制结果

Question

相对较新的R和第一次发布，所以如果我的问题中缺少某些东西，我很抱歉。 我使用纯素的ordistep函数进行变量选择，使用默认的“双向”方向方法。 我使用Hellinger转换的物种丰度数据框作为响应变量和一个12列独立变量数据框。 我似乎得到不同的包含变量，如果我多次运行相同的代码（如下）。我假设有一个伪代码生成器涉及set.seed不能约束。有没有办法重现一致的结果？

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

谢谢。

Answer 1

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

当前版本的ordistep没有pstep参数。虽然我们现在有了一个新的实施方案，但我们一直很荣幸R set.seed()。