如何将R语言for循环正确转换为foreach循环并保持输出结构一致？

原代码逻辑 #

原for循环的代码如下：

x1 <- array(0, c(10, 20))
for(indv in 1:10)
{
  x1[indv,] <- cbind(c(1:10) * indv * 1.0, c(31:40))
}
array(x1)

这里的核心细节是：每次循环生成的是10行2列的矩阵，但赋值给x1的某一行（1行20列）时，R会自动将这个矩阵按列展平为1行20列的向量。最终x1是10行20列的矩阵，调用array(x1)时，R会按默认的列优先顺序把矩阵拉成一维数组——这就是原输出后半部分会出现31重复10次、32重复10次……直到40重复10次的原因。

原foreach代码的问题 #

你尝试的foreach代码：

x2 <- foreach(indv = 1:10, .combine=rbind) %do% {
  cbind(c(1:10) * indv * 1.0, c(31:40))
}
array(x2)

问题出在：每次循环返回的还是10行2列的矩阵，用rbind合并后得到的是100行2列的矩阵。再调用array(x2)时，按列优先展平，后半部分就变成了31-40循环10次，和原逻辑完全不符。

修正后的foreach代码 #

要匹配原逻辑，需要先把每次循环生成的矩阵展平成1行20列的向量，再按行合并。修改后的代码如下：

library(foreach)
x2 <- foreach(indv = 1:10, .combine=rbind) %do% {
  # 显式将10×2矩阵按列展平为1×20的向量
  as.vector(cbind(c(1:10) * indv * 1.0, c(31:40)))
}
array(x2)

运行这段代码后，array(x2)的输出结构与原array(x1)完全一致。

另一种等效写法：先按列合并所有向量，再转置矩阵，最后展平：

x2 <- foreach(indv = 1:10, .combine=cbind) %do% {
  as.vector(cbind(c(1:10) * indv * 1.0, c(31:40)))
}
array(t(x2))

内容的提问来源于stack exchange，提问作者Deamonpog