* using log directory ‘/data/gannet/ripley/R/packages/tests-LENGTH1/bimixt.Rcheck’ * using R Under development (unstable) (2022-04-26 r82260) * using platform: x86_64-pc-linux-gnu (64-bit) * using session charset: UTF-8 * using option ‘--no-stop-on-test-error’ * checking for file ‘bimixt/DESCRIPTION’ ... OK * checking extension type ... Package * this is package ‘bimixt’ version ‘1.0’ * checking package namespace information ... OK * checking package dependencies ... OK * checking if this is a source package ... OK * checking if there is a namespace ... OK * checking for executable files ... OK * checking for hidden files and directories ... OK * checking for portable file names ... OK * checking for sufficient/correct file permissions ... OK * checking whether package ‘bimixt’ can be installed ... [8s/19s] OK * checking package directory ... OK * checking DESCRIPTION meta-information ... OK * checking top-level files ... OK * checking for left-over files ... OK * checking index information ... OK * checking package subdirectories ... OK * checking R files for non-ASCII characters ... OK * checking R files for syntax errors ... OK * checking whether the package can be loaded ... OK * checking whether the package can be loaded with stated dependencies ... OK * checking whether the package can be unloaded cleanly ... OK * checking whether the namespace can be loaded with stated dependencies ... OK * checking whether the namespace can be unloaded cleanly ... OK * checking loading without being on the library search path ... OK * checking use of S3 registration ... OK * checking dependencies in R code ... OK * checking S3 generic/method consistency ... OK * checking replacement functions ... OK * checking foreign function calls ... OK * checking R code for possible problems ... [17s/47s] OK * checking Rd files ... OK * checking Rd metadata ... OK * checking Rd line widths ... OK * checking Rd cross-references ... OK * checking for missing documentation entries ... OK * checking for code/documentation mismatches ... OK * checking Rd \usage sections ... OK * checking Rd contents ... OK * checking for unstated dependencies in examples ... OK * checking examples ... [34s/86s] ERROR Running examples in ‘bimixt-Ex.R’ failed The error most likely occurred in: > ### Name: lr.test > ### Title: lr.test > ### Aliases: lr.test > > ### ** Examples > > case=rmix(50,10,1.2,15,1,.7) > control=rmix(50,10,1.2,15,1,.95) > model1=bimixt.model(case=case,control=control, type="4c") > model2=bimixt.model(case=case,control=control, type="binorm") > lr.test(model1, model2) ----------- FAILURE REPORT -------------- --- failure: length > 1 in coercion to logical --- --- srcref --- : --- package (from environment) --- bimixt --- call from context --- lr.test(model1, model2) --- call from argument --- model1$control != model2$control || model1$case != model2$case --- R stacktrace --- where 1: lr.test(model1, model2) --- value of length: 50 type: logical --- [1] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE [13] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE [25] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE [37] FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE FALSE [49] FALSE FALSE --- function from context --- function (model1, model2) { if (model1$control != model2$control || model1$case != model2$case) stop("The data are not consistent across models. LR test is not applicable.") if (model1$type == model2$type) warning("models are equivalent") my.chi <- function(ll0, lla, df) { test.stat = 2 * abs(lla - ll0) p.val = 1 - pchisq(test.stat, df) return(p.val) } if (model1$type == "binorm") { if (model2$type == "binorm") { p = my.chi(model1$max.loglike, model2$max.loglike, 0) } if (model2$type == "2cc") { if (model1$max.loglike > model2$max.loglike) { p = NA } else { p = my.chi(model1$max.loglike, model2$max.loglike, 1) } } if (model2$type == "2cu") { if (model1$max.loglike > model2$max.loglike) { p = NA } else { p = my.chi(model1$max.loglike, model2$max.loglike, 2) } } if (model2$type == "4c") { if (model1$max.loglike > model2$max.loglike) { p = NA } else { p = my.chi(model1$max.loglike, model2$max.loglike, 6) } } } if (model1$type == "2cc") { if (model2$type == "binorm") { if (model1$max.loglike < model2$max.loglike) { p = NA } else { p = my.chi(model2$max.loglike, model1$max.loglike, 1) } } if (model2$type == "2cc") { p = my.chi(model1$max.loglike, model2$max.loglike, 0) } if (model2$type == "2cu") { if (model1$max.loglike > model2$max.loglike) { p = NA } else { p = my.chi(model1$max.loglike, model2$max.loglike, 1) } } if (model2$type == "4c") { if (model1$max.loglike > model2$max.loglike) { p = NA } else { p = my.chi(model1$max.loglike, model2$max.loglike, 5) } } } if (model1$type == "2cu") { if (model2$type == "binorm") { if (model1$max.loglike < model2$max.loglike) { p = NA } else { p = my.chi(model2$max.loglike, model1$max.loglike, 2) } } if (model2$type == "2cc") { if (model1$max.loglike < model2$max.loglike) { p = NA } else { p = my.chi(model2$max.loglike, model1$max.loglike, 1) } } if (model2$type == "2cu") { p = my.chi(model2$max.loglike, model1$max.loglike, 0) } if (model2$type == "4c") { if (model1$max.loglike > model2$max.loglike) { p = NA } else { p = my.chi(model1$max.loglike, model2$max.loglike, 4) } } } if (model1$type == "4c") { if (model2$type == "binorm") { if (model1$max.loglike < model2$max.loglike) { p = NA } p = my.chi(model2$max.loglike, model1$max.loglike, 6) } if (model2$type == "2cc") { if (model1$max.loglike < model2$max.loglike) { p = NA } else { p = my.chi(model2$max.loglike, model1$max.loglike, 5) } } if (model2$type == "2cu") { if (model1$max.loglike < model2$max.loglike) { p = NA } else { p = my.chi(model2$max.loglike, model1$max.loglike, 4) } } if (model2$type == "4c") { p = my.chi(model1$max.loglike, model2$max.loglike, 0) } } if (is.na(p)) { warning("The smaller model has a larger maximum likelihood value than the larger model.") } names(p) = "p value" return(p) } --- function search by body --- Function lr.test in namespace bimixt has this body. ----------- END OF FAILURE REPORT -------------- Fatal error: length > 1 in coercion to logical * checking PDF version of manual ... OK * checking for non-standard things in the check directory ... OK * checking for detritus in the temp directory ... OK * DONE Status: 1 ERROR See ‘/data/gannet/ripley/R/packages/tests-LENGTH1/bimixt.Rcheck/00check.log’ for details. Command exited with non-zero status 1 Time 4:01.24, 83.91 + 7.05