循环通过R中的矩阵

Question

我尝试通过矩阵循环，但找不到一种简单而优雅的方法，而不是编写许多（> 10）方程......有人可以帮我吗？

我的矩阵看起来像这样：

我想计算以下内容：

(0 * 0 * 4/24) + (0 * 1 * 6/24) + (0 * 2 * 3/24) + (1 * 0 * 3/24) + (1 * 1 * 4/24) + (1 * 2 * 4/24)

而不是使用

__ 顺便说一句：我的矩阵代码

vals<- c(4/24, 6/24, 3/24, 3/24, 4/24, 4/24)
x <- c(0,1)
y <- c(0,1,2)
df <- matrix(vals, byrow = TRUE, nrow = 2, ncol = 3,
             dimnames = list(x,y))

而不是手动计算每个步骤，我认为应该有一个for循环方法，但无法弄清楚..

Answer 1

一个可能的解决方案：

c(x %*% df %*% y) 

#> [1] 0.5

---

另一种可能的解决方案，基于outer ：

sum(outer(x, y, Vectorize(\(x,y) x*y*df[x+1,y+1])))

#> [1] 0.5

Answer 2

x <- c(0, 1)
y <- c(0, 1, 2)
vals<- c(4/24, 6/24, 3/24, 3/24, 4/24, 4/24)
mat <- matrix(vals, byrow = TRUE, nrow = 2, ncol = 3,
              dimnames = list(x,y))  ## not a data frame; don't call it "df"

还有比for循环更好的方法：

sum(tcrossprod(x, y) * mat)
#[1] 0.5

Answer 3

sum((x %o% y) * df)

解释：

x %o% y得到向量x和y的外积，即：

#>      [,1] [,2] [,3]
#> [1,]    0    0    0
#> [2,]    0    1    2

由于它与df具有相同的维度，因此您可以将相应的元素相乘并得到总和： sum((x %o% y) * df)

Answer 4

如果您是 R 新手（就像我一样），这里是循环方法。

result = 0
for (i in 1:length(x)) {
  for (j in 1:length(y)) {
    result = result + x[i] * y[j] * df[i, j]
  }
}
result