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