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


* using log directory ‘/Users/biocbuild/bbs-3.18-bioc/meat/BUSseq.Rcheck’
* using R version 4.3.2 Patched (2023-11-01 r85457)
* using platform: x86_64-apple-darwin20 (64-bit)
* R was compiled by
    Apple clang version 14.0.3 (clang-1403.0.22.14.1)
    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 ‘BUSseq/DESCRIPTION’ ... OK
* checking extension type ... Package
* this is package ‘BUSseq’ version ‘1.8.0’
* checking package namespace information ... OK
* checking package dependencies ... OK
* checking if this is a source package ... OK
* checking if there is a namespace ... OK
* checking for hidden files and directories ... OK
* checking for portable file names ... OK
* checking for sufficient/correct file permissions ... OK
* checking whether package ‘BUSseq’ can be installed ... OK
* used C++ compiler: ‘Apple clang version 14.0.0 (clang-1400.0.29.202)’
* used SDK: ‘MacOSX11.3.sdk’
* 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 R files for non-ASCII characters ... OK
* checking R files for syntax errors ... OK
* checking whether the package can be loaded ... OK
* checking whether the package can be loaded with stated dependencies ... OK
* checking whether the package can be unloaded cleanly ... OK
* checking whether the namespace can be loaded with stated dependencies ... OK
* checking whether the namespace can be unloaded cleanly ... OK
* checking 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 LazyData ... 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 sizes of PDF files under ‘inst/doc’ ... OK
* checking files in ‘vignettes’ ... OK
* checking examples ... OK
* checking for unstated dependencies in ‘tests’ ... OK
* checking tests ...
  Running ‘BUSseq_example.R’
 OK
* checking for unstated dependencies in vignettes ... OK
* checking package vignettes in ‘inst/doc’ ... OK
* checking running R code from vignettes ... SKIPPED
* checking re-building of vignette outputs ... SKIPPED
* checking PDF version of manual ... OK
* DONE

Status: 1 NOTE
See
  ‘/Users/biocbuild/bbs-3.18-bioc/meat/BUSseq.Rcheck/00check.log’
for details.

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


* installing to library ‘/Library/Frameworks/R.framework/Versions/4.3-x86_64/Resources/library’
* installing *source* package ‘BUSseq’ ...
** using staged installation
** libs
using C++ compiler: ‘Apple clang version 14.0.0 (clang-1400.0.29.202)’
using SDK: ‘MacOSX11.3.sdk’
/Library/Frameworks/R.framework/Resources/share/make/shlib.mk:10: warning: overriding commands for target `BUSseq.so'
Makevars:9: warning: ignoring old commands for target `BUSseq.so'
clang++ -arch x86_64 -std=gnu++17 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c BUSseq.cpp -o BUSseq.o
BUSseq.cpp:3033:6: warning: variable 'count_gene' set but not used [-Wunused-but-set-variable]
        int count_gene = 0;
            ^
BUSseq.cpp:2674:7: warning: variable 'All_Drop' set but not used [-Wunused-but-set-variable]
        bool All_Drop = true;
             ^
2 warnings generated.
clang++ -arch x86_64 -std=gnu++17 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c rngstream.cpp -o rngstream.o
clang++ -arch x86_64 -std=gnu++17 -dynamiclib -Wl,-headerpad_max_install_names -undefined dynamic_lookup -L/Library/Frameworks/R.framework/Resources/lib -L/opt/R/x86_64/lib -o BUSseq.so BUSseq.o rngstream.o -L/Library/Frameworks/R.framework/Resources/lib -lRlapack -L/Library/Frameworks/R.framework/Resources/lib -lRblas -L/opt/gfortran/lib/gcc/x86_64-apple-darwin20.0/12.2.0 -L/opt/gfortran/lib -lgfortran -lquadmath -F/Library/Frameworks/R.framework/.. -framework R -Wl,-framework -Wl,CoreFoundation
installing to /Library/Frameworks/R.framework/Versions/4.3-x86_64/Resources/library/00LOCK-BUSseq/00new/BUSseq/libs
** R
** data
*** moving datasets to lazyload DB
** 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 (BUSseq)

R version 4.3.2 Patched (2023-11-01 r85457) -- "Eye Holes"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin20 (64-bit)

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.

> #######################################
> # Apply BUSseq to the Simulation Data #
> #######################################
> library(BUSseq)
> library(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, sort, 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

> RawCountData <- assay(BUSseqfits_example, "counts")
> batch_ind <- unlist(colData(BUSseqfits_example))
> sce <- SingleCellExperiment(assays = list(counts = RawCountData),
+                             colData = DataFrame(Batch_ind = batch_ind))
> BUSseqfits_res <- BUSseq_MCMC(ObservedData = sce, 
+                               seed = 1234, n.cores = 2,
+                               n.celltypes = 4, n.iterations = 500)
   conducting the posterior sampling...

[>                                       ] Finish 0.00k/0.50k iterations.
[====>                                   ] Finish 0.05k/0.50k iterations.
[========>                               ] Finish 0.10k/0.50k iterations.
[============>                           ] Finish 0.15k/0.50k iterations.
[================>                       ] Finish 0.20k/0.50k iterations.
[====================>                   ] Finish 0.25k/0.50k iterations.
[========================>               ] Finish 0.30k/0.50k iterations.
[============================>           ] Finish 0.35k/0.50k iterations.
[================================>       ] Finish 0.40k/0.50k iterations.
[====================================>   ] Finish 0.45k/0.50k iterations.
[========================================] Finish 0.50k/0.50k iterations.

   The MCMC sampling takes: 0.523 mins

   conducting the posterior inference...

   Posterior inference takes: 0.019 mins

> 
> ################################################
> # Extract Estimates from the BUSseqfits Object #
> ################################################
> 
> #return cell type indicators
> w.est <- celltypes(BUSseqfits_res)
Batch 1 cells' cell type indicators: 1,1,1... ...

Batch 2 cells' cell type indicators: 1,1,1... ...

The output format is an N-dimensional verctor.

> 
> #return the intercept and odds ratio of the logistic regression
> #for dropout events
> gamma.est <- dropout_coefficient_values(BUSseqfits_res)
The output format is a matrix.

Each row represents a batch, the first column corresponds to intercept and the second column is the odd ratio.

> 
> #return the log-scale baseline expression values
> alpha.est <-  baseline_expression_values(BUSseqfits_res)
The output format is a vector.

> 
> #return the cell-type effects
> beta.est <- celltype_effects(BUSseqfits_res)
The output format is a matrix.

Each row represents a gene, and each column corresponds to a cell type.

> 
> #return the mean expression levels
> mu.est <- celltype_mean_expression(BUSseqfits_res)
The output format is a matrix.

Each row represents a gene, and each column corresponds to a cell type.

> 
> #return the cell-specific global effects
> delta.est <- cell_effect_values(BUSseqfits_res)
The output format is an N-dimensional vector.

> 
> #return the location batch effects
> nu.est <- location_batch_effects(BUSseqfits_res)
The output format is a matrix.

Each row represents a gene, and each column corresponds to a batch.

> 
> #return the overdispersion parameters
> phi.est <- overdispersions(BUSseqfits_res)
The output format is a matrix.

Each row represents a gene, and each column corresponds to a batch.

> 
> #return the intrinsic gene indices
> D.est <- intrinsic_genes_BUSseq(BUSseqfits_res)
> 
> #return the BIC value
> BIC <- BIC_BUSseq(BUSseqfits_res)
BIC is 460939.471305478

The output is a scalar.

> 
> #return the raw read count matrix
> CountData_raw <- raw_read_counts(BUSseqfits_res)
The output format is a matrix, in which each row represents a gene and each column does a cell.

> 
> #return the imputed read count matrix
> CountData_imputed <- imputed_read_counts(BUSseqfits_res)
The output format is a matrix, in which each row represents a gene and each column does a cell.

> 
> #return the corrected read count matrix
> BUSseqfits_res <- corrected_read_counts(BUSseqfits_res)
   correcting read counts...

The corrected read count matrix is added into the output "SingleCellExperiment" object.

> 
> #################
> # Visualization #
> #################
> #generate the heatmap of raw read count data
> heatmap_data_BUSseq(BUSseqfits_res, project_name="Heatmap_raw")
null device 
          1 
> 
> #generate the heatmap of imputed read count data
> heatmap_data_BUSseq(BUSseqfits_res, data_type = "Imputed",
+                     project_name="Heatmap_imputed")
null device 
          1 
> 
> #generate the heatmap of corrected read count data
> heatmap_data_BUSseq(BUSseqfits_res, data_type = "Corrected", 
+                     project_name="Heatmap_corrected")
null device 
          1 
> 
> proc.time()
   user  system elapsed 
 42.520   1.155  44.548