[Rd] This mail was found after problem with mailserver Oct 3rd - postmaster@hh.umu.se (PR#681)
maechler@stat.math.ethz.ch
maechler@stat.math.ethz.ch
Tue, 3 Oct 2000 19:09:33 +0200 (MET DST)
>>>>> "MM" == Martin Maechler <maechler@stat.math.ethz.ch> writes:
>>>>> "david" == david <david@orion-10.cs.byu.edu> writes:
david> When I issue the command: if ( c(2.8,3)[1] < pretty
david> (c(2.8,3))[1]) print ("problem")
david> the result is "problem" prints. I have traced it to bits and
david> the 2.8 of the returned vector from pretty is slightly larger.
david> This seems to be due to the conversion to base 10 with the ns
david> and unit variables in src/appl/pretty.c.
david> Anyway if I do something like: hist(c(2.8,3)) it fails as well
david> because it uses pretty.
david> Thanks for you attention, DavidMcLaughlin@byu.edu
MM> Thank you for the bug report. This bug is ``new'' (a bit more than
MM> a year!) and was introduced when someone fixed another problem in
MM> pretty.c by adding "1e-5" fuzz a bit too generously hmm,.......
MM> I'll have a look.
david> Version: .... R 1.1.1 ...
Ok, I've fixed pretty() for the development version
such that
p <- pretty(x)
all(p[1] <= x & x <= p[2])
should be guaranteed (again).
Note however, that this really also changes some "results",
that may be unexpected,
e.g.
R 1.1.1 :
---------
> data(women); (hT <- hist(women$weight, nc= 7, plot = FALSE))[2]
$counts
[1] 3 1 2 2 1 1 2 1 1 1
--------------
R-devel [to be 1.2.0]:
---------------------
> data(women); (hT <- hist(women$weight, nc= 7, plot = FALSE))[2]
$counts
[1] 2 2 2 1 2 1 2 1 1 1
>
------
Martin Maechler <maechler@stat.math.ethz.ch> http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum LEO D10 Leonhardstr. 27
ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND
phone: x-41-1-632-3408 fax: ...-1228 <><
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._