Đảo ngược ma trận trong R

Question

Tôi đã tự hỏi cách bạn đề xuất để tính toán nghịch đảo của ma trận là gì?

Cách tôi thấy dường như không thỏa đáng. Ví dụ:

> c=rbind(c(1, -1/4), c(-1/4, 1)) 
> c 
     [,1] [,2] 
[1,] 1.00 -0.25 
[2,] -0.25 1.00 
> inv(c) 
Error: could not find function "inv" 
> solve(c)  
      [,1]  [,2] 
[1,] 1.0666667 0.2666667 
[2,] 0.2666667 1.0666667 
> solve(c)*c 
      [,1]  [,2] 
[1,] 1.06666667 -0.06666667 
[2,] -0.06666667 1.06666667 
> qr.solve(c)*c 
      [,1]  [,2] 
[1,] 1.06666667 -0.06666667 
[2,] -0.06666667 1.06666667 

Cảm ơn!

Answer 1

solve(c) cung cấp sai chính xác. Vấn đề với mã của bạn là bạn đang sử dụng toán tử sai cho phép nhân ma trận. Bạn nên sử dụng solve(c) %*% c để gọi phép nhân ma trận trong R.

R thực hiện phần tử theo phép nhân tố khi bạn gọi solve(c) * c.

Answer 2

Bạn có thể sử dụng chức năng ginv() (Moore-Penrose tổng quát nghịch đảo) trong gói MASS

Answer 3

Lưu ý rằng nếu bạn quan tâm đến tốc độ và không cần phải lo lắng về kỳ dị, solve() nên được ưa thích để ginv() vì nó là nhanh hơn nhiều, như bạn có thể kiểm tra:

require(MASS) 
mat <- matrix(rnorm(1e6),nrow=1e3,ncol=1e3) 

t0 <- proc.time() 
inv0 <- ginv(mat) 
proc.time() - t0 

t1 <- proc.time() 
inv1 <- solve(mat) 
proc.time() - t1 

Answer 4

Trong ký hiệu ma trận nó làm cho một sự khác biệt lớn các nhà điều hành "* "và toán tử" %*% ". Phần tử nhân tố đầu tiên theo phần tử, phần tử thứ hai là công thức chính xác cho phép nhân ma trận. Điều bạn nên làm là:

c = rbind(c(1, -1/4), c(-1/4, 1)) 

solve(c) %*% c