最簡單的做法並行重複

Question

我喜歡的parallel包R和是多麼容易和直觀做並行版本apply，sapply等

是否有replicate類似的並行功能？

Answer 1

你可以只使用lapply或sapply並行版本，而不是說要複製這個表達式n時候，你做1:n，取而代之的是應用給出一個表達式，你將該表達式包裝在忽略發送給它的參數的函數中。

可能是這樣的：

#create cluster 
library(parallel) 
cl <- makeCluster(detectCores()-1) 
#get library support needed to run the code 
clusterEvalQ(cl,library(repsych)) 
#put objects in place that might be needed for the code 
clusterExport(cl,c("myData")) 
#... then parallel replicate... 
parSapply(cl, 1:10000, function(i,...) { x <- rnorm(10); mean(x)/sd(x) }) 
#stop the cluster 
stopCluster(cl) 

爲並聯等效：

replicate(10000, {x <- rnorm(10); mean(x)/sd(x) }) 

Answer 2

這是最好的我能想出：

cl <- makeCluster(getOption("cl.cores", 4)) 
clusterCall(cl, replicate(50, simulate_fxns())) 
stopCluster(cl) 

Answer 3

使用clusterEvalQ作爲一種模式，我想我會實現一個並行replicate爲：

parReplicate <- function(cl, n, expr, simplify=TRUE, USE.NAMES=TRUE) 
    parSapply(cl, integer(n), function(i, ex) eval(ex, envir=.GlobalEnv), 
      substitute(expr), simplify=simplify, USE.NAMES=USE.NAMES) 

參數simplify和USE.NAMES與sapply兼容，而不是replicate，但在我看來，它們使其成爲parSapply的更好包裝。

下面是來自replicate手冊頁導出一個例子：

library(parallel) 
cl <- makePSOCKcluster(3) 
hist(parReplicate(cl, 100, mean(rexp(10))))