从 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()
.
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()
.