r - 在 R 和 Julia 中生成相同的随机数

Question

我想在 R 和 Julia 中生成相同的随机数。默认情况下，这两种语言似乎都使用 Mersenne-Twister 库，但在 Julia 1.0.0 中:

julia> using Random
julia> Random.seed!(3)
julia> rand()
0.8116984049958615

生产 0.811... , 而在 R 中:

set.seed(3)
runif(1)

生产 0.168 .

有任何想法吗？

相关 SO 问题 here和 here .

我对那些感兴趣的人的用例:通过将输出与 R 中等效库的输出进行比较来测试需要随机数生成(例如统计引导)的新 Julia 代码。

Answer 1

那是个老问题。

Paul Gilbert 在 1990 年代后期 (!!) 试图断言 R 中的模拟(当时是新人)给出的结果与 S-Plus(当时的现任者)中的结果相同时，他解决了同样的问题。

他的解决方案，仍然是 AFAICT 的黄金方法:在两种语言中重新实现新代码，因为这是确保相同的种子、状态......以及其他任何影响它的唯一方法。