R中的并行随机数发生器

Question

R中有Mersenne-Twister和Wichmann-Hill并行随机数发生器的实现吗？我正在使用parallel包，但似乎只有L‘厄瓜多尔并行随机数生成器的实现。

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万个随机数，每一半使用不同的算法。