[R-SIG-Mac] Unexpected output from round()

Duncan Murdoch murdoch.duncan at gmail.com
Sat Jun 13 04:05:46 CEST 2015


On 12/06/2015 8:53 PM, John Maindonald wrote:
> I cannot reproduce it (R version 3.2.0 Patched (2015-04-19 r68206)) in a terminal window:

If you set options(digits=16) you'll get it.  It's just FAQ 7.31, and
the fact that R only stores 15 or 16 digit accuracy.

Duncan Murdoch

> 
>> dput(pi)
> 3.14159265358979
>> round(dput(pi),2)
> 3.14159265358979
> [1] 3.14
>> Cl <- c(82.9888290534313, 100.544763749464)
>> dput(Cl)
> c(82.9888290534313, 100.544763749464)
>> round(dput(Cl),2)
> c(82.9888290534313, 100.544763749464)
> [1]  82.99 100.54
>> sessionInfo()
> R version 3.2.0 Patched (2015-04-19 r68206)
> Platform: x86_64-apple-darwin13.4.0 (64-bit)
> Running under: OS X 10.10.3 (Yosemite)
> 
> locale:
> [1] en_NZ.UTF-8/en_NZ.UTF-8/en_NZ.UTF-8/C/en_NZ.UTF-8/en_NZ.UTF-8
> 
> attached base packages:
> [1] stats     graphics  grDevices utils     datasets  methods   base  
> 
> John Maindonald             email: john.maindonald at anu.edu.au
> 
>> On 13/06/2015, at 12:39, Dennis Fisher <fisher at plessthan.com> wrote:
>>
>> R 3.2.0
>> OS X
>>
>> Colleagues
>>
>> I received unexpected output from a command issued in a terminal window:
>>
>>> dput(CI)
>> c(82.9888290534313, 100.544763749464)
>>
>>> round(dput(CI), 2)
>> c(82.9888290534313, 100.544763749464)
>> [1]  82.98999999999999 100.54000000000001
>>
>> I expected the result of the second command to be:
>> 	82.99   100.54
>>
>> Interesting, the problem also occurs if I round to 5 or 8 digits.  
>>
>> Of note, the problem does NOT occur in R GUI or in Windows.
>>
>> I am aware of FAQ 7.31 (Why doesn’t R think these numbers are equal?) but I don’t see how this could be manifestation of that issue.
>>
>> I am curious as to the explanation of this issue.
>>
>> Dennis
>>
>>
>> Dennis Fisher MD
>> P < (The "P Less Than" Company)
>> Phone: 1-866-PLessThan (1-866-753-7784)
>> Fax: 1-866-PLessThan (1-866-753-7784)
>> www.PLessThan.com
>>
>> _______________________________________________
>> R-SIG-Mac mailing list
>> R-SIG-Mac at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-mac
> 
> _______________________________________________
> R-SIG-Mac mailing list
> R-SIG-Mac at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-mac
>



More information about the R-SIG-Mac mailing list