Не могу решить эту матрицу

Question

Может кто-нибудь помочь, пожалуйста?Невозможно решить эту матрицу, используя R. Код воспроизводим.

(a <- matrix(c(1, 0, 1, 0, 10, 1, 0, 0, 1, 20, 0, 1, 1, 0, 5, 0, 1, 0, 1, 10), nrow = 5))
(b <- matrix(c(20, 15, 10, 25, 475), nrow = 5))
solve(a) %*% b

Answer 1

У вас есть неквадратная матрица.Поскольку в нем больше строк, чем столбцов, вы можете использовать OLS:

> lm.fit(a, b)$coefficients
x1 x2 x3 x4 
 5 15  5 10 

В качестве альтернативы вы можете использовать обобщенный обратный, например, из MASS:

> MASS::ginv(a) %*% b
     [,1]
[1,]    5
[2,]   15
[3,]    5
[4,]   10