>>>>> "Remy" == Remy Drouilhet <[EMAIL PROTECTED]>
>>>>> 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 <[EMAIL PROTECTED]> 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 <><
______________________________________________
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-devel