[Rd] small bug on qchisq (PR#6442)

maechler at stat.math.ethz.ch maechler at stat.math.ethz.ch
Mon Jan 19 19:03:29 MET 2004


>>>>> "Remy" == Remy Drouilhet <Remy.Drouilhet at 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 at 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			<><



More information about the R-devel mailing list