Hi, I am working with large numbers and identified that R looses precision for such high numbers.
The precision is lost exactly when the number is equal or larger than 53 bits. See the following output which shows that the numbers below 53 bit have proper precision: > 2^53 [1] 9007199254740992 > 2^53-1 [1] 9007199254740991 > 2^53-2 [1] 9007199254740990 Now, see the numbers above 53 bit: > 2^53 [1] 9007199254740992 > 2^53+1 [1] 9007199254740992 > 2^53+2 [1] 9007199254740994 > 2^53+3 [1] 9007199254740996 > 2^53+4 [1] 9007199254740996 Is there a solution to the problem? Thanks a lot Stephan ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
