##############################################################################
##############################################################################
###
### Running command:
###
###   /Library/Frameworks/R.framework/Resources/bin/R CMD check --install=check:POMA.install-out.txt --library=/Library/Frameworks/R.framework/Resources/library --no-vignettes --timings POMA_1.13.26.tar.gz
###
##############################################################################
##############################################################################


* using log directory ‘/Users/biocbuild/bbs-3.19-bioc/meat/POMA.Rcheck’
* using R Under development (unstable) (2024-03-18 r86148)
* using platform: x86_64-apple-darwin20
* R was compiled by
    Apple clang version 14.0.0 (clang-1400.0.29.202)
    GNU Fortran (GCC) 12.2.0
* running under: macOS Monterey 12.7.1
* using session charset: UTF-8
* using option ‘--no-vignettes’
* checking for file ‘POMA/DESCRIPTION’ ... OK
* this is package ‘POMA’ version ‘1.13.26’
* package encoding: UTF-8
* 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 hidden files and directories ... OK
* checking for portable file names ... OK
* checking for sufficient/correct file permissions ... OK
* checking whether package ‘POMA’ can be installed ... OK
* checking installed package size ... OK
* checking package directory ... OK
* checking ‘build’ 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 code 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 whether startup messages can be suppressed ... OK
* checking dependencies in R code ... NOTE
Unexported object imported by a ':::' call: ‘utils:::.getHelpFile’
  See the note in ?`:::` about the use of this operator.
* checking S3 generic/method consistency ... OK
* checking replacement functions ... OK
* checking foreign function calls ... OK
* checking R code for possible problems ... NOTE
PomaBatch: no visible global function definition for ‘is’
PomaBatch: no visible global function definition for ‘validObject’
PomaBoxplots: no visible global function definition for ‘is’
PomaBoxplots: no visible binding for global variable ‘sample_id’
PomaBoxplots: no visible binding for global variable ‘group_factor’
PomaBoxplots: no visible global function definition for ‘median’
PomaBoxplots: no visible binding for global variable ‘value’
PomaBoxplots: no visible global function definition for ‘reorder’
PomaBoxplots: no visible binding for global variable ‘median_rank’
PomaBoxplots: no visible binding for global variable ‘name’
PomaClust: no visible global function definition for ‘is’
PomaClust: no visible binding for global variable ‘clust’
PomaClust: no visible binding for global variable ‘V1’
PomaClust: no visible binding for global variable ‘V2’
PomaClust: no visible binding for global variable ‘Dim1’
PomaClust: no visible binding for global variable ‘Dim2’
PomaCorr: no visible global function definition for ‘is’
PomaCorr: no visible global function definition for ‘cor’
PomaCorr: no visible binding for global variable ‘column’
PomaCorr: no visible binding for global variable ‘cor’
PomaCorr: no visible binding for global variable ‘p’
PomaCorr: no visible global function definition for ‘p.adjust’
PomaCorr: no visible binding for global variable ‘pvalue’
PomaCorr: no visible binding for global variable ‘corr’
PomaCorr: no visible binding for global variable ‘feature1’
PomaCorr: no visible binding for global variable ‘name’
PomaCorr: no visible binding for global variable ‘value’
PomaCreateObject: no visible global function definition for
  ‘validObject’
PomaDESeq: no visible global function definition for ‘is’
PomaDESeq: no visible binding for global variable ‘padj’
PomaDESeq: no visible binding for global variable ‘adj_pvalue’
PomaDensity: no visible global function definition for ‘is’
PomaDensity: no visible binding for global variable ‘sample_id’
PomaDensity: no visible binding for global variable ‘group_factor’
PomaDensity: no visible binding for global variable ‘name’
PomaDensity: no visible binding for global variable ‘value’
PomaHeatmap: no visible global function definition for ‘is’
PomaImpute: no visible global function definition for ‘is’
PomaImpute: no visible global function definition for ‘aggregate’
PomaImpute: no visible binding for global variable ‘group_factor’
PomaImpute: no visible global function definition for ‘validObject’
PomaLM: no visible global function definition for ‘is’
PomaLM: no visible global function definition for ‘lm’
PomaLM: no visible global function definition for ‘p.adjust’
PomaLM: no visible binding for global variable ‘p.value’
PomaLM: no visible binding for global variable ‘term’
PomaLM: no visible binding for global variable ‘estimate’
PomaLM: no visible binding for global variable ‘std.error’
PomaLM: no visible binding for global variable ‘statistic’
PomaLM: no visible binding for global variable ‘adj_pvalue’
PomaLM: no visible binding for global variable ‘pvalue’
PomaLM: no visible global function definition for ‘reorder’
PomaLM: no visible binding for global variable ‘feature’
PomaLM: no visible binding for global variable ‘std_err’
PomaLMM: no visible global function definition for ‘is’
PomaLMM : lmm_fun: no visible global function definition for ‘vcov’
PomaLMM : lmm_fun: no visible binding for global variable ‘vcov’
PomaLMM : lmm_fun: no visible binding for global variable ‘grp’
PomaLMM : lmm_fun: no visible binding for global variable
  ‘variance_percent’
PomaLMM: no visible binding for global variable ‘feature’
PomaLMM: no visible binding for global variable ‘name’
PomaLMM: no visible binding for global variable ‘value’
PomaLasso: no visible global function definition for ‘is’
PomaLasso: no visible binding for global variable ‘estimate’
PomaLasso: no visible binding for global variable ‘conf.low’
PomaLasso: no visible binding for global variable ‘conf.high’
PomaLasso: no visible global function definition for ‘predict’
PomaLasso: no visible binding for global variable ‘term’
PomaLimma: no visible global function definition for ‘is’
PomaLimma: no visible global function definition for ‘as.formula’
PomaLimma: no visible binding for global variable ‘P.Value’
PomaLimma: no visible binding for global variable ‘adj.P.Val’
PomaNorm: no visible global function definition for ‘is’
PomaNorm: no visible binding for global variable ‘var’
PomaNorm : <anonymous>: no visible global function definition for ‘sd’
PomaNorm: no visible global function definition for ‘validObject’
PomaOddsRatio: no visible global function definition for ‘is’
PomaOddsRatio: no visible binding for global variable ‘group’
PomaOddsRatio: no visible binding for global variable ‘OddsRatio’
PomaOddsRatio: no visible binding for global variable ‘feature’
PomaOddsRatio: no visible binding for global variable ‘upr’
PomaOddsRatio: no visible binding for global variable ‘lwr’
PomaOutliers: no visible global function definition for ‘is’
PomaOutliers: no visible binding for global variable ‘groups’
PomaOutliers: no visible global function definition for ‘quantile’
PomaOutliers: no visible global function definition for ‘IQR’
PomaOutliers: no visible binding for global variable ‘limit’
PomaOutliers: no visible binding for global variable ‘out’
PomaOutliers : find_hull: no visible global function definition for
  ‘chull’
PomaOutliers: no visible binding for global variable ‘group’
PomaOutliers: no visible binding for global variable ‘.’
PomaOutliers: no visible binding for global variable ‘PCoA1’
PomaOutliers: no visible binding for global variable ‘PCoA2’
PomaOutliers: no visible global function definition for ‘validObject’
PomaPCA: no visible global function definition for ‘is’
PomaPCA: no visible global function definition for ‘prcomp’
PomaPCA: no visible binding for global variable ‘PC1’
PomaPCA: no visible binding for global variable ‘PC2’
PomaPCA: no visible binding for global variable ‘group’
PomaPCA: no visible binding for global variable ‘sample_id’
PomaPCA: no visible global function definition for ‘reorder’
PomaPCA: no visible binding for global variable ‘comp’
PomaPCA: no visible binding for global variable ‘var_exp’
PomaPCA: no visible binding for global variable ‘feature’
PomaPCA: no visible binding for global variable ‘value’
PomaPCA: no visible binding for global variable ‘name’
PomaPCA: no visible binding for global variable ‘to_x’
PomaPCA: no visible binding for global variable ‘to_y’
PomaPCR: no visible global function definition for ‘is’
PomaPCR: no visible binding for global variable ‘PC1’
PomaPCR: no visible global function definition for ‘lm’
PomaPCR: no visible global function definition for ‘p.adjust’
PomaPCR: no visible binding for global variable ‘p.value’
PomaPCR: no visible binding for global variable ‘term’
PomaPCR: no visible binding for global variable ‘estimate’
PomaPCR: no visible binding for global variable ‘std.error’
PomaPCR: no visible binding for global variable ‘statistic’
PomaPCR: no visible binding for global variable ‘adj_pvalue’
PomaPCR: no visible binding for global variable ‘pvalue’
PomaPLS: no visible global function definition for ‘is’
PomaPLS: no visible binding for global variable ‘comp1’
PomaPLS: no visible binding for global variable ‘comp2’
PomaPLS: no visible binding for global variable ‘sample_id’
PomaPLS: no visible binding for global variable ‘feature’
PomaPLS: no visible global function definition for ‘reorder’
PomaPLS: no visible binding for global variable ‘value’
PomaPLS: no visible binding for global variable ‘name’
PomaPLS: no visible binding for global variable ‘component’
PomaPLS: no visible binding for global variable ‘error’
PomaPLS: no visible binding for global variable ‘feature_sd’
PomaPLS: no visible binding for global variable ‘sd’
PomaRandForest: no visible global function definition for ‘is’
PomaRandForest: no visible binding for global variable ‘OOB’
PomaRandForest: no visible binding for global variable
  ‘MeanDecreaseGini’
PomaRandForest: no visible global function definition for ‘reorder’
PomaRandForest: no visible binding for global variable ‘feature’
PomaRankProd: no visible global function definition for ‘is’
PomaRankProd: no visible binding for global variable ‘P.value’
PomaRankProd: no visible binding for global variable ‘gene.index’
PomaUMAP: no visible global function definition for ‘is’
PomaUMAP: no visible binding for global variable ‘clust’
PomaUMAP: no visible binding for global variable ‘UMAP1’
PomaUMAP: no visible binding for global variable ‘UMAP2’
PomaUnivariate: no visible global function definition for ‘is’
PomaUnivariate: no visible binding for global variable ‘group’
PomaUnivariate : <anonymous>: no visible global function definition for
  ‘t.test’
PomaUnivariate: no visible global function definition for ‘p.adjust’
PomaUnivariate: no visible binding for global variable ‘pvalue’
PomaUnivariate: no visible binding for global variable ‘feature’
PomaUnivariate: no visible binding for global variable ‘fold_change’
PomaUnivariate: no visible binding for global variable ‘diff_means’
PomaUnivariate: no visible binding for global variable ‘adj_pvalue’
PomaUnivariate : <anonymous>: no visible global function definition for
  ‘anova’
PomaUnivariate : <anonymous>: no visible global function definition for
  ‘aov’
PomaUnivariate: no visible global function definition for ‘TukeyHSD’
PomaUnivariate: no visible global function definition for ‘aov’
PomaUnivariate: no visible binding for global variable ‘adj.p.value’
PomaUnivariate: no visible global function definition for ‘as.formula’
PomaUnivariate: no visible binding for global variable ‘term’
PomaUnivariate: no visible binding for global variable ‘p.value’
PomaUnivariate: no visible binding for global variable ‘contrast’
PomaUnivariate: no visible binding for global variable ‘p adj’
PomaUnivariate : <anonymous>: no visible global function definition for
  ‘wilcox.test’
PomaUnivariate : <anonymous>: no visible global function definition for
  ‘kruskal.test’
PomaUnivariate: no visible binding for global variable ‘kw_rank_sum’
PomaUnivariate: no visible binding for global variable ‘Comparison’
PomaUnivariate: no visible binding for global variable ‘P.adj’
PomaVolcano: no visible global function definition for ‘quantile’
PomaVolcano: no visible binding for global variable ‘logFC’
PomaVolcano: no visible binding for global variable ‘pvalue’
PomaVolcano: no visible binding for global variable ‘feature’
create_mock_data: no visible global function definition for ‘runif’
create_mock_summarized_experiment: no visible global function
  definition for ‘runif’
help_extract: no visible global function definition for
  ‘capture.output’
make_legend: no visible binding for global variable ‘POMA’
title_extract: no visible global function definition for
  ‘capture.output’
Undefined global functions or variables:
  . Comparison Dim1 Dim2 IQR MeanDecreaseGini OOB OddsRatio P.Value
  P.adj P.value PC1 PC2 PCoA1 PCoA2 POMA TukeyHSD UMAP1 UMAP2 V1 V2
  adj.P.Val adj.p.value adj_pvalue aggregate anova aov as.formula
  capture.output chull clust column comp comp1 comp2 component
  conf.high conf.low contrast cor corr diff_means error estimate
  feature feature1 feature_sd fold_change gene.index group group_factor
  groups grp is kruskal.test kw_rank_sum limit lm logFC lwr median
  median_rank name out p p adj p.adjust p.value padj prcomp predict
  pvalue quantile reorder runif sample_id sd statistic std.error
  std_err t.test term to_x to_y upr validObject value var var_exp
  variance_percent vcov wilcox.test
Consider adding
  importFrom("grDevices", "chull")
  importFrom("methods", "is", "validObject")
  importFrom("stats", "IQR", "TukeyHSD", "aggregate", "anova", "aov",
             "as.formula", "cor", "kruskal.test", "lm", "median",
             "p.adjust", "prcomp", "predict", "quantile", "reorder",
             "runif", "sd", "t.test", "var", "vcov", "wilcox.test")
  importFrom("utils", "capture.output")
to your NAMESPACE file (and ensure that your DESCRIPTION Imports field
contains 'methods').
* checking Rd files ... OK
* checking Rd metadata ... 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 contents of ‘data’ directory ... OK
* checking data for non-ASCII characters ... OK
* checking LazyData ... OK
* checking data for ASCII and uncompressed saves ... OK
* checking files in ‘vignettes’ ... OK
* checking examples ...sh: line 1: 20604 Segmentation fault: 11  LANGUAGE=en _R_CHECK_INTERNALS2_=1 R_LIBS=/tmp/Rtmp0CQ1JI/RLIBS_45af3fc9f12b R_ENVIRON_USER='' R_LIBS_USER='NULL' R_LIBS_SITE='NULL' '/Library/Frameworks/R.framework/Resources/bin/R' --vanilla > 'POMA-Ex.Rout' 2>&1 < 'POMA-Ex.R'
 ERROR
Running examples in ‘POMA-Ex.R’ failed
The error most likely occurred in:

> base::assign(".ptime", proc.time(), pos = "CheckExEnv")
> ### Name: PomaLMM
> ### Title: Linear Mixed Models
> ### Aliases: PomaLMM
> 
> ### ** Examples
> 
> data("st000284")
> 
> # Perform linear mixed model with all features
> st000284 %>% 
+ PomaLMM()

 *** caught segfault ***
address 0x0, cause 'unknown'

Traceback:
 1: .Call(merPredDCreate, as(X, "matrix"), Lambdat, LamtUt, Lind,     RZX, Ut, Utr, V, VtV, Vtr, Xwts, Zt, beta0, delb, delu, theta,     u0)
 2: initializePtr()
 3: .Object$initialize(...)
 4: initialize(value, ...)
 5: initialize(value, ...)
 6: methods::new(def, ...)
 7: (new("refMethodDef", .Data = function (...) {    methods::new(def, ...)}, mayCall = c("methods", "new"), name = "new", refClassName = "refGeneratorSlot",     superClassMethod = ""))(Zt = new("dgCMatrix", i = c(2L, 3L, 6L, 10L, 1L, 4L, 6L, 10L, 2L, 3L, 6L, 9L, 2L, 3L, 6L, 10L, 1L, 4L, 6L, 9L, 1L, 4L, 6L, 10L, 0L, 4L, 6L, 9L, 2L, 3L, 6L, 10L, 1L, 4L, 6L, 9L, 1L, 4L, 6L, 9L, 2L, 4L, 6L, 9L, 1L, 4L, 6L, 9L, 2L, 5L, 6L, 10L, 1L, 4L, 6L, 10L, 2L, 4L, 6L, 9L, 2L, 3L, 6L, 10L, 1L, 4L, 6L, 10L, 1L, 4L, 6L, 9L, 2L, 5L, 6L, 10L, 2L, 3L, 6L, 10L, 2L, 4L, 6L, 9L, 1L, 4L, 6L, 10L, 2L, 3L, 6L, 9L, 1L, 4L, 6L, 10L, 2L, 3L, 6L, 9L, 1L, 4L, 6L, 10L, 2L, 4L, 6L, 9L, 1L, 4L, 6L, 9L, 1L, 4L, 6L, 10L, 2L, 4L, 6L, 10L, 2L, 3L, 6L, 9L, 1L, 4L, 6L, 9L, 2L, 4L, 6L, 10L, 1L, 4L, 6L, 9L, 2L, 3L, 6L, 10L, 1L, 4L, 6L, 9L, 2L, 5L, 6L, 10L, 2L, 3L, 6L, 10L, 2L, 3L, 6L, 9L, 2L, 3L, 6L, 10L, 2L, 3L, 6L, 9L, 2L, 3L, 6L, 10L, 0L, 4L, 6L, 9L, 2L, 4L, 6L, 9L, 2L, 3L, 6L, 9L, 1L, 4L, 6L, 9L, 2L, 4L, 6L, 10L, 0L, 4L, 6L, 10L, 1L, 4L, 6L, 9L, 1L, 4L, 6L, 10L, 1L, 4L, 6L, 10L, 1L, 4L, 6L, 9L, 2L, 4L, 6L, 9L, 1L, 4L, 6L, 10L, 2L, 5L, 6L, 10L, 2L, 3L, 6L, 9L, 1L, 4L, 6L, 9L, 1L, 4L, 6L, 10L, 1L, 4L, 6L, 9L, 1L, 4L, 6L, 10L, 1L, 4L, 6L, 10L, 2L, 3L, 6L, 10L, 2L, 4L, 6L, 10L, 2L, 4L, 6L, 10L, 2L, 3L, 7L, 9L, 0L, 4L, 7L, 9L, 1L, 4L, 7L, 9L, 1L, 4L, 7L, 9L, 1L, 4L, 7L, 9L, 0L, 3L, 7L, 10L, 1L, 3L, 7L, 10L, 1L, 4L, 7L, 10L, 2L, 5L, 7L, 9L, 1L, 3L, 7L, 9L, 1L, 3L, 7L, 10L, 1L, 5L, 7L, 9L, 2L, 3L, 7L, 10L, 2L, 4L, 7L, 9L, 1L, 3L, 7L, 10L, 1L, 3L, 7L, 10L, 2L, 4L, 7L, 9L, 2L, 3L, 7L, 9L, 2L, 3L, 7L, 9L, 2L, 4L, 7L, 10L, 2L, 3L, 7L, 10L, 1L, 4L, 7L, 10L, 1L, 4L, 7L, 9L, 2L, 3L, 7L, 10L, 1L, 4L, 7L, 9L, 2L, 3L, 7L, 10L, 2L, 4L, 7L, 9L, 0L, 4L, 7L, 9L, 2L, 4L, 7L, 10L, 1L, 3L, 7L, 9L, 2L, 3L, 7L, 9L, 2L, 3L, 7L, 9L, 1L, 3L, 7L, 10L, 0L, 4L, 7L, 9L, 2L, 4L, 7L, 9L, 2L, 3L, 7L, 9L, 2L, 4L, 7L, 9L, 0L, 4L, 7L, 9L, 1L, 3L, 7L, 10L, 1L, 4L, 7L, 9L, 1L, 3L, 7L, 10L, 1L, 4L, 7L, 10L, 2L, 4L, 7L, 10L, 1L, 4L, 7L, 10L, 2L, 4L, 7L, 9L, 1L, 4L, 7L, 10L, 2L, 4L, 7L, 9L, 2L, 4L, 7L, 10L, 1L, 4L, 7L, 10L, 2L, 3L, 7L, 9L, 2L, 4L, 7L, 10L, 1L, 5L, 7L, 9L, 2L, 3L, 7L, 10L, 2L, 3L, 7L, 10L, 0L, 4L, 7L, 10L, 0L, 4L, 7L, 10L, 1L, 3L, 7L, 10L, 2L, 4L, 7L, 9L, 2L, 4L, 7L, 10L, 1L, 4L, 7L, 10L, 1L, 3L, 7L, 10L, 1L, 3L, 7L, 10L, 0L, 4L, 7L, 9L, 1L, 5L, 7L, 9L, 1L, 4L, 7L, 10L, 2L, 3L, 7L, 9L, 2L, 5L, 7L, 9L, 2L, 3L, 7L, 10L, 1L, 5L, 7L, 9L, 1L, 4L, 7L, 9L, 2L, 4L, 7L, 9L, 1L, 4L, 7L, 9L, 0L, 4L, 7L, 10L, 1L, 4L, 7L, 9L, 1L, 4L, 7L, 9L, 1L, 5L, 7L, 9L, 2L, 3L, 7L, 9L, 2L, 4L, 7L, 10L, 1L, 3L, 7L, 9L, 2L, 3L, 7L, 10L, 2L, 5L, 7L, 10L, 1L, 4L, 7L, 9L, 2L, 3L, 7L, 10L, 2L, 4L, 7L, 10L, 2L, 4L, 8L, 10L, 2L, 3L, 8L, 9L, 1L, 4L, 8L, 9L, 0L, 3L, 8L, 10L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 10L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 2L, 4L, 8L, 9L, 2L, 3L, 8L, 10L, 1L, 4L, 8L, 10L, 2L, 3L, 8L, 10L, 2L, 3L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 10L, 2L, 3L, 8L, 10L, 2L, 4L, 8L, 10L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 2L, 3L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 10L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 3L, 8L, 10L, 1L, 4L, 8L, 10L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 2L, 4L, 8L, 9L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 3L, 8L, 10L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 10L, 2L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 2L, 4L, 8L, 9L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 5L, 8L, 10L, 2L, 3L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 10L, 2L, 3L, 8L, 10L, 2L, 4L, 8L, 10L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 3L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 9L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 2L, 4L, 8L, 10L, 2L, 3L, 8L, 10L, 0L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 2L, 3L, 8L, 10L, 1L, 4L, 8L, 9L), p = c(0L, 4L, 8L, 12L, 16L, 20L, 24L, 28L, 32L, 36L, 40L, 44L, 48L, 52L, 56L, 60L, 64L, 68L, 72L, 76L, 80L, 84L, 88L, 92L, 96L, 100L, 104L, 108L, 112L, 116L, 120L, 124L, 128L, 132L, 136L, 140L, 144L, 148L, 152L, 156L, 160L, 164L, 168L, 172L, 176L, 180L, 184L, 188L, 192L, 196L, 200L, 204L, 208L, 212L, 216L, 220L, 224L, 228L, 232L, 236L, 240L, 244L, 248L, 252L, 256L, 260L, 264L, 268L, 272L, 276L, 280L, 284L, 288L, 292L, 296L, 300L, 304L, 308L, 312L, 316L, 320L, 324L, 328L, 332L, 336L, 340L, 344L, 348L, 352L, 356L, 360L, 364L, 368L, 372L, 376L, 380L, 384L, 388L, 392L, 396L, 400L, 404L, 408L, 412L, 416L, 420L, 424L, 428L, 432L, 436L, 440L, 444L, 448L, 452L, 456L, 460L, 464L, 468L, 472L, 476L, 480L, 484L, 488L, 492L, 496L, 500L, 504L, 508L, 512L, 516L, 520L, 524L, 528L, 532L, 536L, 540L, 544L, 548L, 552L, 556L, 560L, 564L, 568L, 572L, 576L, 580L, 584L, 588L, 592L, 596L, 600L, 604L, 608L, 612L, 616L, 620L, 624L, 628L, 632L, 636L, 640L, 644L, 648L, 652L, 656L, 660L, 664L, 668L, 672L, 676L, 680L, 684L, 688L, 692L, 696L, 700L, 704L, 708L, 712L, 716L, 720L, 724L, 728L, 732L, 736L, 740L, 744L, 748L, 752L, 756L, 760L, 764L, 768L, 772L, 776L, 780L, 784L, 788L, 792L, 796L, 800L, 804L, 808L, 812L, 816L, 820L, 824L, 828L, 832L, 836L, 840L, 844L, 848L, 852L, 856L, 860L, 864L, 868L, 872L, 876L, 880L, 884L, 888L, 892L, 896L), Dim = c(11L, 224L), Dimnames = list(c("0", "1", "2", "0", "1", "2", "CRC", "Healthy", "Polyp", "0", "1"),     c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11",     "12", "13", "14", "15", "16", "17", "18", "19", "20", "21",     "22", "23", "24", "25", "26", "27", "28", "29", "30", "31",     "32", "33", "34", "35", "36", "37", "38", "39", "40", "41",     "42", "43", "44", "45", "46", "47", "48", "49", "50", "51",     "52", "53", "54", "55", "56", "57", "58", "59", "60", "61",     "62", "63", "64", "65", "66", "67", "68", "69", "70", "71",     "72", "73", "74", "75", "76", "77", "78", "79", "80", "81",     "82", "83", "84", "85", "86", "87", "88", "89", "90", "91",     "92", "93", "94", "95", "96", "97", "98", "99", "100", "101",     "102", "103", "104", "105", "106", "107", "108", "109", "110",     "111", "112", "113", "114", "115", "116", "117", "118", "119",     "120", "121", "122", "123", "124", "125", "126", "127", "128",     "129", "130", "131", "132", "133", "134", "135", "136", "137",     "138", "139", "140", "141", "142", "143", "144", "145", "146",     "147", "148", "149", "150", "151", "152", "153", "154", "155",     "156", "157", "158", "159", "160", "161", "162", "163", "164",     "165", "166", "167", "168", "169", "170", "171", "172", "173",     "174", "175", "176", "177", "178", "179", "180", "181", "182",     "183", "184", "185", "186", "187", "188", "189", "190", "191",     "192", "193", "194", "195", "196", "197", "198", "199", "200",     "201", "202", "203", "204", "205", "206", "207", "208", "209",     "210", "211", "212", "213", "214", "215", "216", "217", "218",     "219", "220", "221", "222", "223", "224")), x = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), factors = list()), theta = c(1, 1, 1, 1), Lambdat = new("dgCMatrix", i = 0:10, p = 0:11, Dim = c(11L, 11L), Dimnames = list(NULL, NULL), x = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), factors = list()), Lind = c(1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L), n = 224L, X = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 45, 56, 53, 55, 55, 29, 51, 76, 66, 46, 52, 86, 65, 67, 66, 53, 71, 63, 37, 73, 68, 27, 47, 68, 74, 40, 48, 50, 66, 51, 42, 45, 50, 45, 60, 59, 58, 62, 86, 75, 30, 50, 60, 67, 88, 45, 70, 62, 68, 51, 63, 61, 42, 49, 62, 78, 50, 77, 56, 32, 70, 46, 67, 56, 59, 31, 44, 51, 58, 60, 60, 62, 48, 57, 50, 57, 63, 52, 44, 52, 56, 64, 63, 54, 57, 65, 49, 61, 64, 45, 44, 52, 49, 52, 70, 63, 38, 59, 62, 63, 41, 56, 57, 71, 68, 36, 43, 54, 18, 36, 25, 48, 75, 69, 54, 57, 64, 70, 54, 73, 80, 66, 50, 60, 60, 78, 74, 76, 71, 66, 44, 59, 56, 60, 21, 21, 24, 57, 47, 60, 64, 52, 53, 61, 49, 46, 53, 51, 56, 63, 58, 56, 51, 51, 53, 56, 63, 44, 46, 61, 61, 54, 59, 56, 57, 55, 54, 46, 63, 57, 64, 42, 86, 50, 63, 52, 62, 51, 55, 59, 52, 54, 52, 53, 60, 53, 48, 58, 50, 55, 50, 52, 57, 52, 54, 58, 58, 64, 62, 52, 55, 56, 58, 37, 46, 62, 61, 57, 49, 53, 63, 59, 51, 57, 41, 50, 61, 54, 62, 39, 64, 52, 63, 58))
 8: do.call(merPredD$new, c(reTrms[c("Zt", "theta", "Lambdat", "Lind")],     n = nrow(X), list(X = X)))
 9: (function (fr, X, reTrms, REML = TRUE, start = NULL, verbose = 0,     control = lmerControl(), ...) {    p <- ncol(X)    rho <- new.env(parent = parent.env(environment()))    rho$pp <- do.call(merPredD$new, c(reTrms[c("Zt", "theta",         "Lambdat", "Lind")], n = nrow(X), list(X = X)))    REMLpass <- if (REML)         p    else 0L    rho$resp <- if (missing(fr))         mkRespMod(REML = REMLpass, ...)    else mkRespMod(fr, REML = REMLpass)    pp <- resp <- NULL    rho$lmer_Deviance <- lmer_Deviance    devfun <- function(theta) .Call(lmer_Deviance, pp$ptr(),         resp$ptr(), as.double(theta))    environment(devfun) <- rho    if (is.null(start) && all(reTrms$cnms == "(Intercept)") &&         length(reTrms$flist) == length(reTrms$lower) && !is.null(y <- model.response(fr))) {        v <- sapply(reTrms$flist, function(f) var(ave(y, f)))        v.e <- var(y) - sum(v)        if (!is.na(v.e) && v.e > 0) {            v.rel <- v/v.e            if (all(v.rel >= reTrms$lower^2))                 rho$pp$setTheta(sqrt(v.rel))        }    }    if (length(rho$resp$y) > 0)         devfun(rho$pp$theta)    rho$lower <- reTrms$lower    devfun})(fr = list(response = c(390953.9725, 567984.4359, 558842.7545, 476949, 398317, 411224, 412250.6387, 414501.1114, 416166, 384549, 567128.8402, 458708, 368820, 363294.2595, 427550, 452920.4295, 572292.4009, 485945.981, 449289, 496170, 367490, 258236.7898, 357025.7389, 344961, 636235.025, 176265.8095, 585747.3231, 475436.0544, 496706.6782, 347051, 531463.2532, 406965, 433397.7677, 419417, 644673.6716, 295481.054, 424800.3286, 308847.5372, 399763.7731, 547712.3037, 364836.3348, 491016.6547, 356220.4178, 406035.9345, 537138, 441152, 464832, 407801.8058, 485609.0257, 411321.0878, 475664.8346, 489104.4686, 502646.3458, 380903.8611, 656531.0073, 357126.4125, 467118, 498204.9177, 578935.0375, 276282, 454961.8572, 354310.5058, 427676.4167, 335439.426, 477828.7284, 607468.3174, 525454.1885, 466164.2199, 388374.0707, 280601.4948, 538426.2521, 327845.8298, 310835.6962, 255581.1249, 310008.1058, 389023, 298260.2186, 253327.1314, 300536.256, 486057.4105, 327751.7548, 390423, 237308.4491, 1523682.102, 421348.114, 318941.1581, 310440.3791, 413432.0079, 361152.9157, 351938.1531, 343043.8915, 234086.0594, 431558.5135, 353378.2722, 662683, 619711, 393067, 436226, 357746.4897, 548544.9683, 446176.8347, 542946.7913, 370335, 395990.2168, 480289, 390791.7625, 405983.9203, 437346.4258, 411781.7759, 374353.0144, 348753.8564, 334510.9635, 455288.3659, 543543.8499, 475513.5542, 404896.4106, 577284.8669, 376892.3394, 612692.1242, 697588.0669, 608678.025, 480417.1678, 344548.6902, 377705.1731, 468790.538, 481168.2502, 349056.0956, 457475.515, 289511.447, 674400.4058, 442579.9069, 371751.64, 368035.0699, 364312.3182, 366940, 350912.5973, 407689.5782, 380417.6583, 367871.6673, 497267.9214, 350756.2318, 481736.4942, 359270.3287, 421373.045, 300508.5606, 241816.791, 353131.7319, 287105.2006, 597283.4463, 631177.7867, 474671.1668, 251142.6019, 334706.3165, 285859.1257, 302899.837, 679116.5808, 471889.8585, 294215.9893, 316314.4425, 430424.2632, 440289.1202, 375136.4106, 364348.4453, 350617.6839, 379959.6248, 263543.4117, 325594.1963, 254864.304, 374745.7934, 769700.8879, 265546.2757, 288690.5904, 448223.4843, 285781.3607, 478075.6061, 374421.358, 384377.8503, 425791.7744, 335740.0346, 372416.8626, 284134.6466, 242757.7857, 423237.2564, 449263.049, 369001.8556, 429189.2298, 488419.1339, 393047.9881, 436593.5514, 542624.9502, 518005.0948, 369574.1131, 395788.661, 411113.5836, 555784.8948, 615489.9839, 473038.7806, 350876.9488, 456311.1378, 438896.1988, 520021.259, 419016.127, 620559.6054, 353877.1194, 402095.5672, 335185.5304, 577251.1366, 484518.6279, 441758.2892, 349607.1891, 423897.7969, 410682.2447, 393578.3778, 429855.8436, 390493.7447, 375952.1076, 399195.9585, 1018447.846, 503210.6702, 311604.8089, 514713.0726, 407366.5548, 463101.8944, 510398.3294), age_at_consent = c(45, 56, 53, 55, 55, 29, 51, 76, 66, 46, 52, 86, 65, 67, 66, 53, 71, 63, 37, 73, 68, 27, 47, 68, 74, 40, 48, 50, 66, 51, 42, 45, 50, 45, 60, 59, 58, 62, 86, 75, 30, 50, 60, 67, 88, 45, 70, 62, 68, 51, 63, 61, 42, 49, 62, 78, 50, 77, 56, 32, 70, 46, 67, 56, 59, 31, 44, 51, 58, 60, 60, 62, 48, 57, 50, 57, 63, 52, 44, 52, 56, 64, 63, 54, 57, 65, 49, 61, 64, 45, 44, 52, 49, 52, 70, 63, 38, 59, 62, 63, 41, 56, 57, 71, 68, 36, 43, 54, 18, 36, 25, 48, 75, 69, 54, 57, 64, 70, 54, 73, 80, 66, 50, 60, 60, 78, 74, 76, 71, 66, 44, 59, 56, 60, 21, 21, 24, 57, 47, 60, 64, 52, 53, 61, 49, 46, 53, 51, 56, 63, 58, 56, 51, 51, 53, 56, 63, 44, 46, 61, 61, 54, 59, 56, 57, 55, 54, 46, 63, 57, 64, 42, 86, 50, 63, 52, 62, 51, 55, 59, 52, 54, 52, 53, 60, 53, 48, 58, 50, 55, 50, 52, 57, 52, 54, 58, 58, 64, 62, 52, 55, 56, 58, 37, 46, 62, 61, 57, 49, 53, 63, 59, 51, 57, 41, 50, 61, 54, 62, 39, 64, 52, 63, 58), factors = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L),     gender = c(2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L,     2L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L,     1L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L,     1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L,     2L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L,     1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L,     2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L,     2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L,     2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L,     1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L,     2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L,     2L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 1L,     1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L,     2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 1L), smoking_condition = c(1L, 2L, 1L, 1L, 2L, 2L, 2L,     1L, 2L, 2L, 2L, 2L, 3L, 2L, 2L, 1L, 2L, 2L, 3L, 1L, 2L, 2L,     1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 3L,     1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 2L, 3L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 2L,     2L, 2L, 1L, 1L, 2L, 3L, 1L, 1L, 3L, 1L, 2L, 1L, 1L, 2L, 1L,     1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,     2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 1L, 2L, 3L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 2L,     3L, 2L, 1L, 3L, 1L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 1L, 2L,     1L, 1L, 3L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L,     2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L,     1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 1L, 2L, 2L, 2L,     2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 2L, 1L, 2L, 2L, 1L, 2L), alcohol_consumption = c(3L,     2L, 3L, 3L, 2L, 2L, 1L, 3L, 2L, 2L, 3L, 2L, 3L, 2L, 3L, 3L,     2L, 2L, 3L, 3L, 3L, 2L, 3L, 2L, 3L, 2L, 3L, 2L, 2L, 3L, 3L,     2L, 3L, 2L, 3L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 3L, 3L, 2L,     3L, 1L, 2L, 2L, 2L, 2L, 3L, 2L, 3L, 3L, 2L, 2L, 2L, 2L, 2L,     3L, 3L, 3L, 3L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 3L, 2L, 2L, 2L,     3L, 3L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 3L, 2L, 3L, 3L,     1L, 3L, 2L, 3L, 3L, 2L, 1L, 3L, 3L, 3L, 1L, 2L, 2L, 2L, 2L,     3L, 2L, 3L, 2L, 3L, 3L, 2L, 3L, 3L, 2L, 3L, 3L, 1L, 1L, 2L,     3L, 3L, 2L, 2L, 2L, 1L, 2L, 2L, 3L, 3L, 3L, 2L, 2L, 3L, 2L,     1L, 2L, 2L, 2L, 3L, 3L, 2L, 3L, 3L, 2L, 3L, 3L, 3L, 3L, 2L,     1L, 3L, 2L, 3L, 2L, 2L, 2L, 2L, 3L, 3L, 2L, 3L, 3L, 2L, 3L,     3L, 3L, 2L, 2L, 3L, 2L, 3L, 2L, 2L, 3L, 3L, 2L, 2L, 3L, 2L,     2L, 3L, 2L, 2L, 3L, 3L, 2L, 2L, 3L, 2L, 3L, 3L, 2L, 3L, 3L,     2L, 3L, 3L, 2L, 2L, 2L, 3L, 2L, 2L, 3L, 3L, 3L, 2L, 3L, 2L,     2L, 3L, 3L, 2L, 2L, 2L, 2L, 3L, 3L, 1L, 2L, 3L, 2L)), X = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 45, 56, 53, 55, 55, 29, 51, 76, 66, 46, 52, 86, 65, 67, 66, 53, 71, 63, 37, 73, 68, 27, 47, 68, 74, 40, 48, 50, 66, 51, 42, 45, 50, 45, 60, 59, 58, 62, 86, 75, 30, 50, 60, 67, 88, 45, 70, 62, 68, 51, 63, 61, 42, 49, 62, 78, 50, 77, 56, 32, 70, 46, 67, 56, 59, 31, 44, 51, 58, 60, 60, 62, 48, 57, 50, 57, 63, 52, 44, 52, 56, 64, 63, 54, 57, 65, 49, 61, 64, 45, 44, 52, 49, 52, 70, 63, 38, 59, 62, 63, 41, 56, 57, 71, 68, 36, 43, 54, 18, 36, 25, 48, 75, 69, 54, 57, 64, 70, 54, 73, 80, 66, 50, 60, 60, 78, 74, 76, 71, 66, 44, 59, 56, 60, 21, 21, 24, 57, 47, 60, 64, 52, 53, 61, 49, 46, 53, 51, 56, 63, 58, 56, 51, 51, 53, 56, 63, 44, 46, 61, 61, 54, 59, 56, 57, 55, 54, 46, 63, 57, 64, 42, 86, 50, 63, 52, 62, 51, 55, 59, 52, 54, 52, 53, 60, 53, 48, 58, 50, 55, 50, 52, 57, 52, 54, 58, 58, 64, 62, 52, 55, 56, 58, 37, 46, 62, 61, 57, 49, 53, 63, 59, 51, 57, 41, 50, 61, 54, 62, 39, 64, 52, 63, 58), reTrms = list(Zt = new("dgCMatrix",     i = c(2L, 3L, 6L, 10L, 1L, 4L, 6L, 10L, 2L, 3L, 6L, 9L, 2L,     3L, 6L, 10L, 1L, 4L, 6L, 9L, 1L, 4L, 6L, 10L, 0L, 4L, 6L,     9L, 2L, 3L, 6L, 10L, 1L, 4L, 6L, 9L, 1L, 4L, 6L, 9L, 2L,     4L, 6L, 9L, 1L, 4L, 6L, 9L, 2L, 5L, 6L, 10L, 1L, 4L, 6L,     10L, 2L, 4L, 6L, 9L, 2L, 3L, 6L, 10L, 1L, 4L, 6L, 10L, 1L,     4L, 6L, 9L, 2L, 5L, 6L, 10L, 2L, 3L, 6L, 10L, 2L, 4L, 6L,     9L, 1L, 4L, 6L, 10L, 2L, 3L, 6L, 9L, 1L, 4L, 6L, 10L, 2L,     3L, 6L, 9L, 1L, 4L, 6L, 10L, 2L, 4L, 6L, 9L, 1L, 4L, 6L,     9L, 1L, 4L, 6L, 10L, 2L, 4L, 6L, 10L, 2L, 3L, 6L, 9L, 1L,     4L, 6L, 9L, 2L, 4L, 6L, 10L, 1L, 4L, 6L, 9L, 2L, 3L, 6L,     10L, 1L, 4L, 6L, 9L, 2L, 5L, 6L, 10L, 2L, 3L, 6L, 10L, 2L,     3L, 6L, 9L, 2L, 3L, 6L, 10L, 2L, 3L, 6L, 9L, 2L, 3L, 6L,     10L, 0L, 4L, 6L, 9L, 2L, 4L, 6L, 9L, 2L, 3L, 6L, 9L, 1L,     4L, 6L, 9L, 2L, 4L, 6L, 10L, 0L, 4L, 6L, 10L, 1L, 4L, 6L,     9L, 1L, 4L, 6L, 10L, 1L, 4L, 6L, 10L, 1L, 4L, 6L, 9L, 2L,     4L, 6L, 9L, 1L, 4L, 6L, 10L, 2L, 5L, 6L, 10L, 2L, 3L, 6L,     9L, 1L, 4L, 6L, 9L, 1L, 4L, 6L, 10L, 1L, 4L, 6L, 9L, 1L,     4L, 6L, 10L, 1L, 4L, 6L, 10L, 2L, 3L, 6L, 10L, 2L, 4L, 6L,     10L, 2L, 4L, 6L, 10L, 2L, 3L, 7L, 9L, 0L, 4L, 7L, 9L, 1L,     4L, 7L, 9L, 1L, 4L, 7L, 9L, 1L, 4L, 7L, 9L, 0L, 3L, 7L, 10L,     1L, 3L, 7L, 10L, 1L, 4L, 7L, 10L, 2L, 5L, 7L, 9L, 1L, 3L,     7L, 9L, 1L, 3L, 7L, 10L, 1L, 5L, 7L, 9L, 2L, 3L, 7L, 10L,     2L, 4L, 7L, 9L, 1L, 3L, 7L, 10L, 1L, 3L, 7L, 10L, 2L, 4L,     7L, 9L, 2L, 3L, 7L, 9L, 2L, 3L, 7L, 9L, 2L, 4L, 7L, 10L,     2L, 3L, 7L, 10L, 1L, 4L, 7L, 10L, 1L, 4L, 7L, 9L, 2L, 3L,     7L, 10L, 1L, 4L, 7L, 9L, 2L, 3L, 7L, 10L, 2L, 4L, 7L, 9L,     0L, 4L, 7L, 9L, 2L, 4L, 7L, 10L, 1L, 3L, 7L, 9L, 2L, 3L,     7L, 9L, 2L, 3L, 7L, 9L, 1L, 3L, 7L, 10L, 0L, 4L, 7L, 9L,     2L, 4L, 7L, 9L, 2L, 3L, 7L, 9L, 2L, 4L, 7L, 9L, 0L, 4L, 7L,     9L, 1L, 3L, 7L, 10L, 1L, 4L, 7L, 9L, 1L, 3L, 7L, 10L, 1L,     4L, 7L, 10L, 2L, 4L, 7L, 10L, 1L, 4L, 7L, 10L, 2L, 4L, 7L,     9L, 1L, 4L, 7L, 10L, 2L, 4L, 7L, 9L, 2L, 4L, 7L, 10L, 1L,     4L, 7L, 10L, 2L, 3L, 7L, 9L, 2L, 4L, 7L, 10L, 1L, 5L, 7L,     9L, 2L, 3L, 7L, 10L, 2L, 3L, 7L, 10L, 0L, 4L, 7L, 10L, 0L,     4L, 7L, 10L, 1L, 3L, 7L, 10L, 2L, 4L, 7L, 9L, 2L, 4L, 7L,     10L, 1L, 4L, 7L, 10L, 1L, 3L, 7L, 10L, 1L, 3L, 7L, 10L, 0L,     4L, 7L, 9L, 1L, 5L, 7L, 9L, 1L, 4L, 7L, 10L, 2L, 3L, 7L,     9L, 2L, 5L, 7L, 9L, 2L, 3L, 7L, 10L, 1L, 5L, 7L, 9L, 1L,     4L, 7L, 9L, 2L, 4L, 7L, 9L, 1L, 4L, 7L, 9L, 0L, 4L, 7L, 10L,     1L, 4L, 7L, 9L, 1L, 4L, 7L, 9L, 1L, 5L, 7L, 9L, 2L, 3L, 7L,     9L, 2L, 4L, 7L, 10L, 1L, 3L, 7L, 9L, 2L, 3L, 7L, 10L, 2L,     5L, 7L, 10L, 1L, 4L, 7L, 9L, 2L, 3L, 7L, 10L, 2L, 4L, 7L,     10L, 2L, 4L, 8L, 10L, 2L, 3L, 8L, 9L, 1L, 4L, 8L, 9L, 0L,     3L, 8L, 10L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 10L, 2L, 4L, 8L,     9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 1L,     4L, 8L, 10L, 2L, 4L, 8L, 9L, 2L, 3L, 8L, 10L, 1L, 4L, 8L,     10L, 2L, 3L, 8L, 10L, 2L, 3L, 8L, 9L, 1L, 4L, 8L, 9L, 2L,     4L, 8L, 10L, 2L, 3L, 8L, 10L, 2L, 4L, 8L, 10L, 1L, 4L, 8L,     9L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 2L,     3L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 10L,     2L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 3L, 8L,     10L, 1L, 4L, 8L, 10L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 10L, 1L,     4L, 8L, 10L, 1L, 4L, 8L, 10L, 2L, 4L, 8L, 9L, 2L, 4L, 8L,     9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L, 3L, 8L, 10L, 1L,     4L, 8L, 9L, 2L, 4L, 8L, 10L, 2L, 4L, 8L, 10L, 1L, 4L, 8L,     10L, 2L, 4L, 8L, 9L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 2L,     5L, 8L, 10L, 2L, 3L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L,     9L, 1L, 4L, 8L, 9L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L,     8L, 10L, 2L, 3L, 8L, 10L, 2L, 4L, 8L, 10L, 2L, 4L, 8L, 9L,     1L, 4L, 8L, 9L, 2L, 3L, 8L, 9L, 1L, 4L, 8L, 9L, 1L, 4L, 8L,     9L, 2L, 4L, 8L, 9L, 2L, 4L, 8L, 9L, 1L, 4L, 8L, 10L, 1L,     4L, 8L, 10L, 1L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 2L, 4L, 8L,     10L, 2L, 3L, 8L, 10L, 0L, 4L, 8L, 10L, 1L, 4L, 8L, 10L, 2L,     3L, 8L, 10L, 1L, 4L, 8L, 9L), p = c(0L, 4L, 8L, 12L, 16L,     20L, 24L, 28L, 32L, 36L, 40L, 44L, 48L, 52L, 56L, 60L, 64L,     68L, 72L, 76L, 80L, 84L, 88L, 92L, 96L, 100L, 104L, 108L,     112L, 116L, 120L, 124L, 128L, 132L, 136L, 140L, 144L, 148L,     152L, 156L, 160L, 164L, 168L, 172L, 176L, 180L, 184L, 188L,     192L, 196L, 200L, 204L, 208L, 212L, 216L, 220L, 224L, 228L,     232L, 236L, 240L, 244L, 248L, 252L, 256L, 260L, 264L, 268L,     272L, 276L, 280L, 284L, 288L, 292L, 296L, 300L, 304L, 308L,     312L, 316L, 320L, 324L, 328L, 332L, 336L, 340L, 344L, 348L,     352L, 356L, 360L, 364L, 368L, 372L, 376L, 380L, 384L, 388L,     392L, 396L, 400L, 404L, 408L, 412L, 416L, 420L, 424L, 428L,     432L, 436L, 440L, 444L, 448L, 452L, 456L, 460L, 464L, 468L,     472L, 476L, 480L, 484L, 488L, 492L, 496L, 500L, 504L, 508L,     512L, 516L, 520L, 524L, 528L, 532L, 536L, 540L, 544L, 548L,     552L, 556L, 560L, 564L, 568L, 572L, 576L, 580L, 584L, 588L,     592L, 596L, 600L, 604L, 608L, 612L, 616L, 620L, 624L, 628L,     632L, 636L, 640L, 644L, 648L, 652L, 656L, 660L, 664L, 668L,     672L, 676L, 680L, 684L, 688L, 692L, 696L, 700L, 704L, 708L,     712L, 716L, 720L, 724L, 728L, 732L, 736L, 740L, 744L, 748L,     752L, 756L, 760L, 764L, 768L, 772L, 776L, 780L, 784L, 788L,     792L, 796L, 800L, 804L, 808L, 812L, 816L, 820L, 824L, 828L,     832L, 836L, 840L, 844L, 848L, 852L, 856L, 860L, 864L, 868L,     872L, 876L, 880L, 884L, 888L, 892L, 896L), Dim = c(11L, 224L    ), Dimnames = list(c("0", "1", "2", "0", "1", "2", "CRC",     "Healthy", "Polyp", "0", "1"), c("1", "2", "3", "4", "5",     "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16",     "17", "18", "19", "20", "21", "22", "23", "24", "25", "26",     "27", "28", "29", "30", "31", "32", "33", "34", "35", "36",     "37", "38", "39", "40", "41", "42", "43", "44", "45", "46",     "47", "48", "49", "50", "51", "52", "53", "54", "55", "56",     "57", "58", "59", "60", "61", "62", "63", "64", "65", "66",     "67", "68", "69", "70", "71", "72", "73", "74", "75", "76",     "77", "78", "79", "80", "81", "82", "83", "84", "85", "86",     "87", "88", "89", "90", "91", "92", "93", "94", "95", "96",     "97", "98", "99", "100", "101", "102", "103", "104", "105",     "106", "107", "108", "109", "110", "111", "112", "113", "114",     "115", "116", "117", "118", "119", "120", "121", "122", "123",     "124", "125", "126", "127", "128", "129", "130", "131", "132",     "133", "134", "135", "136", "137", "138", "139", "140", "141",     "142", "143", "144", "145", "146", "147", "148", "149", "150",     "151", "152", "153", "154", "155", "156", "157", "158", "159",     "160", "161", "162", "163", "164", "165", "166", "167", "168",     "169", "170", "171", "172", "173", "174", "175", "176", "177",     "178", "179", "180", "181", "182", "183", "184", "185", "186",     "187", "188", "189", "190", "191", "192", "193", "194", "195",     "196", "197", "198", "199", "200", "201", "202", "203", "204",     "205", "206", "207", "208", "209", "210", "211", "212", "213",     "214", "215", "216", "217", "218", "219", "220", "221", "222",     "223", "224")), x = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,     1, 1, 1, 1, 1, 1, 1, 1, 1, 1), factors = list()), theta = c(1, 1, 1, 1), Lind = c(1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 4L, 4L), Gp = c(0L, 3L, 6L, 9L, 11L), lower = c(0, 0, 0, 0), Lambdat = new("dgCMatrix",     i = 0:10, p = 0:11, Dim = c(11L, 11L), Dimnames = list(NULL,         NULL), x = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), factors = list()),     flist = list(alcohol_consumption = c(3L, 2L, 3L, 3L, 2L,     2L, 1L, 3L, 2L, 2L, 3L, 2L, 3L, 2L, 3L, 3L, 2L, 2L, 3L, 3L,     3L, 2L, 3L, 2L, 3L, 2L, 3L, 2L, 2L, 3L, 3L, 2L, 3L, 2L, 3L,     2L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 3L, 3L, 2L, 3L, 1L, 2L, 2L,     2L, 2L, 3L, 2L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L,     1L, 2L, 2L, 2L, 1L, 2L, 2L, 3L, 2L, 2L, 2L, 3L, 3L, 2L, 2L,     3L, 3L, 3L, 3L, 3L, 2L, 2L, 3L, 2L, 3L, 3L, 1L, 3L, 2L, 3L,     3L, 2L, 1L, 3L, 3L, 3L, 1L, 2L, 2L, 2L, 2L, 3L, 2L, 3L, 2L,     3L, 3L, 2L, 3L, 3L, 2L, 3L, 3L, 1L, 1L, 2L, 3L, 3L, 2L, 2L,     2L, 1L, 2L, 2L, 3L, 3L, 3L, 2L, 2L, 3L, 2L, 1L, 2L, 2L, 2L,     3L, 3L, 2L, 3L, 3L, 2L, 3L, 3L, 3L, 3L, 2L, 1L, 3L, 2L, 3L,     2L, 2L, 2L, 2L, 3L, 3L, 2L, 3L, 3L, 2L, 3L, 3L, 3L, 2L, 2L,     3L, 2L, 3L, 2L, 2L, 3L, 3L, 2L, 2L, 3L, 2L, 2L, 3L, 2L, 2L,     3L, 3L, 2L, 2L, 3L, 2L, 3L, 3L, 2L, 3L, 3L, 2L, 3L, 3L, 2L,     2L, 2L, 3L, 2L, 2L, 3L, 3L, 3L, 2L, 3L, 2L, 2L, 3L, 3L, 2L,     2L, 2L, 2L, 3L, 3L, 1L, 2L, 3L, 2L), smoking_condition = c(1L,     2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 3L, 2L, 2L, 1L,     2L, 2L, 3L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 1L,     2L, 2L, 2L, 1L, 2L, 3L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L,     2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 1L, 2L, 2L, 2L, 2L, 2L,     1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 3L, 1L, 1L, 3L,     1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L,     2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L,     2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 3L, 1L, 1L, 2L, 2L, 1L,     2L, 2L, 2L, 1L, 1L, 2L, 3L, 2L, 1L, 3L, 1L, 3L, 2L, 2L, 2L,     2L, 2L, 2L, 3L, 1L, 2L, 1L, 1L, 3L, 2L, 1L, 2L, 2L, 1L, 2L,     1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L,     1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L,     2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 3L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 2L,     2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L), factors = c(1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,     2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L,     3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,     3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,     3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,     3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L,     3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), gender = c(2L,     2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L,     2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 1L,     1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L,     2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 2L,     2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 1L,     2L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 1L,     1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L,     2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L,     1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L,     2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 1L,     2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 2L,     2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 2L,     1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 1L,     1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,     1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L)), cnms = list(        alcohol_consumption = "(Intercept)", smoking_condition = "(Intercept)",         factors = "(Intercept)", gender = "(Intercept)"), Ztlist = list(        `1 | alcohol_consumption` = new("dgCMatrix", i = c(2L,         1L, 2L, 2L, 1L, 1L, 0L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L,         2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L,         2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 0L,         2L, 2L, 1L, 2L, 0L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 1L,         1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 0L, 1L, 1L, 1L, 0L, 1L,         1L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L,         1L, 1L, 2L, 1L, 2L, 2L, 0L, 2L, 1L, 2L, 2L, 1L, 0L, 2L,         2L, 2L, 0L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L,         2L, 2L, 1L, 2L, 2L, 0L, 0L, 1L, 2L, 2L, 1L, 1L, 1L, 0L,         1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 0L, 1L, 1L, 1L, 2L,         2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 0L, 2L, 1L, 2L,         1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 1L,         1L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L,         1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 1L,         2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L,         1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 0L, 1L, 2L, 1L),             p = 0:224, Dim = c(3L, 224L), Dimnames = list(c("0",             "1", "2"), c("1", "2", "3", "4", "5", "6", "7", "8",             "9", "10", "11", "12", "13", "14", "15", "16", "17",             "18", "19", "20", "21", "22", "23", "24", "25", "26",             "27", "28", "29", "30", "31", "32", "33", "34", "35",             "36", "37", "38", "39", "40", "41", "42", "43", "44",             "45", "46", "47", "48", "49", "50", "51", "52", "53",             "54", "55", "56", "57", "58", "59", "60", "61", "62",             "63", "64", "65", "66", "67", "68", "69", "70", "71",             "72", "73", "74", "75", "76", "77", "78", "79", "80",             "81", "82", "83", "84", "85", "86", "87", "88", "89",             "90", "91", "92", "93", "94", "95", "96", "97", "98",             "99", "100", "101", "102", "103", "104", "105", "106",             "107", "108", "109", "110", "111", "112", "113",             "114", "115", "116", "117", "118", "119", "120",             "121", "122", "123", "124", "125", "126", "127",             "128", "129", "130", "131", "132", "133", "134",             "135", "136", "137", "138", "139", "140", "141",             "142", "143", "144", "145", "146", "147", "148",             "149", "150", "151", "152", "153", "154", "155",             "156", "157", "158", "159", "160", "161", "162",             "163", "164", "165", "166", "167", "168", "169",             "170", "171", "172", "173", "174", "175", "176",             "177", "178", "179", "180", "181", "182", "183",             "184", "185", "186", "187", "188", "189", "190",             "191", "192", "193", "194", "195", "196", "197",             "198", "199", "200", "201", "202", "203", "204",             "205", "206", "207", "208", "209", "210", "211",             "212", "213", "214", "215", "216", "217", "218",             "219", "220", "221", "222", "223", "224")), x = c(1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1), factors = list()), `1 | smoking_condition` = new("dgCMatrix",             i = c(0L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L,             1L, 2L, 1L, 1L, 0L, 1L, 1L, 2L, 0L, 1L, 1L, 0L, 1L,             0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 2L,             0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L,             1L, 1L, 1L, 1L, 2L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L,             1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 2L, 0L, 0L, 2L,             0L, 1L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 1L,             0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 1L,             0L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L,             2L, 0L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 2L,             1L, 0L, 2L, 0L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 0L,             1L, 0L, 0L, 2L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 1L, 1L,             1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 1L, 0L,             1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L,             1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L,             1L, 1L, 1L, 1L, 2L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L,             1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,             0L, 1L, 1L, 0L, 1L), p = 0:224, Dim = c(3L, 224L),             Dimnames = list(c("0", "1", "2"), c("1", "2", "3",             "4", "5", "6", "7", "8", "9", "10", "11", "12", "13",             "14", "15", "16", "17", "18", "19", "20", "21", "22",             "23", "24", "25", "26", "27", "28", "29", "30", "31",             "32", "33", "34", "35", "36", "37", "38", "39", "40",             "41", "42", "43", "44", "45", "46", "47", "48", "49",             "50", "51", "52", "53", "54", "55", "56", "57", "58",             "59", "60", "61", "62", "63", "64", "65", "66", "67",             "68", "69", "70", "71", "72", "73", "74", "75", "76",             "77", "78", "79", "80", "81", "82", "83", "84", "85",             "86", "87", "88", "89", "90", "91", "92", "93", "94",             "95", "96", "97", "98", "99", "100", "101", "102",             "103", "104", "105", "106", "107", "108", "109",             "110", "111", "112", "113", "114", "115", "116",             "117", "118", "119", "120", "121", "122", "123",             "124", "125", "126", "127", "128", "129", "130",             "131", "132", "133", "134", "135", "136", "137",             "138", "139", "140", "141", "142", "143", "144",             "145", "146", "147", "148", "149", "150", "151",             "152", "153", "154", "155", "156", "157", "158",             "159", "160", "161", "162", "163", "164", "165",             "166", "167", "168", "169", "170", "171", "172",             "173", "174", "175", "176", "177", "178", "179",             "180", "181", "182", "183", "184", "185", "186",             "187", "188", "189", "190", "191", "192", "193",             "194", "195", "196", "197", "198", "199", "200",             "201", "202", "203", "204", "205", "206", "207",             "208", "209", "210", "211", "212", "213", "214",             "215", "216", "217", "218", "219", "220", "221",             "222", "223", "224")), x = c(1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), factors = list()),         `1 | factors` = new("dgCMatrix", i = c(0L, 0L, 0L, 0L,         0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,         0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,         0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,         0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,         0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,         1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,         1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,         1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,         1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,         1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,         1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,         2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,         2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,         2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,         2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,         2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), p = 0:224, Dim = c(3L,         224L), Dimnames = list(c("CRC", "Healthy", "Polyp"),             c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10",             "11", "12", "13", "14", "15", "16", "17", "18", "19",             "20", "21", "22", "23", "24", "25", "26", "27", "28",             "29", "30", "31", "32", "33", "34", "35", "36", "37",             "38", "39", "40", "41", "42", "43", "44", "45", "46",             "47", "48", "49", "50", "51", "52", "53", "54", "55",             "56", "57", "58", "59", "60", "61", "62", "63", "64",             "65", "66", "67", "68", "69", "70", "71", "72", "73",             "74", "75", "76", "77", "78", "79", "80", "81", "82",             "83", "84", "85", "86", "87", "88", "89", "90", "91",             "92", "93", "94", "95", "96", "97", "98", "99", "100",             "101", "102", "103", "104", "105", "106", "107",             "108", "109", "110", "111", "112", "113", "114",             "115", "116", "117", "118", "119", "120", "121",             "122", "123", "124", "125", "126", "127", "128",             "129", "130", "131", "132", "133", "134", "135",             "136", "137", "138", "139", "140", "141", "142",             "143", "144", "145", "146", "147", "148", "149",             "150", "151", "152", "153", "154", "155", "156",             "157", "158", "159", "160", "161", "162", "163",             "164", "165", "166", "167", "168", "169", "170",             "171", "172", "173", "174", "175", "176", "177",             "178", "179", "180", "181", "182", "183", "184",             "185", "186", "187", "188", "189", "190", "191",             "192", "193", "194", "195", "196", "197", "198",             "199", "200", "201", "202", "203", "204", "205",             "206", "207", "208", "209", "210", "211", "212",             "213", "214", "215", "216", "217", "218", "219",             "220", "221", "222", "223", "224")), x = c(1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,         1, 1, 1, 1, 1, 1), factors = list()), `1 | gender` = new("dgCMatrix",             i = c(1L, 1L, 0L, 1L, 0L, 1L, 0L, 1L, 0L, 0L, 0L,             0L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 1L,             0L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 1L,             1L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L,             1L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 1L,             1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 0L,             1L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L,             1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L,             1L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 1L,             0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L,             1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L,             1L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 1L,             0L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 1L,             1L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L,             1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L,             1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L,             1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L,             1L, 1L, 1L, 1L, 0L), p = 0:224, Dim = c(2L, 224L),             Dimnames = list(c("0", "1"), c("1", "2", "3", "4",             "5", "6", "7", "8", "9", "10", "11", "12", "13",             "14", "15", "16", "17", "18", "19", "20", "21", "22",             "23", "24", "25", "26", "27", "28", "29", "30", "31",             "32", "33", "34", "35", "36", "37", "38", "39", "40",             "41", "42", "43", "44", "45", "46", "47", "48", "49",             "50", "51", "52", "53", "54", "55", "56", "57", "58",             "59", "60", "61", "62", "63", "64", "65", "66", "67",             "68", "69", "70", "71", "72", "73", "74", "75", "76",             "77", "78", "79", "80", "81", "82", "83", "84", "85",             "86", "87", "88", "89", "90", "91", "92", "93", "94",             "95", "96", "97", "98", "99", "100", "101", "102",             "103", "104", "105", "106", "107", "108", "109",             "110", "111", "112", "113", "114", "115", "116",             "117", "118", "119", "120", "121", "122", "123",             "124", "125", "126", "127", "128", "129", "130",             "131", "132", "133", "134", "135", "136", "137",             "138", "139", "140", "141", "142", "143", "144",             "145", "146", "147", "148", "149", "150", "151",             "152", "153", "154", "155", "156", "157", "158",             "159", "160", "161", "162", "163", "164", "165",             "166", "167", "168", "169", "170", "171", "172",             "173", "174", "175", "176", "177", "178", "179",             "180", "181", "182", "183", "184", "185", "186",             "187", "188", "189", "190", "191", "192", "193",             "194", "195", "196", "197", "198", "199", "200",             "201", "202", "203", "204", "205", "206", "207",             "208", "209", "210", "211", "212", "213", "214",             "215", "216", "217", "218", "219", "220", "221",             "222", "223", "224")), x = c(1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,             1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), factors = list())),     nl = c(alcohol_consumption = 3L, smoking_condition = 3L,     factors = 3L, gender = 2L)), REML = FALSE, wmsgs = character(0),     start = NULL, verbose = 0L, control = list(optimizer = "nloptwrap",         restart_edge = TRUE, boundary.tol = 1e-05, calc.derivs = TRUE,         use.last.params = FALSE, checkControl = list(check.nobs.vs.rankZ = "ignore",             check.nobs.vs.nlev = "stop", check.nlev.gtreq.5 = "ignore",             check.nlev.gtr.1 = "stop", check.nobs.vs.nRE = "stop",             check.rankX = "message+drop.cols", check.scaleX = "warning",             check.formula.LHS = "stop"), checkConv = list(check.conv.grad = list(            action = "warning", tol = 0.002, relTol = NULL),             check.conv.singular = list(action = "message", tol = 1e-04),             check.conv.hess = list(action = "warning", tol = 1e-06)),         optCtrl = list()))
10: do.call(mkLmerDevfun, c(lmod, list(start = start, verbose = verbose,     control = control)))
11: lme4::lmer(formula_str, data = data.frame(response = x, fixed_effects,     random_effects), REML = FALSE)
12: lmm_fun(x)
13: FUN(newX[, i], ...)
14: apply(to_lmm, 2, function(x) {    lmm_fun(x)})
15: list2(...)
16: dplyr::bind_rows(apply(to_lmm, 2, function(x) {    lmm_fun(x)}))
17: dplyr::mutate(., feature = colnames(to_lmm))
18: dplyr::select(., feature, dplyr::everything())
19: dplyr::as_tibble(.)
20: dplyr::bind_rows(apply(to_lmm, 2, function(x) {    lmm_fun(x)})) %>% dplyr::mutate(feature = colnames(to_lmm)) %>% dplyr::select(feature,     dplyr::everything()) %>% dplyr::as_tibble()
21: withCallingHandlers(expr, message = function(c) if (inherits(c,     classes)) tryInvokeRestart("muffleMessage"))
22: suppressMessages({    res_lmm <- dplyr::bind_rows(apply(to_lmm, 2, function(x) {        lmm_fun(x)    })) %>% dplyr::mutate(feature = colnames(to_lmm)) %>% dplyr::select(feature,         dplyr::everything()) %>% dplyr::as_tibble()})
23: PomaLMM(.)
24: st000284 %>% PomaLMM()
An irrecoverable exception occurred. R is aborting now ...
* checking for unstated dependencies in ‘tests’ ... OK
* checking tests ...
  Running ‘testthat.R’
 OK
* checking for unstated dependencies in vignettes ... OK
* checking package vignettes ... OK
* checking running R code from vignettes ... SKIPPED
* checking re-building of vignette outputs ... SKIPPED
* checking PDF version of manual ... OK
* DONE

Status: 1 ERROR, 2 NOTEs
See
  ‘/Users/biocbuild/bbs-3.19-bioc/meat/POMA.Rcheck/00check.log’
for details.

##############################################################################
##############################################################################
###
### Running command:
###
###   /Library/Frameworks/R.framework/Resources/bin/R CMD INSTALL POMA
###
##############################################################################
##############################################################################


* installing to library ‘/Library/Frameworks/R.framework/Versions/4.4-x86_64/Resources/library’
* installing *source* package ‘POMA’ ...
** using staged installation
** R
** data
*** moving datasets to lazyload DB
** inst
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
*** copying figures
** building package indices
** installing vignettes
** testing if installed package can be loaded from temporary location
** testing if installed package can be loaded from final location
** testing if installed package keeps a record of temporary installation path
* DONE (POMA)

R Under development (unstable) (2024-03-18 r86148) -- "Unsuffered Consequences"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin20

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(testthat)
> library(POMA)
Welcome to POMA!
Version 1.13.26
POMAShiny app: https://github.com/pcastellanoescuder/POMAShiny
> 
> test_check("POMA")
Rank Product analysis for unpaired case 
 

 done  No genes called significant under class1 < class2 

No genes called significant under class1 > class2 

Rank Product analysis for unpaired case 
 

 done  No genes called significant under class1 < class2 

No genes called significant under class1 > class2 

Rank Product analysis for unpaired case 
 

 done  No genes called significant under class1 < class2 

No genes called significant under class1 > class2 

Rank Product analysis for unpaired case 
 

 done  No genes called significant under class1 < class2 

Table2: Genes called significant under class1 > class2 

[ FAIL 0 | WARN 45 | SKIP 6 | PASS 209 ]

══ Skipped tests (6) ═══════════════════════════════════════════════════════════
• This test is skipped. (6): 'test-PomaLMM.R:3:3', 'test-PomaLMM.R:12:3',
  'test-PomaLMM.R:18:3', 'test-PomaLMM.R:27:3', 'test-PomaLMM.R:34:3',
  'test-PomaLMM.R:40:3'

[ FAIL 0 | WARN 45 | SKIP 6 | PASS 209 ]
> 
> proc.time()
   user  system elapsed 
 75.682   1.944  78.307