我被提供了一些数据以向量形式在 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
最佳答案
我们可以使用 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
关于R/从下三角形元素的向量中创建一个对称矩阵,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/35721068/