Создание симметричной матрицы в R

Question

У меня есть матрица в R, которая должна быть симметричной, однако из-за точности машины матрица никогда не симметрична (значения отличаются примерно на 10^-16). Так как я знаю, что матрица симметрична я делал это до сих пор, чтобы обойти эту проблему:

s.diag = diag(s) 
s[lower.tri(s,diag=T)] = 0 
s = s + t(s) + diag(s.diag,S) 

Есть ли лучшая команда одна строка для этого?

Answer 1

Это обходное решение действительно необходимо, если значения только сильно отличаются друг от друга?

Кто-то указал, что мой предыдущий ответ был неправильным. Мне нравятся некоторые из тех, что лучше, но так как я не могу удалить этот (принят пользователем, который оставил), вот еще одно решение, используя micEcon пакет:

symMatrix(s[upper.tri(s, TRUE)], nrow=nrow(s), byrow=TRUE) 

Answer 2

s<-matrix(1:25,5) 
pmean <- function(x,y) (x+y)/2 
s[] <- pmean(s, matrix(s, nrow(s), byrow=TRUE)) 
s 
#------- 
    [,1] [,2] [,3] [,4] [,5] 
[1,] 1 4 7 10 13 
[2,] 4 7 10 13 16 
[3,] 7 10 13 16 19 
[4,] 10 13 16 19 22 
[5,] 13 16 19 22 25 

Answer 3

Вы можете заставить матрица симметричной с использованием forceSymmetric функции в Matrix пакета в R:

library(Matrix) 
x<-Matrix(rnorm(9), 3) 
> x 
3 x 3 Matrix of class "dgeMatrix" 
      [,1]  [,2]  [,3] 
[1,] -1.3484514 -0.4460452 -0.2828216 
[2,] 0.7076883 -1.0411563 0.4324291 
[3,] -0.4108909 -0.3292247 -0.3076071 

A <- forceSymmetric(x) 
> A 
3 x 3 Matrix of class "dsyMatrix" 
      [,1]  [,2]  [,3] 
[1,] -1.3484514 -0.4460452 -0.2828216 
[2,] -0.4460452 -1.0411563 0.4324291 
[3,] -0.2828216 0.4324291 -0.3076071 

Answer 4

s<-matrix(1:25,5) 
s[lower.tri(s)] = t(s)[lower.tri(s)]