R中的并行随机数发生器

Question

R中是否有Mersenne-Twister和Wichmann-Hill并行随机数发生器的实现？ 我使用parallel包，但似乎只有L＆＃39; Ecuyer并行随机数生成器的实现。

Answer 1

R中的默认随机数生成器是Mersenne-Twister。您可以使用

在它们之间进行更改

setRNG('Wichmann-Hill')
setRNG('default')#or setRNG('Mersenne-Twister')

如果您想并行生成数字，可以使用foreach包。

require(foreach)
require(doParallel)
c1 <- makeCluster(2)
registerDoParallel(c1)
generateRandom <- function(rng='default',n) {
  setRNG(rng)
  runif(n)
}
result = foreach(i = 1:2,rng = c('default','Wichmann-Hill'), .combine = 'c') %dopar%
  generateRandom(rng,1e7)
stopCluster(c1)

此示例从U（0,1）分布生成2000万个随机数，每一半使用不同的算法。