>>>>> "Adam" == Adam Witney <[EMAIL PROTECTED]>
>>>>> 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" <[EMAIL PROTECTED]> wrote:
>>>>>>> "UweL" == Uwe Ligges <[EMAIL PROTECTED]>
>>>>>>> on Thu, 24 May 2007 17:34:16 +0200 writes:
>>
UweL> Some of these test are expected from time to time, since they are
>> using
UweL> random numbers. Just re-run.
>>
>> eehm, "some of these", yes, but not the ones Adam mentioned,
>> d-p-q-r-tests.R.
>>
>> 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...
[1] 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
Execution halted
Ok, thanks;
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.
>>
>> Indeed!
>>
>> 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
Adam> adam
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
