将下三角元素向量创建成对称矩阵

Question

我得到了一些用R语言处理的向量形式的数据。这个大小为n（n + 1）/2的向量的元素应该被重新排列成对称矩阵。

例子：

n<-4
x<-seq(from=1,to=n*(n+1)/2) 

如何优雅地获得

mat<- 1 2 3 4
      2 5 6 7
      3 6 8 9  
      4 7 9 10 

- Stefan Voigt

Answer 1

我们可以使用lower.tri。

 m1 <- matrix(, n, n)
 m1[lower.tri(m1, diag=TRUE)] <- x
 m2 <- t(m1)
 m2[lower.tri(m2, diag=TRUE)] <- x

或者，您可以跳过最后两个步骤。

 pmax(m1, t(m1), na.rm=TRUE)
 #      [,1] [,2] [,3] [,4]
 #[1,]    1    2    3    4
 #[2,]    2    5    6    7
 #[3,]    3    6    8    9
 #[4,]    4    7    9   10

- akrun