2

正方対称の正定行列を逆にしたい。Rには2つの関数があることは知っていますが、結果は異なりsolve()ます。chol2inv()なぜこれが起こるのかを知る必要がありますか？

ありがとうございました。

4

2 に答える 2

6

とを含む、逆行列を計算するいくつかの方法がsolve()ありchol2inv()ます。

> A <- matrix(c(2, -1, 0, -1, 2, -1, 0, -1, 2), 3)

> solve(A)
     [,1] [,2] [,3]
[1,] 0.75  0.5 0.25
[2,] 0.50  1.0 0.50
[3,] 0.25  0.5 0.75

> chol2inv(chol(A))
     [,1] [,2] [,3]
[1,] 0.75  0.5 0.25
[2,] 0.50  1.0 0.50
[3,] 0.25  0.5 0.75

> library(MASS)
> ginv(A)
     [,1] [,2] [,3]
[1,] 0.75  0.5 0.25
[2,] 0.50  1.0 0.50
[3,] 0.25  0.5 0.75
于 2013-03-11T10:30:38.080 に答える
5

solve元の行列を与える必要がありますが、chol2inv事前に計算されたコレスキー分解を使用する場合：

set.seed(1)
a<-crossprod(matrix(rnorm(9),3,3))
a_chol<-chol(a)
solve(a)
            [,1]        [,2]       [,3]
[1,]  1.34638151 -0.02957435  0.8010735
[2,] -0.02957435  0.32780020 -0.1786295
[3,]  0.80107345 -0.17862950  1.4533671
chol2inv(a_chol)
            [,1]        [,2]       [,3]
[1,]  1.34638151 -0.02957435  0.8010735
[2,] -0.02957435  0.32780020 -0.1786295
[3,]  0.80107345 -0.17862950  1.4533671
于 2013-03-11T10:27:41.390 に答える