[Rd] Re: [R] Bug in qchisq? (PR#875)

p.dalgaard@biostat.ku.dk p.dalgaard@biostat.ku.dk
Sat, 10 Mar 2001 02:24:18 +0100 (MET)


Kenneth Cabrera <krcabrer@epm.net.co> writes:

> Hello developers and users:
> 
> My system fails (the computer freezes) when I use the ncp  parameter,
> with the lower.tail=FALSE option in the qchisq function.
> 
> qchisq(0.025,31,ncp=1,lower.tail=FALSE)

Yup, that's a bug. We have in pnchisq.c

48          for (ux = 1.0; pnchisq(ux, n, lambda, lower_tail, log_p) <
p; ux *= 2);
49          for (lx = ux;  pnchisq(lx, n, lambda, lower_tail, log_p) >
p; lx *= 0.5);

but if we look at the opposite tail, we also need to reverse that
logic since pnchisq(x,...) is now decreasing in x. Otherwise the
algorithm will get stuck in one of the for loops.

(Mailed to r-bugs)
-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk)             FAX: (+45) 35327907

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._