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


* using log directory ‘/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck’
* using R version 4.4.0 Patched (2024-04-24 r86482)
* using platform: aarch64-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 Ventura 13.6.6
* using session charset: UTF-8
* using option ‘--no-vignettes’
* checking for file ‘BASiCS/DESCRIPTION’ ... OK
* checking extension type ... Package
* this is package ‘BASiCS’ version ‘2.17.0’
* 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 ‘BASiCS’ can be installed ... WARNING
Found the following significant warnings:
  ./utils.h:357:9: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  ./utils.h:358:10: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  ./updates.h:163:6: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  ./updates.h:258:11: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  ./updates.h:266:13: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  ./updatesNoSpikes.h:76:11: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  ./updatesNoSpikes.h:102:9: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  ./updatesRegNoSpikes.h:107:10: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  ./updatesRegNoSpikes.h:139:9: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  ./MCMCReg.h:383:9: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  ./MCMCRegNoSpikes.h:382:8: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
See ‘/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/00install.out’ for details.
* used C++ compiler: ‘Apple clang version 15.0.0 (clang-1500.0.40.1)’
* used SDK: ‘MacOSX11.3.sdk’
* checking C++ specification ... NOTE
  Specified C++11: please drop specification unless essential
* 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 ... OK
* checking S3 generic/method consistency ... OK
* checking replacement functions ... OK
* checking foreign function calls ... OK
* checking R code for possible problems ... OK
* 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 data for ASCII and uncompressed saves ... OK
* checking line endings in C/C++/Fortran sources/headers ... OK
* checking line endings in Makefiles ... OK
* checking compilation flags in Makevars ... OK
* checking for GNU extensions in Makefiles ... OK
* checking for portable use of $(BLAS_LIBS) and $(LAPACK_LIBS) ... OK
* checking use of PKG_*FLAGS in Makefiles ... OK
* checking compiled code ... NOTE
Note: information on .o files is not available
* checking files in ‘vignettes’ ... OK
* checking examples ... OK
Examples with CPU (user + system) or elapsed time > 5s
                         user system elapsed
BASiCS_MCMC             7.261  0.117   7.412
BASiCS_TestDE           5.829  0.083   5.914
BASiCS_DivideAndConquer 3.921  0.070  31.786
* 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 WARNING, 2 NOTEs
See
  ‘/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/00check.log’
for details.

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


* installing to library ‘/Library/Frameworks/R.framework/Versions/4.4-arm64/Resources/library’
* installing *source* package ‘BASiCS’ ...
** using staged installation
** libs
using C++ compiler: ‘Apple clang version 15.0.0 (clang-1500.0.40.1)’
using C++11
using SDK: ‘MacOSX11.3.sdk’
clang++ -arch arm64 -std=gnu++11 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG  -I'/Library/Frameworks/R.framework/Versions/4.4-arm64/Resources/library/Rcpp/include' -I'/Library/Frameworks/R.framework/Versions/4.4-arm64/Resources/library/RcppArmadillo/include' -I/opt/R/arm64/include    -fPIC  -falign-functions=64 -Wall -g -O2   -c BASiCS.cpp -o BASiCS.o
In file included from BASiCS.cpp:1:
In file included from ./updates.h:4:
./utils.h:357:9: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
    if( arma::is_finite(log_aux(i)) & arma::is_finite(y(i)) ) {
        ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                                    &&
./utils.h:357:9: note: cast one or both operands to int to silence this warning
./utils.h:358:10: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
      if((log(u(i)) < log_aux(i)) & (y(i) > threshold)) { ind(i) = 1; }
         ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                                  &&
./utils.h:358:10: note: cast one or both operands to int to silence this warning
In file included from BASiCS.cpp:1:
./updates.h:163:6: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  if(all(prop_var * phi1 < 2.5327372760800758e+305)  & 
     ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
./updates.h:163:6: note: cast one or both operands to int to silence this warning
./updates.h:163:6: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  if(all(prop_var * phi1 < 2.5327372760800758e+305)  & 
     ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
./updates.h:163:6: note: cast one or both operands to int to silence this warning
./updates.h:163:6: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
  if(all(prop_var * phi1 < 2.5327372760800758e+305)  & 
     ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                                                     &&
./updates.h:163:6: note: cast one or both operands to int to silence this warning
./updates.h:258:11: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
      if (!R_IsNA(chi(j)) & (chi(j) > 0)) {
          ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                          &&
./updates.h:258:11: note: cast one or both operands to int to silence this warning
./updates.h:266:13: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
        if (!(chi(j) < 0) & (lambda(j) > 0)) {
            ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                          &&
./updates.h:266:13: note: cast one or both operands to int to silence this warning
In file included from BASiCS.cpp:3:
./updatesNoSpikes.h:76:11: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
      if ((log(u(iAux)) < log_aux(iAux)) & (mu1(iAux) > mintol)) {
          ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                                         &&
./updatesNoSpikes.h:76:11: note: cast one or both operands to int to silence this warning
./updatesNoSpikes.h:102:9: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
    if ((log(u(iAux)) < log_aux(iAux)) & (mu1(iAux) > mintol)) {
        ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                                       &&
./updatesNoSpikes.h:102:9: note: cast one or both operands to int to silence this warning
In file included from BASiCS.cpp:4:
./updatesRegNoSpikes.h:107:10: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
      if((log(u(iAux)) < log_aux(iAux)) & (mu1(iAux) > mintol)) {
         ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                                        &&
./updatesRegNoSpikes.h:107:10: note: cast one or both operands to int to silence this warning
./updatesRegNoSpikes.h:139:9: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
    if ((log(u(iAux)) < log_aux(iAux)) & (mu1(iAux) > mintol)) {
        ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                                       &&
./updatesRegNoSpikes.h:139:9: note: cast one or both operands to int to silence this warning
In file included from BASiCS.cpp:6:
./MCMCReg.h:383:9: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
    if ((det(V1) != 0) & all(arma::eig_sym(sigma2Aux * VAux) > 0)) {
        ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                       &&
./MCMCReg.h:383:9: note: cast one or both operands to int to silence this warning
In file included from BASiCS.cpp:8:
./MCMCRegNoSpikes.h:382:8: warning: use of bitwise '&' with boolean operands [-Wbitwise-instead-of-logical]
    if((det(V1)!=0) & all(arma::eig_sym(sigma2Aux * VAux) > 0)) {
       ^~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
                    &&
./MCMCRegNoSpikes.h:382:8: note: cast one or both operands to int to silence this warning
13 warnings generated.
clang++ -arch arm64 -std=gnu++11 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG  -I'/Library/Frameworks/R.framework/Versions/4.4-arm64/Resources/library/Rcpp/include' -I'/Library/Frameworks/R.framework/Versions/4.4-arm64/Resources/library/RcppArmadillo/include' -I/opt/R/arm64/include    -fPIC  -falign-functions=64 -Wall -g -O2   -c RcppExports.cpp -o RcppExports.o
clang++ -arch arm64 -std=gnu++11 -dynamiclib -Wl,-headerpad_max_install_names -undefined dynamic_lookup -L/Library/Frameworks/R.framework/Resources/lib -L/opt/R/arm64/lib -o BASiCS.so BASiCS.o RcppExports.o -L/Library/Frameworks/R.framework/Resources/lib -lRlapack -L/Library/Frameworks/R.framework/Resources/lib -lRblas -L/opt/gfortran/lib/gcc/aarch64-apple-darwin20.0/12.2.0 -L/opt/gfortran/lib -lgfortran -lemutls_w -lquadmath -F/Library/Frameworks/R.framework/.. -framework R -Wl,-framework -Wl,CoreFoundation
installing to /Library/Frameworks/R.framework/Versions/4.4-arm64/Resources/library/00LOCK-BASiCS/00new/BASiCS/libs
** R
** data
** inst
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
** building package indices
** installing vignettes
** testing if installed package can be loaded from temporary location
** checking absolute paths in shared objects and dynamic libraries
** testing if installed package can be loaded from final location
** testing if installed package keeps a record of temporary installation path
* DONE (BASiCS)

R version 4.4.0 Patched (2024-04-24 r86482) -- "Puppy Cup"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: aarch64-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(BASiCS)
Loading required package: SingleCellExperiment
Loading required package: SummarizedExperiment
Loading required package: MatrixGenerics
Loading required package: matrixStats

Attaching package: 'MatrixGenerics'

The following objects are masked from 'package:matrixStats':

    colAlls, colAnyNAs, colAnys, colAvgsPerRowSet, colCollapse,
    colCounts, colCummaxs, colCummins, colCumprods, colCumsums,
    colDiffs, colIQRDiffs, colIQRs, colLogSumExps, colMadDiffs,
    colMads, colMaxs, colMeans2, colMedians, colMins, colOrderStats,
    colProds, colQuantiles, colRanges, colRanks, colSdDiffs, colSds,
    colSums2, colTabulates, colVarDiffs, colVars, colWeightedMads,
    colWeightedMeans, colWeightedMedians, colWeightedSds,
    colWeightedVars, rowAlls, rowAnyNAs, rowAnys, rowAvgsPerColSet,
    rowCollapse, rowCounts, rowCummaxs, rowCummins, rowCumprods,
    rowCumsums, rowDiffs, rowIQRDiffs, rowIQRs, rowLogSumExps,
    rowMadDiffs, rowMads, rowMaxs, rowMeans2, rowMedians, rowMins,
    rowOrderStats, rowProds, rowQuantiles, rowRanges, rowRanks,
    rowSdDiffs, rowSds, rowSums2, rowTabulates, rowVarDiffs, rowVars,
    rowWeightedMads, rowWeightedMeans, rowWeightedMedians,
    rowWeightedSds, rowWeightedVars

Loading required package: GenomicRanges
Loading required package: stats4
Loading required package: BiocGenerics

Attaching package: 'BiocGenerics'

The following objects are masked from 'package:stats':

    IQR, mad, sd, var, xtabs

The following objects are masked from 'package:base':

    Filter, Find, Map, Position, Reduce, anyDuplicated, aperm, append,
    as.data.frame, basename, cbind, colnames, dirname, do.call,
    duplicated, eval, evalq, get, grep, grepl, intersect, is.unsorted,
    lapply, mapply, match, mget, order, paste, pmax, pmax.int, pmin,
    pmin.int, rank, rbind, rownames, sapply, setdiff, table, tapply,
    union, unique, unsplit, which.max, which.min

Loading required package: S4Vectors

Attaching package: 'S4Vectors'

The following object is masked from 'package:utils':

    findMatches

The following objects are masked from 'package:base':

    I, expand.grid, unname

Loading required package: IRanges
Loading required package: GenomeInfoDb
Loading required package: Biobase
Welcome to Bioconductor

    Vignettes contain introductory material; view with
    'browseVignettes()'. To cite Bioconductor, see
    'citation("Biobase")', and for packages 'citation("pkgname")'.


Attaching package: 'Biobase'

The following object is masked from 'package:MatrixGenerics':

    rowMedians

The following objects are masked from 'package:matrixStats':

    anyMissing, rowMedians

    Welcome to 'BASiCS'. If you used 'BASiCS' before its release in
    Bioconductor, please visit:
    https://github.com/catavallejos/BASiCS/wiki.
> 
> test_check("BASiCS")
-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6714
(ratio SC vs P&S).
-------------------------------------------------------------

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6714
(ratio SC vs P&S).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6714
(ratio Group1 vs Group2).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

For Differential mean task:
the posterior probability threshold chosen via EFDR calibration is too low.
Probability threshold automatically set equal to 'ProbThresholdM'.
EFDR calibration failed for Differential dispersion task. 
Probability threshold automatically set equal to 'ProbThresholdD'.
-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6714
(ratio SC vs P&S).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6714
(ratio Group1 vs Group2).
-------------------------------------------------------------

For Differential mean task:
the posterior probability threshold chosen via EFDR calibration is too low.
Probability threshold automatically set equal to 'ProbThresholdM'.
EFDR calibration failed for Differential dispersion task. 
Probability threshold automatically set equal to 'ProbThresholdD'.
-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6714
(ratio SC vs P&S).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6714
(ratio SC vs P&S).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6787
(ratio SC vs P&S).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6787
(ratio Group1 vs Group2).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6787
(ratio Group1 vs Group2).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

EFDR calibration failed for Differential mean task. 
Probability threshold automatically set equal to 'ProbThresholdM'.
EFDR calibration failed for Differential dispersion task. 
Probability threshold automatically set equal to 'ProbThresholdD'.
For Differential residual dispersion task:
the posterior probability threshold chosen via EFDR calibration is too low.
Probability threshold automatically set equal to 'ProbThresholdR'.
-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6787
(ratio SC vs P&S).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6787
(ratio Group1 vs Group2).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

EFDR calibration failed for Differential mean task. 
Probability threshold automatically set equal to 'ProbThresholdM'.
EFDR calibration failed for Differential dispersion task. 
Probability threshold automatically set equal to 'ProbThresholdD'.
For Differential residual dispersion task:
the posterior probability threshold chosen via EFDR calibration is too low.
Probability threshold automatically set equal to 'ProbThresholdR'.
altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.01

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.011
system: 0
elapsed: 0.01

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.009

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.01
system: 0
elapsed: 0.011

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for phi
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
Combining batch posteriors for beta
Combining batch posteriors for sigma2
Combining batch posteriors for epsilon
Combining batch posteriors for RBFLocations
altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.008

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.009

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for phi
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.009

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.01

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
Combining batch posteriors for beta
Combining batch posteriors for sigma2
Combining batch posteriors for epsilon
Combining batch posteriors for RBFLocations
-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.007
system: 0
elapsed: 0.007

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.008

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.01
system: 0
elapsed: 0.011

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.003
system: 0
elapsed: 0.003

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for phi
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
Combining batch posteriors for beta
Combining batch posteriors for sigma2
Combining batch posteriors for epsilon
Combining batch posteriors for RBFLocations
altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.008
system: 0
elapsed: 0.007

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.009

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for phi
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.01
system: 0
elapsed: 0.01

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.009

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
Combining batch posteriors for beta
Combining batch posteriors for sigma2
Combining batch posteriors for epsilon
Combining batch posteriors for RBFLocations
-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.007
system: 0
elapsed: 0.006

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.008

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.014
system: 0
elapsed: 0.014

-------------------------------------------------------------
Output 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.059
system: 0
elapsed: 0.058

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.063
system: 0
elapsed: 0.063

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for phi
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
Combining batch posteriors for beta
Combining batch posteriors for sigma2
Combining batch posteriors for epsilon
Combining batch posteriors for RBFLocations
Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for phi
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
Combining batch posteriors for beta
Combining batch posteriors for sigma2
Combining batch posteriors for epsilon
Combining batch posteriors for RBFLocations
Combining batch posteriors for mu
Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for phi
Combining batch posteriors for s
Combining batch posteriors for nu
Combining batch posteriors for theta
Cannot find a balanced split with 9 quantiles, trying again with 8
altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.063
system: 0
elapsed: 0.064

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.055
system: 0
elapsed: 0.055

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

Generating partitions...
Starting MCMC...
-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-----------------------------------------------------
MCMC sampler has been started: 200 iterations to go.
-----------------------------------------------------
RefGene: 13
RefGene: 44
RefGene: 31
RefGene: 43
RefGene: 0
RefGene: 16
RefGene: 13
RefGene: 17
RefGene: 15
RefGene: 17
RefGene: 16
RefGene: 44
RefGene: 16
RefGene: 13
RefGene: 43
RefGene: 39
RefGene: 11
RefGene: 43
RefGene: 39
RefGene: 31
RefGene: 17
RefGene: 44
RefGene: 44
RefGene: 39
RefGene: 11
RefGene: 13
RefGene: 13
RefGene: 31
RefGene: 16
RefGene: 43
RefGene: 31
RefGene: 43
RefGene: 0
RefGene: 43
RefGene: 43
RefGene: 11
RefGene: 17
RefGene: 44
RefGene: 13
RefGene: 11
RefGene: 11
RefGene: 16
RefGene: 13
RefGene: 11
RefGene: 0
RefGene: 11
RefGene: 13
RefGene: 39
RefGene: 44
RefGene: 16
RefGene: 39
RefGene: 0
RefGene: 11
RefGene: 17
RefGene: 44
RefGene: 15
RefGene: 11
RefGene: 17
RefGene: 31
RefGene: 15
RefGene: 15
RefGene: 39
RefGene: 15
RefGene: 0
RefGene: 43
RefGene: 13
RefGene: 39
RefGene: 31
RefGene: 44
RefGene: 11
RefGene: 17
RefGene: 39
RefGene: 0
RefGene: 44
RefGene: 13
RefGene: 39
RefGene: 43
RefGene: 43
RefGene: 16
RefGene: 31
RefGene: 15
RefGene: 31
RefGene: 31
RefGene: 44
RefGene: 39
RefGene: 16
RefGene: 43
RefGene: 44
RefGene: 0
RefGene: 11
RefGene: 17
RefGene: 31
RefGene: 0
RefGene: 31
RefGene: 0
RefGene: 16
RefGene: 43
RefGene: 31
RefGene: 31
RefGene: 16
-----------------------------------------------------
End of Burn-in period.
-----------------------------------------------------
RefGene: 44
RefGene: 17
RefGene: 16
RefGene: 11
RefGene: 16
RefGene: 0
RefGene: 44
RefGene: 44
RefGene: 13
RefGene: 13
RefGene: 15
RefGene: 11
RefGene: 17
RefGene: 16
RefGene: 43
RefGene: 39
RefGene: 0
RefGene: 11
RefGene: 39
RefGene: 16
RefGene: 0
RefGene: 0
RefGene: 15
RefGene: 13
RefGene: 11
RefGene: 0
RefGene: 15
RefGene: 0
RefGene: 0
RefGene: 0
RefGene: 0
RefGene: 39
RefGene: 15
RefGene: 16
RefGene: 17
RefGene: 13
RefGene: 11
RefGene: 39
RefGene: 11
RefGene: 39
RefGene: 43
RefGene: 44
RefGene: 11
RefGene: 16
RefGene: 15
RefGene: 11
RefGene: 0
RefGene: 17
RefGene: 16
RefGene: 44
RefGene: 17
RefGene: 43
RefGene: 0
RefGene: 44
RefGene: 17
RefGene: 0
RefGene: 0
RefGene: 16
RefGene: 31
RefGene: 13
RefGene: 44
RefGene: 15
RefGene: 44
RefGene: 0
RefGene: 0
RefGene: 16
RefGene: 31
RefGene: 13
RefGene: 31
RefGene: 11
RefGene: 16
RefGene: 11
RefGene: 0
RefGene: 16
RefGene: 0
RefGene: 0
RefGene: 43
RefGene: 17
RefGene: 39
RefGene: 31
RefGene: 39
RefGene: 43
RefGene: 39
RefGene: 15
RefGene: 39
RefGene: 13
RefGene: 17
RefGene: 15
RefGene: 44
RefGene: 31
RefGene: 16
RefGene: 15
RefGene: 39
RefGene: 17
RefGene: 16
RefGene: 0
RefGene: 13
RefGene: 15
RefGene: 44
RefGene: 39
 
-----------------------------------------------------
-----------------------------------------------------
All 200 MCMC iterations have been completed.
-----------------------------------------------------
-----------------------------------------------------
 
-----------------------------------------------------
Please see below a summary of the overall acceptance rates.
-----------------------------------------------------
 
Minimum acceptance rate among mu[i]'s: 0.16
Average acceptance rate among mu[i]'s: 0.496
Maximum acceptance rate among mu[i]'s: 0.76
 
 
Minimum acceptance rate among delta[i]'s: 0.54
Average acceptance rate among delta[i]'s: 0.681
Maximum acceptance rate among delta[i]'s: 0.84
 
 
Minimum acceptance rate among nu[jk]'s: 0.84
Average acceptance rate among nu[jk]'s: 0.9528
Maximum acceptance rate among nu[jk]'s: 1
 
 
Minimum acceptance rate among theta[k]'s: 1
Average acceptance rate among theta[k]'s: 1
Maximum acceptance rate among theta[k]'s: 1
 
-----------------------------------------------------
 
-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.181
system: 0
elapsed: 0.181

-------------------------------------------------------------
Output 
-------------------------------------------------------------


-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-----------------------------------------------------
MCMC sampler has been started: 200 iterations to go.
-----------------------------------------------------
RefGene: 14
RefGene: 49
RefGene: 14
RefGene: 9
RefGene: 14
RefGene: 47
RefGene: 36
RefGene: 10
RefGene: 49
RefGene: 9
RefGene: 11
RefGene: 49
RefGene: 19
RefGene: 14
RefGene: 23
RefGene: 10
RefGene: 14
RefGene: 22
RefGene: 10
RefGene: 9
RefGene: 47
RefGene: 9
RefGene: 19
RefGene: 14
RefGene: 11
RefGene: 23
RefGene: 22
RefGene: 19
RefGene: 23
RefGene: 23
RefGene: 22
RefGene: 10
RefGene: 11
RefGene: 49
RefGene: 11
RefGene: 22
RefGene: 47
RefGene: 47
RefGene: 47
RefGene: 19
RefGene: 49
RefGene: 10
RefGene: 19
RefGene: 19
RefGene: 22
RefGene: 11
RefGene: 14
RefGene: 10
RefGene: 10
RefGene: 19
RefGene: 36
RefGene: 47
RefGene: 23
RefGene: 19
RefGene: 19
RefGene: 11
RefGene: 22
RefGene: 36
RefGene: 11
RefGene: 14
RefGene: 49
RefGene: 19
RefGene: 36
RefGene: 36
RefGene: 10
RefGene: 10
RefGene: 47
RefGene: 23
RefGene: 47
RefGene: 11
RefGene: 19
RefGene: 19
RefGene: 36
RefGene: 47
RefGene: 14
RefGene: 49
RefGene: 9
RefGene: 22
RefGene: 9
RefGene: 36
RefGene: 23
RefGene: 9
RefGene: 23
RefGene: 19
RefGene: 47
RefGene: 10
RefGene: 14
RefGene: 10
RefGene: 19
RefGene: 23
RefGene: 9
RefGene: 11
RefGene: 49
RefGene: 19
RefGene: 19
RefGene: 47
RefGene: 49
RefGene: 23
RefGene: 14
RefGene: 47
-----------------------------------------------------
End of Burn-in period.
-----------------------------------------------------
RefGene: 14
RefGene: 19
RefGene: 10
RefGene: 14
RefGene: 47
RefGene: 23
RefGene: 19
RefGene: 49
RefGene: 49
RefGene: 47
RefGene: 49
RefGene: 10
RefGene: 47
RefGene: 9
RefGene: 9
RefGene: 36
RefGene: 49
RefGene: 10
RefGene: 47
RefGene: 36
RefGene: 19
RefGene: 10
RefGene: 14
RefGene: 11
RefGene: 47
RefGene: 10
RefGene: 22
RefGene: 23
RefGene: 11
RefGene: 49
RefGene: 9
RefGene: 14
RefGene: 10
RefGene: 19
RefGene: 36
RefGene: 9
RefGene: 49
RefGene: 49
RefGene: 23
RefGene: 49
RefGene: 19
RefGene: 47
RefGene: 22
RefGene: 36
RefGene: 36
RefGene: 47
RefGene: 9
RefGene: 10
RefGene: 36
RefGene: 22
RefGene: 11
RefGene: 47
RefGene: 19
RefGene: 36
RefGene: 14
RefGene: 14
RefGene: 49
RefGene: 19
RefGene: 19
RefGene: 14
RefGene: 47
RefGene: 36
RefGene: 47
RefGene: 11
RefGene: 23
RefGene: 19
RefGene: 22
RefGene: 14
RefGene: 22
RefGene: 10
RefGene: 11
RefGene: 19
RefGene: 11
RefGene: 11
RefGene: 19
RefGene: 19
RefGene: 47
RefGene: 19
RefGene: 23
RefGene: 49
RefGene: 49
RefGene: 14
RefGene: 10
RefGene: 11
RefGene: 14
RefGene: 36
RefGene: 36
RefGene: 19
RefGene: 14
RefGene: 22
RefGene: 19
RefGene: 47
RefGene: 22
RefGene: 23
RefGene: 22
RefGene: 47
RefGene: 47
RefGene: 10
RefGene: 47
RefGene: 47
 
-----------------------------------------------------
-----------------------------------------------------
All 200 MCMC iterations have been completed.
-----------------------------------------------------
-----------------------------------------------------
 
-----------------------------------------------------
Please see below a summary of the overall acceptance rates.
-----------------------------------------------------
 
Minimum acceptance rate among mu[i]'s: 0.13
Average acceptance rate among mu[i]'s: 0.5164
Maximum acceptance rate among mu[i]'s: 0.81
 
 
Minimum acceptance rate among delta[i]'s: 0.55
Average acceptance rate among delta[i]'s: 0.6708
Maximum acceptance rate among delta[i]'s: 0.84
 
 
Minimum acceptance rate among nu[jk]'s: 0.84
Average acceptance rate among nu[jk]'s: 0.9603
Maximum acceptance rate among nu[jk]'s: 1
 
 
Minimum acceptance rate among theta[k]'s: 0.59
Average acceptance rate among theta[k]'s: 0.795
Maximum acceptance rate among theta[k]'s: 1
 
-----------------------------------------------------
 
-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.176
system: 0
elapsed: 0.176

-------------------------------------------------------------
Output 
-------------------------------------------------------------


Combining batch posteriors for mu
Combining batch posteriors for delta
Combining batch posteriors for s

Combining batch posteriors for nu
Combining batch posteriors for theta

sumAux (line 117): 146.424
sumAux (line 117): 136.459
For LVG detection task:
the posterior probability threshold chosen via EFDR calibration is too low.
Probability threshold automatically set equal to 'ProbThreshold'.
For LVG detection task:
the posterior probability threshold chosen via EFDR calibration is too low.
Probability threshold automatically set equal to 'ProbThreshold'.
For LVG detection task:
the posterior probability threshold chosen via EFDR calibration is too low.
Probability threshold automatically set equal to 'ProbThreshold'.
Evaluating variance contribution threshold = 25 % ... 

Evaluating variance contribution threshold = 50 % ... 

Evaluating variance contribution threshold = 75 % ... 

Evaluating variance contribution threshold = 25 % ... 

Evaluating variance contribution threshold = 50 % ... 

Evaluating variance contribution threshold = 75 % ... 

For LVG detection task:
the posterior probability threshold chosen via EFDR calibration is too low.
Probability threshold automatically set equal to 'ProbThreshold'.
altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.006

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.005

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.01
system: 0
elapsed: 0.01

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.004
system: 0
elapsed: 0.004

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.004
system: 0
elapsed: 0.004

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.007
system: 0
elapsed: 0.007

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.007
system: 0
elapsed: 0.007

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.004
system: 0
elapsed: 0.003

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.003
system: 0
elapsed: 0.004

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.004
system: 0
elapsed: 0.004

-------------------------------------------------------------
Output 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.005

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.001
system: 0
elapsed: 0.002

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.002
system: 0
elapsed: 0.002

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.005

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.005

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.006

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.006
system: 0
elapsed: 0.006

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.006
system: 0
elapsed: 0.007

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.007
system: 0
elapsed: 0.007

-------------------------------------------------------------
Output 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.004
system: 0
elapsed: 0.004

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.006

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.006
system: 0
elapsed: 0.006

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.007
system: 0
elapsed: 0.008

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS_Chain object stored as chain_.Rds file in
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
Storing trajectories of adaptive proposal variances (log-scale) aschain_ls_.Rds file in 
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.009
system: 0
elapsed: 0.01

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS_Chain object stored as chain_.Rds file in
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
Storing trajectories of adaptive proposal variances (log-scale) aschain_ls_.Rds file in 
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.007
system: 0
elapsed: 0.006

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS_Chain object stored as chain_.Rds file in
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
Storing trajectories of adaptive proposal variances (log-scale) aschain_ls_.Rds file in 
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : horizontal integration (no-spikes case)
-------------------------------------------------------------

Randomly, 1 out of 10
genes was left as reference at each iteration 
List of reference genes and frequencies stored as a .txt file in
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.01
system: 0
elapsed: 0.009

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS_Chain object stored as chain_.Rds file in
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
Storing trajectories of adaptive proposal variances (log-scale) aschain_ls_.Rds file in 
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : horizontal integration (no-spikes case)
-------------------------------------------------------------

Randomly, 1 out of 10
genes was left as reference at each iteration 
List of reference genes and frequencies stored as a .txt file in
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.008
system: 0
elapsed: 0.007

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

`BASiCS_Chain` class definition was outdated. 
Object updated to be compatible with BASiCS version 2.17.0.
Set 'StoreUpdatedChain' = TRUE to save updated object.

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.008
system: 0
elapsed: 0.007

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

`BASiCS_Chain` class definition was outdated. 
Object updated to be compatible with BASiCS version 2.17.0.
Set 'StoreUpdatedChain' = TRUE to save updated object.
(this overwrites original input file).

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 1
(ratio Group1 vs Group2).
-------------------------------------------------------------

EFDR calibration failed for Differential mean task. 
Probability threshold automatically set equal to 'ProbThresholdM'.
EFDR calibration failed for Differential dispersion task. 
Probability threshold automatically set equal to 'ProbThresholdD'.
EFDR calibration failed for Differential residual dispersion task. 
Probability threshold automatically set equal to 'ProbThresholdR'.
Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.022
system: 0.001
elapsed: 0.022

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 1
(ratio Group1 vs Group2).
-------------------------------------------------------------

EFDR calibration failed for Differential mean task. 
Probability threshold automatically set equal to 'ProbThresholdM'.
EFDR calibration failed for Differential dispersion task. 
Probability threshold automatically set equal to 'ProbThresholdD'.
EFDR calibration failed for Differential residual dispersion task. 
Probability threshold automatically set equal to 'ProbThresholdR'.
altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.335
system: 0.001
elapsed: 0.336

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.003
system: 0
elapsed: 0.003

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
Storing trajectories of adaptive proposal variances (log-scale) aschain_ls_.Rds file in 
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.34
system: 0
elapsed: 0.341

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.307
system: 0
elapsed: 0.307

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.307
system: 0.001
elapsed: 0.307

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.003
system: 0
elapsed: 0.003

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
Storing trajectories of adaptive proposal variances (log-scale) aschain_ls_.Rds file in 
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.307
system: 0.004
elapsed: 0.311

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.305
system: 0.001
elapsed: 0.305

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.56
system: 0.012
elapsed: 0.572

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.556
system: 0.011
elapsed: 0.566

-------------------------------------------------------------
Output 
-------------------------------------------------------------

Running no spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.022
system: 0
elapsed: 0.022

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
Storing trajectories of adaptive proposal variances (log-scale) aschain_ls_.Rds file in 
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.445
system: 0.001
elapsed: 0.446

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.004

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
Storing trajectories of adaptive proposal variances (log-scale) aschain_ls_.Rds file in 
'/Users/biocbuild/bbs-3.20-bioc/meat/BASiCS.Rcheck/tests/testthat' directory ... 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6714
(ratio Group1 vs Group2).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

For Differential mean task:
the posterior probability threshold chosen via EFDR calibration is too low.
Probability threshold automatically set equal to 'ProbThresholdM'.
EFDR calibration failed for Differential dispersion task. 
Probability threshold automatically set equal to 'ProbThresholdD'.
-------------------------------------------------------------
Log-fold change thresholds are now set in a log2 scale. 
Original BASiCS release used a natural logarithm scale.
-------------------------------------------------------------
Offset estimate: 0.6787
(ratio Group1 vs Group2).
To visualise its effect, please use 'PlotOffset = TRUE'.
-------------------------------------------------------------

EFDR calibration failed for Differential mean task. 
Probability threshold automatically set equal to 'ProbThresholdM'.
EFDR calibration failed for Differential dispersion task. 
Probability threshold automatically set equal to 'ProbThresholdD'.
For Differential residual dispersion task:
the posterior probability threshold chosen via EFDR calibration is too low.
Probability threshold automatically set equal to 'ProbThresholdR'.
altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.005

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.006

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.001
system: 0
elapsed: 0.002

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

Running with spikes BASiCS sampler (regression case) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.005
system: 0
elapsed: 0.004

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

'BatchInfo' was treated as a factor with 2 levels: 
1  2  
'BatchInfo' was treated as a factor with 2 levels: 
1  2  
-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running no spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.004
system: 0
elapsed: 0.004

-------------------------------------------------------------
Output 
-------------------------------------------------------------

'BatchInfo' was treated as a factor with 2 levels: 
1  2  
altExp 'spike-ins' is assumed to contain spike-in genes.
see help(altExp) for details. 

-------------------------------------------------------------
NOTE: default choice PriorDelta = 'log-normal'  (recommended value). 
Vallejos et al (2015) used a 'gamma' prior instead.
-------------------------------------------------------------

Running with spikes BASiCS sampler (no regression) ... 

-------------------------------------------------------------
MCMC running time 
-------------------------------------------------------------
user: 0.004
system: 0
elapsed: 0.004

-------------------------------------------------------------
Output 
-------------------------------------------------------------

-------------------------------------------------------------
BASiCS version 2.17.0 : 
vertical integration (spikes case) 
-------------------------------------------------------------

[ FAIL 0 | WARN 0 | SKIP 0 | PASS 459 ]
> 
> proc.time()
   user  system elapsed 
241.992   4.891 297.772