web-dev-qa-db-ja.com

Rで対称行列を作成する

Rには対称であるはずの行列がありますが、マシンの精度により、行列は決して対称ではありません（値は約10 ^ -16異なります）。私は行列が対称であることを知っているので、私はこれまで問題を回避するためにこれを行ってきました：

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

これに適した1行のコマンドはありますか？

rsymmetric
27
user2005253

値がその分だけ異なる場合、回避策は本当に必要ですか？

私の以前の答えは間違っていたと誰かが指摘しました。私は他のいくつかのほうが好きですが、これを削除することはできないので（去った​​ユーザーに受け入れられます）、micEconパッケージを使用する別のソリューションがあります：

symMatrix(s[upper.tri(s, TRUE)], nrow=nrow(s), byrow=TRUE)
10
Peyton
s<-matrix(1:25,5)
s[lower.tri(s)] = t(s)[lower.tri(s)]
52
user3318600

RのforceSymmetricパッケージのMatrix関数を使用して、マトリックスを強制的に対称にすることができます。

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
13
Metrics
 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
8
42-