从Foreach获取输出等于R中for循环的结果

Question

我目前正在将for 循环转换为foreach 循环，以便我可以并行化内部循环的计算。我能够做到这一点并获得正确的输出。但是，我希望foreach 的输出格式与我在旧for 循环中的格式相同。这是我想要做的一个可重复的例子：

FOR 循环：

x <- c(1:10)
y <- c(5)
z <- c(3)

result_list <- list()
for(i in length(x)){
        for(j in 1:y){
                for(k in 1:z){
                        
                result <- rnorm(n = 250)
                result_list[[paste(i, j, k)]] <- result
                }
        }
}

这就是我得到的输出和我想要实现的目标：

library(dplyr)
glimpse(result_list)
List of 15
 $ 10 1 1: num [1:250] -0.382 -0.156 -0.747 1.139 -1.824 ...
 $ 10 1 2: num [1:250] -0.181 0.526 0.255 -0.369 -1.517 ...
 $ 10 1 3: num [1:250] 0.3621 1.3634 -0.0507 -0.3943 -1.2183 ...
 $ 10 2 1: num [1:250] -1.453 0.731 2.761 -0.586 -1.631 ...
 $ 10 2 2: num [1:250] -0.663 -0.641 -1.852 1.58 0.133 ...
 $ 10 2 3: num [1:250] -1.334 0.803 -0.116 0.618 -1.339 ...
 $ 10 3 1: num [1:250] 0.158 -1.296 -0.947 -0.515 -0.208 ...
 $ 10 3 2: num [1:250] -0.604 1.956 0.127 1.846 -0.549 ...
 $ 10 3 3: num [1:250] 0.365 -0.467 -0.589 -1.251 0.491 ...
 $ 10 4 1: num [1:250] -1.138 0.883 1 0.729 -1.566 ...
 $ 10 4 2: num [1:250] -0.0461 2.3096 -1.5347 0.3722 0.3125 ...
 $ 10 4 3: num [1:250] 0.127 -0.728 0.402 1.783 -1.457 ...
 $ 10 5 1: num [1:250] 1.855 2.224 1.301 0.166 -0.28 ...
 $ 10 5 2: num [1:250] 0.463 -1.011 1.067 -1.305 -0.51 ...
 $ 10 5 3: num [1:250] 1.937 0.651 -0.424 -0.714 -0.225 ...

这是 FOREACH 循环的代码：

library(foreach)
result_list <- list()
resultado_foreach <- foreach::foreach(i = 1:10, .inorder = TRUE) %do% {
        
                foreach::foreach(j = 1:5, .inorder = TRUE) %do% {
                        foreach::foreach((k = 1:3)) %dopar% {
                                
                                result <- rnorm(n = 250)
                                result_list[[paste(i, j, k)]] <- result
                                
                        }
                
        }
}

虽然，我得到了一个列表列表（嵌套列表）：

glimpse(resultado_foreach)
List of 10
 $ :List of 5
  ..$ :List of 3
  .. ..$ : num [1:250] 0.911 0.594 -0.453 0.651 2.303 ...
  .. ..$ : num [1:250] -0.664 -0.696 0.741 -2.78 -0.992 ...
  .. ..$ : num [1:250] -0.6877 -1.1266 -1.7784 0.473 0.0185 ...
  ..$ :List of 3
  .. ..$ : num [1:250] 1.273 0.129 0.902 2.47 0.177 ...
  .. ..$ : num [1:250] 0.705 0.519 1.219 -1.682 -0.355 ...
  .. ..$ : num [1:250] 1.138 0.422 -1.025 -0.237 0.418 ...
  ..$ :List of 3
  .. ..$ : num [1:250] -1.636 -1.297 -1.115 -0.138 0.174 ...
  .. ..$ : num [1:250] 0.56 -1.311 0.641 0.861 -0.601 ...
  .. ..$ : num [1:250] 0.198 -1.197 0.781 -0.571 -0.141 ...
  ..$ :List of 3
  .. ..$ : num [1:250] -0.355 -0.649 -1.046 -0.717 -0.97 ...
  .. ..$ : num [1:250] -1.086 0.912 -0.996 0.303 1.418 ...
  .. ..$ : num [1:250] 0.8827 -0.0761 1.3504 -0.5301 0.2267 ...
  ..$ :List of 3
  .. ..$ : num [1:250] -1.826 1.286 1.585 -0.359 -0.955 ...
  .. ..$ : num [1:250] -0.588 -0.314 -0.223 -0.779 0.569 ...
  .. ..$ : num [1:250] 1.047 -0.242 -0.345 0.27 -0.158 ...

foreach 输出比我放在这里的要长得多。

我已经尝试了很多组合，并在foreach 函数中为.combine 参数设置了很多函数。那么，如何获得与for loop相同格式的输出呢？

Answer 1

目前在你的 for 循环中

for(i in length(x)) {...}

但是length(x) 是 10，所以你只是在做for(i in 10)，它只循环一次。这与for (i in 1:10) 不同。更安全的选择是for (i in seq_along(x))。

假设您期望有 10 个外循环，每个循环有 5 个中间循​​环，每个循环有 3 个内循环。这应该是 10x5x3 = 150 个总结果。

如果你想使用 for each 的嵌套循环，你应该使用%:% 操作符。另外，您还需要使用.combine=c 将内部循环连接到最后的单个列表中。这些选项在?foreach 帮助页面上进行了讨论。

更好的版本是

resultado_foreach <- 
  foreach::foreach(i = 1:10, .inorder = TRUE, .combine=c) %:%
    foreach::foreach(j = 1:5, .inorder = TRUE, .combine=c) %:% 
      foreach::foreach((k = 1:3)) %dopar% {
        rnorm(n = 250)
      }

这将返回一个长度为 150 的列表，每个列表带有一个长度为 250 的向量。 请注意，使用foreach 时不应修改全局变量。每个块都应该返回一个值，foreach 将为您收集和组合这些值。