R中的廣義反向

我可以使用ginv函數從MASS庫中獲得矩陣的Moore-Penrose Generalisied Inverse。R中的廣義反向

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

在SAS中，我們確實有多個函數來獲得矩陣的廣義逆。 SVD可以用來找到廣義逆，但這又是一個Moore-Penrose。我想知道除了Moore-Penrose Generalisied Inverse之外，R中是否有任何函數可以得到矩陣（它不是唯一的）的廣義逆。預先感謝您的幫助和時間。

編輯

一個矩陣A的廣義逆被定義爲任何矩陣G的是滿足方程式AGA = A.

這G不是摩爾-Penrose廣義逆Generalisied所以它不是獨一無二的。

2012-02-03 MYaseen208

你會喜歡什麼其他的僞反向？你能更準確地從SAS中模擬出什麼嗎？ – Iterator 2012-02-03 03:36:25

引用R幫助列表中經常出現的簽名：「您嘗試解決什麼問題？」 – 2012-02-03 14:49:03

大多數時候你並不真正想要矩陣的逆矩陣，因爲最終的結果可能會在你完成時被舍入誤差所破壞。

使用部分旋轉和縮放創建LU分解更爲典型。使用它來執行右側向量的前向/後向替換以獲得解決方案。如果您有多個RHS向量，這是特別有用的，因爲您可以重複應用它。

您需要Matrix包才能做到這一點。

2012-02-03 00:54:07 duffymo

感謝您的評論。你想提供一個例子嗎？謝謝 – MYaseen208 2012-02-03 00:57:42

是的，這是一個很大的不便R包不再可用。或者，您可以使用pracma包。和你的摩爾 - 彭羅斯Generalisied逆：

pinv(m)

2014-11-03 19:17:18

回答