求和后将 1xn + 1xn 矩阵 reshape 为对称矩阵

Question

例如，我有两个尺寸为 1x5 的矩阵，我想在执行 c = a+b 后将它们 reshape 为对称的 5x5 矩阵。操作然后添加 diag(c) <- 1

举个例子:

a <- matrix(seq(1,5), byrow = T)
b <- matrix(seq(1,5), byrow = T)

我正在寻找的结果是:

1 3 4 5 6
3 1 5 6 7
4 5 1 7 8
5 6 7 1 9
6 7 8 9 1

请帮忙，提前谢谢你

Answer 1

我们可以使用outer 来做sum，然后将对角线元素赋给1

out <- outer(a[,1], b[,1], FUN = `+`)
diag(out) <- 1
out
#     [,1] [,2] [,3] [,4] [,5]
#[1,]    1    3    4    5    6
#[2,]    3    1    5    6    7
#[3,]    4    5    1    7    8
#[4,]    5    6    7    1    9
#[5,]    6    7    8    9    1