Full_Name: Drouilhet R?my Version: 1.8.1 OS: Linux Submission from: (NULL) (195.221.43.136) qchisq(1,10) works well but qchisq(1,10,ncp=0) does not work whereas ncp=0 is the default value of the function qchisq(1,10). (of course, 10 will be replaced by any integer value). Let us notice that this bug occurs only when applying probability one. (qchisq(seq(0,.9,.1),10,ncp=0) works very well).
>>>>> "Remy" == Remy Drouilhet <Remy.Drouilhet@upmf-grenoble.fr> >>>>> on Mon, 19 Jan 2004 18:10:04 +0100 (CET) writes:Remy> qchisq(1,10) works well but qchisq(1,10,ncp=0) does Remy> not work whereas ncp=0 is the default value of the Remy> function qchisq(1,10). (of course, 10 will be replaced Remy> by any integer value). Remy> Let us notice that this bug occurs only when applying probability one. Remy> (qchisq(seq(0,.9,.1),10,ncp=0) works very well). Merci, Remy, for the report. Note that I've been working on improvements to "pnchisq" (C) and "qnchisq" for R-devel for several days now. You will see fixes there pretty soon, I hope. Regards, Martin Maechler <maechler@stat.math.ethz.ch> http://stat.ethz.ch/~maechler/ Seminar fuer Statistik, ETH-Zentrum LEO C16 Leonhardstr. 27 ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND phone: x-41-1-632-3408 fax: ...-1228 <><
>>>>> "MM" == Martin Maechler <maechler@stat.math.ethz.ch> >>>>> on Mon, 19 Jan 2004 19:03:29 +0100 (CET) writes:>>>>> "Remy" == Remy Drouilhet <Remy.Drouilhet@upmf-grenoble.fr> >>>>> on Mon, 19 Jan 2004 18:10:04 +0100 (CET) writes:Remy> qchisq(1,10) works well but qchisq(1,10,ncp=0) does Remy> not work whereas ncp=0 is the default value of the Remy> function qchisq(1,10). (of course, 10 will be replaced Remy> by any integer value). Remy> Let us notice that this bug occurs only when applying Remy> probability one. (qchisq(seq(0,.9,.1),10,ncp=0) works Remy> very well). MM> Merci, Remy, for the report. MM> Note that I've been working on improvements to "pnchisq" MM> (C) and "qnchisq" for R-devel for several days now. MM> You will see fixes there pretty soon, I hope. As a matter of fact, the bug you mention isn't new. That's been PR#6421 which has been fixed several weeks ago in both "R-patched" and "R-devel". (The other things I said above are true in any case.) Martin