3

我有一个矩阵如下;

   1  2  3  4  5  
1  0  1  1  0  0  
2  0  0  1  1  0  
3  1  0  0  0  1  
4  0  0  1  0  0   
5  0  1  0  0  0  

我需要得到对称矩阵如下,对角线为 0;

   1  2  3  4  5  
1  0  1  1  0  0  
2  1  0  1  1  1  
3  1  1  0  1  1  
4  0  1  1  0  0   
5  0  1  1  0  0  

我用for循环试过了,当矩阵的大小很小的时候很好。但是对于大尺寸,这需要时间。在 R 中有什么有效的方法吗?

编辑:第一行在第二和第三列有 1,第一列在第三行有 1。但是要具有对称的第一列也需要在第二行中有 1 个。同样,每一行和每一列都应该考虑 1 来保持对称矩阵。

4

3 回答 3

5

使用逻辑运算

sym_mat <- (mat | t(mat)) * 1

@chinsoon12 的建议非常好。

sym_mat <- +(mat | t(mat))
于 2018-04-09T02:35:00.680 回答
2
mat=read.table(text="   1  2  3  4  5  
1  0  1  1  0  0  
           2  0  0  1  1  0  
           3  1  0  0  0  1  
           4  0  0  1  0  0   
           5  0  1  0  0  0  ")

mat=as.matrix(mat)
m=ifelse(!t(mat)[lower.tri(mat)], mat[lower.tri(mat)], t(mat)[lower.tri(mat)])
mat[lower.tri(mat)]=m
mat[upper.tri(mat)]=t(mat)[upper.tri(mat)]


 mat
  X1 X2 X3 X4 X5
1  0  1  1  0  0
2  1  0  1  1  1
3  1  1  0  1  1
4  0  1  1  0  0
5  0  1  1  0  0
于 2018-04-09T02:28:05.887 回答
1

怎么样:

myTrues <- which(mat == 1, arr.ind = T)
mat[rbind(myTrues, myTrues[,c(2,1)])] <- 1L

mat
     [,1] [,2] [,3] [,4] [,5]
[1,]    0    1    1    0    0
[2,]    1    0    1    1    1
[3,]    1    1    0    1    1
[4,]    0    1    1    0    0
[5,]    0    1    1    0    0
于 2018-04-09T02:35:05.233 回答