R中的并行随机数生成器

Question

R 中是否有 Mersenne-Twister 和 Wichmann-Hill 并行随机数生成器的实现？ 我正在使用 parallel 包，但似乎只有 L'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 万个随机数，每一半使用不同的算法。