It's not the compiler. pghyper() in SuppDists does the same thing. Its just rounding error. Set the result to 0 if it bothers you.

pghyper(24,514,53, 5961, lower.tail=F)
[1] -3.325965e-12


Roger D. Peng wrote:
Did you compile with gcc-2.96? I think there were some problems with the floating point arithmetic with that compiler (at least for the earlier versions released by Red Hat).

-roger

[EMAIL PROTECTED] wrote:

Hello,

I've come across the following situation in R-1.8.1 (compile + running under
RedHat 7.1):



phyper(24, 514, 5961-514, 53, lower.tail=T)


[1] 1

phyper(24, 514, 5961-514, 53, lower.tail=F)


[1] -1.037310e-11

I'd expect the later to be 0 or some very small positive number. Is this a
number under-flow of the calculation? Do you think I'm safe if I just set the
result to 0 in these cases?


kind regards,

Arne

______________________________________________
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html



______________________________________________
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html




--
Bob Wheeler --- http://www.bobwheeler.com/
        ECHIP, Inc. ---
Randomness comes in bunches.

______________________________________________
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Reply via email to