Tạo ma trận đối xứng trong R

Question

Tôi có ma trận trong R được cho là đối xứng, tuy nhiên, do độ chính xác của máy, ma trận không bao giờ đối xứng (các giá trị khác nhau khoảng 10^-16). Vì tôi biết ma trận là đối xứng nên tôi đã làm điều này cho đến nay để giải quyết vấn đề:

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

Có lệnh nào tốt hơn cho điều này không?

Answer 1

Giải pháp thay thế có thực sự cần thiết không nếu giá trị chỉ khác nhau nhiều?

Ai đó đã chỉ ra rằng câu trả lời trước của tôi là sai. Tôi thích một số những cái khác tốt hơn, nhưng kể từ khi tôi không thể xóa cái này (chấp nhận bởi một người dùng rời), đây là thêm một giải pháp sử dụng các micEcon gói:

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

Bạn có thể buộc các ma trận đối xứng sử dụng forceSymmetric chức năng trong Matrix gói vào 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)]