R的多个实例与MICE产生相同的推断——随机种子问题



多重插补数据包mice使用随机数生成器,该生成器可以使用set.seed参数进行设置。我正在运行R(即3)的多个实例,其中用相同的mice参数估算相同的数据,其中set.seed=NA不指定任何种子值。

然而,R的三个例子产生了完全相同的推论。我对这一观察结果感到困惑,因为我认为在使用set.seed=NA时,R的不同实例应该使用不同的种子值,但它们似乎使用相同的种子值。当我手动设置set.seed时,会产生不同的输入。

我的问题是:这是否合理,或者可能存在计算机(PC windows 7)问题,导致种子相同,即使我设置了set.seed=NA。我欢迎对这个问题有任何想法。

顺便说一句:我使用了R的多个实例,因为MICE函数ibind()可以用于组合mids(乘法估算数据集)对象。因此,它节省了多核系统的计算时间。如果产生同样的指责,这当然是无稽之谈。

如果使用NULL作为set.seed的参数,则RNG将产生不同的结果。

> set.seed(NULL)
> rnorm(5)
[1] -0.45861478 -0.56525262 -0.95349202  0.05073898 -1.01956671
> set.seed(NULL)
> rnorm(5)
[1]  0.57787302  3.27165390  0.53466909 -1.15390604 -0.02624054

我建议使用序列

set.seed(NULL)
mice(...) # without seed argument

最新更新