R /从下三角元素的向量中创建对称矩阵

时间:2016-03-01 10:57:43

标签: r matrix vector

我在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 

1 个答案:

答案 0 :(得分:7)

我们可以使用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