将三个 for 循环的结果保存为向量

Question

[假设我们有三个向量b1、b2 和b3，我需要遍历它们所有可能的值组合并将它们的乘积存储在一个新向量中。例如，

b1 <- seq(0.2, 2, by=0.2)
b2 <- seq(0.2, 2, by=0.2)
b3 <- seq(0.2, 2, by=0.2)

product <- vector()

for (i in 1:length(b1))
{
  for (j in 1:length(b2))
  {
    for (k in 1:length(b3))
    {
      product[i] <- b1[i]*b2[j]*b3[k]
    }
  }
}

现在的问题是 product 只存储 10 个值。如何解决这个索引问题？

Answer 1

要解决循环中的索引问题，您可以添加一个新变量，比如说counter，它会随着循环中的每次迭代而递增。

counter <- 1
for (i in 1:length(b1)) {
   for (j in 1:length(b2))  {
     for (k in 1:length(b3)) {
       product[counter] <- b1[i]*b2[j]*b3[k]
       counter = counter + 1
    }
  }
}

---

您也可以使用outer，它会给您相同的结果。

c(outer(b3, c(outer(b1, b2))))

如果有多个这样的向量将它们放在一个列表中并使用Reduce

c(Reduce(outer, list(b1, b2, b3)))

请注意，outer 期望将函数应用于向量的组合，但由于乘法 (*) 是默认函数，因此我们在此不再提及。