Mir*_*Mir 2 r symmetric matrix
我正在尝试使用 R 制作一个对称矩阵。我已经有了一个矩阵。我的矩阵非常大,所以下面是一个简单的例子。
前任。
1 2 3
4 5 6
7 8 9
我需要让它们像这样。
1 2+4 3+7
4+2 5 6+8
7+3 8+6 9
//所以我尝试了这个。// mat 是我正在使用的矩阵。
lowervector <- square_07[lower.tri(square_07, diag = FALSE)]
uppervector <- square_07[upper.tri(square_07, diag = FALSE)]
lowermat <- square_07
uppermat <- square_07
lowermat[lower.tri(lowermat, diag = FALSE)] <- t(square_07)[lower.tri(square_07, diag = FALSE)]
uppermat[upper.tri(uppermat, diag = FALSE)] <- t(square_07)[upper.tri(square_07, diag = FALSE)]
当我执行最后两行时,出现错误;
下标 'upper.tri(uppermat, diag = FALSE)' 是一个矩阵,数据 't[upper.tri(square_07, diag = FALSE)]' 的大小必须为 1。
你应该知道。上面的矩阵只是一个例子。我的实际矩阵要大得多。它是一个248*248的矩阵。
我怎么解决这个问题?
这更简单:
result <- square_07 + t(square_07)
diag(result) <- diag(square_07)
result
# [,1] [,2] [,3]
# [1,] 1 6 10
# [2,] 6 5 14
# [3,] 10 14 9
Run Code Online (Sandbox Code Playgroud)