[R] Odp: Decimal Accuracy Loss?
petr.pikal at precheza.cz
Wed Apr 6 17:43:48 CEST 2011
r-help-bounces at r-project.org napsal dne 06.04.2011 17:33:48:
> This is hopefully a quick question on decimal accuracy. Is any
> decimal accuracy lost when casting a numeric vector as a matrix? And
> then again casting the result back to a numeric?
> I'm finding that my calculation values are different when I run for
> loops that manually calculate matrix multiplication as compared to
> when I cast the vectors as matrices and multiply them using "%*%".
> (The errors are very small, but the process is run iteratively
> thousands of times, at which point the error between the two
> differences becomes noticeable.)
> I've read FAQ # 7.31 "Why doesn't R think these numbers are equal?",
> but just want to confirm that the differences in values are due to
> differences in the matrix multiplication operator and manual
> calculation via for loops, rather than information that is lost when
> casting a numeric as a matrix and back again.
Without some example it is difficult to see the possible sources of
difference. Clever people may know how %*% operator really works, but only
those who are able mind reading can know what you do inside your for
> Thanks in advance for the help,
> R-help at r-project.org mailing list
> PLEASE do read the posting guide
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help