[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