Question

我有一个矢量说

a = c(1,2,3,4,5,6)

我想将它们组织成一个上三角矩阵（不考虑对角元素，它们都是零）按行。我的目标是获得以下矩阵：

     [,1] [,2] [,3] [,4]
[1,]    0    1    2    3
[2,]    0    0    4    5
[3,]    0    0    0    6
[4,]    0    0    0    0

但我采用以下方式是用这个向量替换对角元素，但是按列分配值。例如，

b= matrix(0, 4, 4)
b[upper.tri(b, diag=FALSE)]=a

它会给我以下矩阵

   [,1] [,2] [,3] [,4]
[1,]    0    1    2    4
[2,]    0    0    3    5
[3,]    0    0    0    6
[4,]    0    0    0    0

原因是当R为矩阵赋值时，默认情况下，它会按列分配。我想知道是否有一种简单的方法可以在不编写for循环的情况下解决我的问题。

我之前发现了与我的问题相关的类似帖子，但它没有解释按行分配上三角矩阵的值：

creating a triangular matrix

提前致谢！

Answer 1

这是一个选项

b[lower.tri(b, diag=FALSE)] <- a
b <- t(b)
b
#      [,1] [,2] [,3] [,4]
# [1,]    0    1    2    3
# [2,]    0    0    4    5
# [3,]    0    0    0    6
# [4,]    0    0    0    0

或者，根据需要重新排序a并将其分配到右上角三角形中：

ut <- upper.tri(b, diag=FALSE)
b[ut] <- a[order(row(ut)[ut], col(ut)[ut])]
b
     [,1] [,2] [,3] [,4]
[1,]    0    1    2    3
[2,]    0    0    4    5
[3,]    0    0    0    6
[4,]    0    0    0    0

Answer 2

请注意，要填充ASYMMETRIC矩阵，您可以先通过上面显示的代码填充上三角形，然后使用不同的矢量填充下部（不需要转置）。

  c <- c(7,8,9,10,11,12)
  b[lower.tri(b, diag=FALSE)] <- c

R中行的上三角矩阵的向量

2 个答案: