Nghịch đảo của ma trận trong R

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

Những cách tôi thấy có vẻ không khả quan. 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!

solve(c)không đưa ra nghịch đảo chính xác. Vấn đề với mã của bạn là bạn đang sử dụng sai toán tử 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ép nhân phần tử khi bạn gọi solve(c) * c.

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

Lưu ý rằng nếu bạn quan tâm đến tốc độ và không cần lo lắng về các điểm kỳ dị, solve()thì nên ưu tiên hơn ginv()vì nó nhanh hơn nhiều, vì 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