From p.c.j.kleiweg at rug.nl  Tue Oct 31 11:47:23 2006
From: p.c.j.kleiweg at rug.nl (p.c.j.kleiweg at rug.nl)
Date: Tue, 31 Oct 2006 11:47:23 +0100 (CET)
Subject: [Rd] 'make check' fails on d-p-q-r-tests (PR#9326)
Message-ID: <20061031104723.0E81FAC02@slim.kubism.ku.dk>


'make check' fails on d-p-q-r-tests:

> ##-- 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


Here is some more testing:

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)); print(pchisq(xB, df=df, ncp=ncp) == 1) }
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] FALSE FALSE FALSE  TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] FALSE FALSE FALSE  TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] FALSE FALSE FALSE  TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE

options(digits=22)
for(df in c(0.1, 1, 10)) for(ncp in c(0, 1, 10, 100)) { print(pchisq(xB, df=df, ncp=ncp)); print(pchisq(xB, df=df, ncp=ncp) == 1) }
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 0.9999999999999998889777 0.9999999999999998889777 0.9999999999999998889777 1.0000000000000000000000
[1] FALSE FALSE FALSE  TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 0.9999999999999998889777 0.9999999999999998889777 0.9999999999999998889777 1.0000000000000000000000
[1] FALSE FALSE FALSE  TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE
[1] 0.9999999999999998889777 0.9999999999999998889777 0.9999999999999998889777 1.0000000000000000000000
[1] FALSE FALSE FALSE  TRUE
[1] 1 1 1 1
[1] TRUE TRUE TRUE TRUE



--please do not edit the information below--

Version:
 platform = i686-pc-linux-gnu
 arch = i686
 os = linux-gnu
 system = i686, linux-gnu
 status = 
 major = 2
 minor = 4.0
 year = 2006
 month = 10
 day = 03
 svn rev = 39566
 language = R
 version.string = R version 2.4.0 (2006-10-03)

Locale:
LC_CTYPE=nl_NL at euro;LC_NUMERIC=C;LC_TIME=nl_NL at euro;LC_COLLATE=C;LC_MONETARY=nl_NL at euro;LC_MESSAGES=nl_NL at euro;LC_PAPER=nl_NL at euro;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=nl_NL at euro;LC_IDENTIFICATION=C

Search Path:
 .GlobalEnv, package:methods, package:stats, package:utils, package:datasets, package:graphics, package:grDevices, Autoloads, package:base




[Rd] 'make check' fails on d-p-q-r-tests (PR#9326)
p.c.j.kleiweg at rug.nl