2012-02-03 20 views
0

MASSライブラリのginv関数を使用して、行列のMoore-Penrose一般逆行列を得ることができます。一般化された逆行列R

m <- matrix(1:9, 3, 3) 
library(MASS) 
ginv(m) 

SASでは、行列の一般化逆行列を得るために複数の関数があります。 SVDを使って一般化された逆行列を見つけることができますが、これもMoore-Penroseです。 Moore-Penrose Generalisied Inverse以外の行列(これはユニークではない)の一般化された逆行列を得るためにRに関数が存在するのだろうかと思います。あなたの助けと時間を事前に感謝します。

編集行列Aの一般化逆行列は、任意のGとして定義される

その

このGがそうムーア・ペンローズ逆Generalisiedない 満たす方程式AGA = A.ユニークではありません。

+1

他の疑似逆行列はどのようにしますか?あなたはSASから模倣したいものをより正確に表現できますか? – Iterator

+0

Rのヘルプリストで頻繁に見られる署名を引用するには、「あなたは解決しようとしている問題は何ですか? –

答えて

6

ほとんどの場合、行列の逆数は必要ありません。最終結果は、完了した時点で誤差を丸めることによって損なわれる可能性があるからです。

部分的なピボットとスケーリングを使用してLU分解を作成することがより一般的です。解決策を得るために、右手側のベクトルで前方/後方の置換を実行するために使用します。これは、繰り返し適用することができるため、複数のRHSベクタがある場合に特に役立ちます。

これを行うには、Matrixパッケージが必要です。

+0

コメントありがとうございます。例を挙げますか?ありがとう – MYaseen208

1

はい、本当に大変なご迷惑をおかけしています。Rパッケージはもう使用できません。あるいは、pracmaパッケージを使用することもできます。 あなたのムーア・ペンローズ一般化逆:

pinv(m) 
関連する問題