0

我是 R 新手。我编写了这个非常简单的脚本来突出我的问题。如果我运行这个常规 for 循环测试数据会在每次迭代时按照我想要的方式更新。

a = 5
b = 4
c = 3
testdata = matrix(nrow=100, ncol=5)
for(j in 1:100){
testdata[j,1] <- a*j
testdata[j,2] <- b*j
testdata[j,3] <- c*j
testdata[j,4] <- (a+b)*j
testdata[j,5] <- (a+c)*j
}

然而，这个使用 foreach 的并行版本完成了计算，但它们没有在 testdata 中更新。 

a = 5
b = 4
c = 3
testdata = matrix(nrow=100, ncol=5)
library(foreach)
library(doParallel)
library(doMC)
registerDoMC()
getDoParWorkers()  # Checking the number of cores. 

foreach(j = 1:100) %dopar% {
  testdata[j,1] <- a*j
  testdata[j,2] <- b*j
  testdata[j,3] <- c*j
  testdata[j,4] <- (a+b)*j
  testdata[j,5] <- (a+c)*j
}

我试图在互联网上和其他地方关注示例，但大多数示例在 R shoptalk 中太深了，我无法理解。我怎样才能让这个并行版本做非并行版本做的事情。谢谢。

4

1 回答 1

0

您应该查看该foreach软件包的文档。在foreach(j = 100)您的代码部分，您可以指定参数.combine来告诉foreach如何编译您的结果。由于您想要一个 5x100 的数据框/矩阵，因此您可以在逻辑上编写一个包含五个参数（即c(a*j, b*j, c*j, (a+b)*j, (a+c)*j)）的向量，然后rbind将它们组成一个数据框。在下面查看我的代码：

a = 5
b = 4
c = 3

library(foreach)
library(doParallel) 
library(parallel)

## Assuming you want to use all of your cores
registerDoParallel(cores = detectCores())

## Specify your .combine argument below
foreach(j = 1:100, .combine = "rbind") %dopar% {
    c(a*j, b*j, c*j, (a+b)*j, (a+c)*j)
}

这吐出：

           [,1] [,2] [,3] [,4] [,5]
result.1      5    4    3    9    8
result.2     10    8    6   18   16
result.3     15   12    9   27   24
result.4     20   16   12   36   32
result.5     25   20   15   45   40
...

然后，您可以通过将其分配给您想要的变量来更进一步：

...
testdata <- foreach(j = 1:100, .combine = "rbind") %dopar% {
                c(a*j, b*j, c*j, (a+b)*j, (a+c)*j)
            }
testdata <- as.data.frame(testdata, row.names = FALSE)

希望这可以帮助！

于 2016-03-23T23:46:12.367 回答