[Rd] Bug in qbinom? (PR#511)
Fri, 7 Apr 2000 18:23:38 +0200 (MET DST)
>>>>> "TL" == Thomas Lumley <firstname.lastname@example.org> writes:
TL> On Fri, 7 Apr 2000 email@example.com wrote:
>> n_10;p_0.5;jjx_0:n;qbinom(pbinom(jjx,n,p),n,p) # This one works as
>> n_100;p_0.5;jjx_0:n;qbinom(pbinom(jjx,n,p),n,p) # This one causes
>> severe problems
TL> The problem is in qbinom(1-epsilon,100,0.5) for epsilon between
TL> 10^-15 and 10^-17. It appears to go for ever.
TL>  -1.110223e-16
TL> An even simpler version is
TL> Presumably a check for zero needs some fuzz.
Note that we (I) also have a considerable improvement for dbinom() in the
pipeline, from Clive Loader's post of a few weeks ago.
But this may be a completely orthogonal issue...
I will have a look over the weekend.
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: firstname.lastname@example.org