多重插补数据包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