R中的广义逆

Question

我可以使用ginv函数从MASS库中得到一个矩阵的摩尔-彭罗斯Generalisied逆.

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

在SAS中,我们确实有多个函数来获得矩阵的广义逆.SVD可用于找到广义逆,但这也是Moore-Penrose.我想知道在R中是否有任何函数可以得到除Moore-Penrose Generalisied Inverse之外的矩阵的广义逆(这不是唯一的).在此先感谢您的帮助和时间.

编辑

矩阵A的广义逆定义为满足等式AGA = A的任何矩阵G.

这个G不是Moore-Penrose Generalisied Inverse,所以它不是唯一的.

Answer 1

大多数情况下,您并不真正想要矩阵的逆矩阵,因为最终结果可能会在您完成时通过舍入误差来破坏.

使用部分旋转和缩放创建LU分解更为典型.使用它在右侧向量上执行向前/向后替换以获得解决方案.如果您有多个RHS向量,这将特别有用,因为您可以重复应用它.

您需要Matrix包才能执行此操作.