3

我有两个向量，我想要一个矩阵，其中元素是向量 1 的每个元素和向量 2 的每个元素的总和。

例如，这个矩阵第一行的第一个元素是向量1的第一个元素和向量2的第一个元素之和；第一行的第二个元素是向量 1 的第一个元素和向量 2 的第二个元素之和，依此类推。

例如，使用这两个向量

u <- c(1,2,3)
v <- c(4,5,6)

期望的结果是：

#       [,1] [,2] [,3]
# [1,]    5    6    7
# [2,]    6    7    8
# [3,]    7    8    9

我试过的：

A <- matrix( c(1:6), 3, 3 )

for(i in 1:3)
{
   for(j in 1:3)
   {
      A[j][i] <- u[i]+v[j]
   }
}

但我收到一些警告：

Warning messages:
1: In A[j][i] <- u[i] + v[j] :
  number of items to replace is not a multiple of replacement length
2: In A[j][i] <- u[i] + v[j] :
  number of items to replace is not a multiple of replacement length
3: In A[j][i] <- u[i] + v[j] :
  number of items to replace is not a multiple of replacement length
4: In A[j][i] <- u[i] + v[j] :
  number of items to replace is not a multiple of replacement length
5: In A[j][i] <- u[i] + v[j] :
  number of items to replace is not a multiple of replacement length
6: In A[j][i] <- u[i] + v[j] :
  number of items to replace is not a multiple of replacement length

有谁能够帮助我？

4

2 回答 2

7

这就是你要做的事情（注意矩阵子集不是两个括号，而是逗号分隔）：

u <- c(1,2,3)
v <- c(4,5,6)

A <- matrix( c(1:6), 3, 3 )

for(i in 1:3)
{
   for(j in 1:3)
   {
      A[i,j] <- u[i]+v[j]
   }
}

但这不是了解 R 的人会接近它的方式。一般来说，在 R 中，有比嵌套 for 循环更好的方法来做事。另一种方法是：

A <- outer(u,v,`+`)
于 2012-09-13T13:48:02.637 回答
1

我们也可以使用sapply

sapply(u, `+`, v)
于 2020-09-06T19:33:00.400 回答