我有一个矩阵的下三角，我试图将其转换为 dissim 矩阵，因此它需要是对称的。
print(rdf)
X0 X1 X2 X3 X4
0 0.0000000 NA NA NA NA
1 0.5340909 0.0000000 NA NA NA
2 0.5340909 0.0000000 0.0000000 NA NA
3 0.3200000 0.5227273 0.5227273 0.0000000 NA
4 0.6263736 0.4945055 0.4945055 0.5384615 0
library(gdata)
upperTriangle(rdf) <- lowerTriangle(rdf)
isSymmetric(rdf)
X0 X1 X2 X3 X4
0 0.0000000 0.5340909 0.5340909 0.6263736 0.4945055
1 0.5340909 0.0000000 0.3200000 0.0000000 0.5227273
2 0.5340909 0.0000000 0.0000000 0.5227273 0.4945055
3 0.3200000 0.5227273 0.5227273 0.0000000 0.5384615
4 0.6263736 0.4945055 0.4945055 0.5384615 0.0000000
[1] FALSE
到底是怎么回事？注意，我不需要使用 gdata
您需要确保要复制的元素顺序正确：
m <- matrix(NA,4,4)
m[lower.tri(m,diag=TRUE)] <- 1:10
[,1] [,2] [,3] [,4]
[1,] 1 NA NA NA
[2,] 2 5 NA NA
[3,] 3 6 8 NA
[4,] 4 7 9 10
makeSymm <- function(m) {
m[upper.tri(m)] <- t(m)[upper.tri(m)]
return(m)
}
makeSymm(m)
或者您可以使用内置的
Matrix::forceSymmetric(m,uplo="L")
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
