Обобщенный обратный в R

Question

Я могу использовать функцию ginv из библиотеки MASS, чтобы получить обобщенную обратную матрицу Мура-Пенроуза.

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

В SAS у нас есть несколько функций для получения обобщенной обратной величиныматрица.СВД можно использовать для поиска обобщенного обратного, но опять-таки это Мур-Пенроуз.Интересно, есть ли в R какая-либо функция для получения обобщенной обратной матрицы (которая не является уникальной), отличной от обобщенной обратной задачи Мура-Пенроуза.Заранее благодарим за помощь и время.

Редактировать

Обобщенная обратная матрица A определяется как любая матрица G, удовлетворяющая уравнению AGA =A.

Этот G не является обобщенным обратным преобразованием Мура-Пенроуза, поэтому он не уникален.

Answer 1

В большинстве случаев вам не нужно, чтобы матрица была инвертирована, потому что конечный результат может быть испорчен ошибками округления к тому времени, когда вы закончите.декомпозиция с использованием частичного поворота и масштабирования.Используйте его, чтобы выполнить прямое / обратное замещение вектора с правой стороны, чтобы получить решение.Это особенно полезно, если у вас есть несколько векторов RHS, потому что вы можете применять его повторно.

Для этого необходим пакет Matrix .

Answer 2

Да, правда, это большое неудобство. Пакеты R больше не доступны.В качестве альтернативы вы можете использовать пакет pracma.И ваш Обобщенный обратный Мур-Пенроуз:

pinv(m)