From maechler at stat.math.ethz.ch  Fri Jan  9 19:04:16 2004
From: maechler at stat.math.ethz.ch (maechler@stat.math.ethz.ch)
Date: Fri Jan  9 19:03:31 2004
Subject: [Rd] qchisq hang (PR#6421)
Message-ID: <20040109180416.4611E10497@slim.kubism.ku.dk>

>>>>> "MM" == Martin Maechler <maechler@stat.math.ethz.ch>
>>>>>     on Thu,  8 Jan 2004 14:09:29 +0100 (CET) writes:

>>>>> "ThomasSt" == statho3  <statho3@web.de>
>>>>>     on Wed,  7 Jan 2004 21:45:40 +0100 (CET) writes:

    ThomasSt> Full_Name: Thomas Stabla Version: 1.8.1 OS: Linux,
    ThomasSt> XP Submission from: (NULL) (80.136.67.200)


    ThomasSt> qchisq(1, df = 1, ncp = 0) seems to hang for me.

    MM> problem confirmed; already for qchisq(1 - 1e-7, df = 1,
    MM> ncp = 0)

    MM> Looking at the source code, the problem is "obvious".

    MM> The analogon on the tail (that hangs) is

    MM> 	qchisq(1e-8, df = 1, ncp = 0, lower.tail = FALSE)

A patch for these has now been committed,
both for R-patched (1.8.2 never to be) 
and R-devel (1.9.0 to be).

Martin Maechler <maechler@stat.math.ethz.ch>	http://stat.ethz.ch/~maechler/


[Rd] qchisq hang (PR#6421)
maechler@stat.math.ethz.ch