Matrix inverse
On Fri, 28 Apr 2000, Thomas Lumley wrote:
On Fri, 28 Apr 2000, Patrik Waldmann wrote:
I haven't found a function that directly calculates the matrix inverse, does it exist? Otherwise what would be the fastest way to do it?
Ainv<-solve(A) Note that if you want to solve Ax=b you can do solve(A,b), which should be quicker than solve(A)%*%b
How about 'Ainv <- qr.solve(A)'? I happened to read the help page for 'qr.solve' the other day, and there I found that qr.solve(A, b) "is == but much better than solve(A) %*% b". (I guess that 'better than' refers to numerical stability?) Is qr.solve generally to be preferred to solve, which seems to be indicated by the cited help page? And what is hidden behind solve? G?ran
G?ran Brostr?m, professor tel: +46 90 786 5223 Department of Statistics fax: +46 90 786 6614 Ume? University SE-90187 Ume?, Sweden e-mail: gb at stat.umu.se -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._