'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: [EMAIL PROTECTED];LC_NUMERIC=C;[EMAIL PROTECTED];LC_COLLATE=C;[EMAIL PROTECTED];[EMAIL PROTECTED];[EMAIL PROTECTED];LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;[EMAIL PROTECTED];LC_IDENTIFICATION=C Search Path: .GlobalEnv, package:methods, package:stats, package:utils, package:datasets, package:graphics, package:grDevices, Autoloads, package:base ______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel