Rの逆行列

Question

行列の逆行列を計算するにはどうすればよいでしょうか？

私が見つけた方法は満足のいくものではないようです。たとえば、

> c=rbind(c(1, -1/4), c(-1/4, 1)) 
> c 
     [,1] [,2] 
[1,] 1.00 -0.25 
[2,] -0.25 1.00 
> inv(c) 
Error: could not find function "inv" 
> solve(c)  
      [,1]  [,2] 
[1,] 1.0666667 0.2666667 
[2,] 0.2666667 1.0666667 
> solve(c)*c 
      [,1]  [,2] 
[1,] 1.06666667 -0.06666667 
[2,] -0.06666667 1.06666667 
> qr.solve(c)*c 
      [,1]  [,2] 
[1,] 1.06666667 -0.06666667 
[2,] -0.06666667 1.06666667 

ありがとうございます！

Answer 1

solve(c)は、正しい逆行列を与えます。あなたのコードの問題は、行列の乗算に間違った演算子を使用していることです。 を使用してRで行列乗算を呼び出す必要があります。

Rは、solve(c) * cを呼び出すときに要素の乗算を実行します。

Answer 2

あなたはMASSパッケージに

Answer 3

注意を機能ginv（）（ムーア・ペンローズ一般逆）を使用することができますが、スピードを気にして特異点を心配する必要がない場合は、solve()をそのあなたがチェックすることができますよう、それは、はるかに高速であるため、ginv()に優先されるべきである。行列表記で

require(MASS) 
mat <- matrix(rnorm(1e6),nrow=1e3,ncol=1e3) 

t0 <- proc.time() 
inv0 <- ginv(mat) 
proc.time() - t0 

t1 <- proc.time() 
inv1 <- solve(mat) 
proc.time() - t1 

Answer 4

それはオペレータ "大きな違いになり* "、演算子" %*% "となります。最初は要素ごとに乗算を行い、2つ目は行列乗算の正しい式です。 何HOUがやるべきことは次のとおりです。

c = rbind(c(1, -1/4), c(-1/4, 1)) 

solve(c) %*% c