[R] Make check failure for R-2.4.1
maechler at stat.math.ethz.ch
Fri May 25 14:34:49 CEST 2007
>>>>> "Adam" == Adam Witney <awitney at sgul.ac.uk>
>>>>> on Fri, 25 May 2007 09:38:29 +0100 writes:
Adam> Thanks for your replies.... Details inline below:
Adam> On 24/5/07 17:12, "Martin Maechler" <maechler at stat.math.ethz.ch> wrote:
>>>>>>> "UweL" == Uwe Ligges <ligges at statistik.uni-dortmund.de>
>>>>>>> on Thu, 24 May 2007 17:34:16 +0200 writes:
UweL> Some of these test are expected from time to time, since they are
UweL> random numbers. Just re-run.
>> eehm, "some of these", yes, but not the ones Adam mentioned,
>> Adam, if you want more info you should report to us the *end*
>> (last dozen of lines) of
>> your d-p-q-r-tests.Rout[.fail] file.
Adam> Ok, here they are...
 TRUE TRUE TRUE TRUE
> ##-- non central Chi^2 :
> xB <- c(2000,1e6,1e50,Inf)
> for(df in c(0.1, 1, 10))
+ for(ncp in c(0, 1, 10, 100)) stopifnot(pchisq(xB, df=df, ncp=ncp) ==1)
Error: pchisq(xB, df = df, ncp = ncp) == 1 is not all TRUE
so, if we want to learn more, we need
the output of something like
xB <- c(2000,1e6,1e50,Inf)
for(df in c(0.1, 1, 10))
for(ncp in c(0, 1, 10, 100))
print(pchisq(xB, df=df, ncp=ncp), digits == 15)
UweL> BTW: We do have R-2.5.0 these days.
>> And gcc 2.95.4 is also very old.
>> Maybe you've recovered an old compiler / math-library bug from
>> that antique compiler suite ?
Adam> Yes, maybe I should start think about upgrading this box!
yes, at least "start" ... ;-)
Adam> Thanks again
More information about the R-help