我被提供了一些数据以向量形式在 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/

10-12 22:32