[R] Matrix multiplication in R is inaccurate!!

Vikas Rawal vikas at mail.jnu.ac.in
Fri Jan 28 06:56:02 CET 2005


If you multiply a matrix by its inverse you should get an identity matrix. In R, you get an answer that is accurate up to about 16 decimal points? Why can't one get a perfect answer?

See for example: 

c(5,3)->x1
c(3,2)->x2
cbind(x1,x2)->x
solve(x)->y
x%*%y

Vikas Rawal

==============================================

 This Mail was Scanned for Virus and found Virus free




More information about the R-help mailing list