矩阵的逆矩阵

我想知道你推荐的方法来计算矩阵的逆矩阵?

我发现的方式似乎不令人满意。例如,

> 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)给出正确的逆。你的代码的问题是,你使用错误的运算符矩阵乘法。你应该使用solve(c)%*%c来调用R中的矩阵乘法。

R当你调用solve(c)* c时执行元素乘元。

http://stackoverflow.com/questions/11995832/inverse-of-matrix-in-r

本站文章除注明转载外,均为本站原创或编译
转载请明显位置注明出处:矩阵的逆矩阵