##############################################################################
##############################################################################
###
### Running command:
###
###   /home/biocbuild/bbs-3.17-bioc/R/bin/R CMD check --install=check:Biobase.install-out.txt --library=/home/biocbuild/bbs-3.17-bioc/R/library --timings Biobase_2.59.0.tar.gz
###
##############################################################################
##############################################################################


* using log directory ‘/home/biocbuild/bbs-3.17-bioc/meat/Biobase.Rcheck’
* using R Under development (unstable) (2023-01-14 r83615)
* using platform: aarch64-unknown-linux-gnu (64-bit)
* R was compiled by
    gcc (Ubuntu 11.3.0-1ubuntu1~22.04) 11.3.0
    GNU Fortran (Ubuntu 11.3.0-1ubuntu1~22.04) 11.3.0
* running under: Ubuntu 22.04.1 LTS
* using session charset: UTF-8
* checking for file ‘Biobase/DESCRIPTION’ ... OK
* this is package ‘Biobase’ version ‘2.59.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 ‘Biobase’ can be installed ... OK
* used C compiler: ‘gcc (Ubuntu 11.3.0-1ubuntu1~22.04) 11.3.0’
* 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 ... WARNING
Invalid citation information in ‘inst/CITATION’:
  Warning message:
  In utils::packageDescription(basename(dir), dirname(dir)) :
    no package 'Biobase' was found
* checking R files for non-ASCII characters ... OK
* checking R files for syntax errors ... OK
* checking whether the package can be loaded ... OK
* checking whether the package can be loaded with stated dependencies ... OK
* checking whether the package can be unloaded cleanly ... OK
* checking whether the namespace can be loaded with stated dependencies ... OK
* checking whether the namespace can be unloaded cleanly ... OK
* checking loading without being on the library search path ... WARNING
Error in library(Biobase, lib.loc = "/home/biocbuild/bbs-3.17-bioc/R/library") : 
  there is no package called ‘Biobase’
Execution halted

It looks like this package has a loading problem when not on .libPaths:
see the messages for details.
* 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 compiled code ... OK
* checking files in ‘vignettes’ ... OK
* checking examples ... SKIPPED
* checking for unstated dependencies in ‘tests’ ... OK
* checking tests ...
  Running ‘test-all.R’
  Running ‘test-rowMedians.R’
 OK
* checking for unstated dependencies in vignettes ... OK
* checking package vignettes in ‘inst/doc’ ... OK
* checking running R code from vignettes ...
  ‘BiobaseDevelopment.Rmd’ using ‘UTF-8’... OK
  ‘esApply.Rmd’ using ‘UTF-8’... OK
  ‘ExpressionSetIntroduction.Rnw’... OK
 OK
* checking re-building of vignette outputs ... OK
* checking PDF version of manual ... OK
* DONE

Status: 2 WARNINGs
See
  ‘/home/biocbuild/bbs-3.17-bioc/meat/Biobase.Rcheck/00check.log’
for details.

##############################################################################
##############################################################################
###
### Running command:
###
###   /home/biocbuild/bbs-3.17-bioc/R/bin/R CMD INSTALL Biobase
###
##############################################################################
##############################################################################


* installing to library ‘/home/biocbuild/bbs-3.17-bioc/R/site-library’
* installing *source* package ‘Biobase’ ...
** using staged installation
** libs
using C compiler: ‘gcc (Ubuntu 11.3.0-1ubuntu1~22.04) 11.3.0’
gcc -I"/home/biocbuild/bbs-3.17-bioc/R/include" -DNDEBUG   -I/usr/local/include    -fPIC  -g -O2  -Wall -c Rinit.c -o Rinit.o
gcc -I"/home/biocbuild/bbs-3.17-bioc/R/include" -DNDEBUG   -I/usr/local/include    -fPIC  -g -O2  -Wall -c anyMissing.c -o anyMissing.o
gcc -I"/home/biocbuild/bbs-3.17-bioc/R/include" -DNDEBUG   -I/usr/local/include    -fPIC  -g -O2  -Wall -c envir.c -o envir.o
gcc -I"/home/biocbuild/bbs-3.17-bioc/R/include" -DNDEBUG   -I/usr/local/include    -fPIC  -g -O2  -Wall -c matchpt.c -o matchpt.o
gcc -I"/home/biocbuild/bbs-3.17-bioc/R/include" -DNDEBUG   -I/usr/local/include    -fPIC  -g -O2  -Wall -c rowMedians.c -o rowMedians.o
gcc -I"/home/biocbuild/bbs-3.17-bioc/R/include" -DNDEBUG   -I/usr/local/include    -fPIC  -g -O2  -Wall -c sublist_extract.c -o sublist_extract.o
gcc -shared -L/home/biocbuild/bbs-3.17-bioc/R/lib -L/usr/local/lib -o Biobase.so Rinit.o anyMissing.o envir.o matchpt.o rowMedians.o sublist_extract.o -L/home/biocbuild/bbs-3.17-bioc/R/lib -lR
installing to /home/biocbuild/bbs-3.17-bioc/R/site-library/00LOCK-Biobase/00new/Biobase/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 (Biobase)

R Under development (unstable) (2023-01-14 r83615) -- "Unsuffered Consequences"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: aarch64-unknown-linux-gnu (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.

> BiocGenerics:::testPackage("Biobase")

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

Welcome to Bioconductor

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



RUNIT TEST PROTOCOL -- Sat Jan 28 18:17:52 2023 
*********************************************** 
Number of test functions: 101 
Number of errors: 0 
Number of failures: 0 

 
1 Test Suite : 
Biobase RUnit Tests - 101 test functions, 0 errors, 0 failures
Number of test functions: 101 
Number of errors: 0 
Number of failures: 0 
> 
> proc.time()
   user  system elapsed 
 13.457   0.122  14.058 

R Under development (unstable) (2023-01-14 r83615) -- "Unsuffered Consequences"
Copyright (C) 2023 The R Foundation for Statistical Computing
Platform: aarch64-unknown-linux-gnu (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.

> library(Biobase)
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

Welcome to Bioconductor

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

> set.seed(1)
> 
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Consistency checks
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> 
> # rowMedians() by rowQ()
> rowMedians2 <- function(imat) {
+   nr <- ncol(imat)
+   half <- (nr + 1)/2
+   if (nr%%2 == 1) {
+     return(rowQ(imat, half))
+   } else {
+     return((rowQ(imat, half) + rowQ(imat, half+1))/2)
+   }
+ }
> 
> cat("Consistency checks:\n")
Consistency checks:
> set.seed(1)
> for (kk in 1:20) {
+   cat("Random test #", kk, "\n", sep="")
+ 
+   # Simulate data in a matrix of any shape
+   nrow <- sample(2000, size=1)
+   ncol <- sample(2000, size=1)
+   x <- rnorm(nrow*ncol)
+   dim(x) <- c(nrow, ncol)
+ 
+   # Add NAs?
+   nas <- sample(c(TRUE,FALSE), size=1)
+   if (nas) {
+     nna <- sample(nrow*ncol, size=1)
+     x[sample(length(x), size=nna)] <- NA
+   }
+ 
+   na.rm <- nas
+   t1 <- system.time({
+     y1 <- rowMedians(x, na.rm=na.rm)
+   })
+   t2 <- system.time({
+     y2 <- apply(x, MARGIN=1, FUN=median, na.rm=na.rm)
+   })
+   # When all values of 'y2' are NA, 'y2' is logical
+   if (is.logical(y2)) y2 <- as.double(y2)
+   stopifnot(all.equal(y1,y2))
+   cat(sprintf("rowMedians()/apply(): %.3g\n", (t1/t2)[3]))
+ 
+   if (!nas) {
+     t3 <- system.time({
+       y3 <- rowMedians2(x)
+     })
+     stopifnot(all.equal(y1,y3))
+     cat(sprintf("rowMedians()/rowMedians2(): %.3g\n", (t1/t3)[3]))
+   }
+ }
Random test #1
rowMedians()/apply(): 0.33
Random test #2
rowMedians()/apply(): 0.13
rowMedians()/rowMedians2(): 0.283
Random test #3
rowMedians()/apply(): 0.0476
Random test #4
rowMedians()/apply(): 0.175
rowMedians()/rowMedians2(): 0.273
Random test #5
rowMedians()/apply(): 0.307
rowMedians()/rowMedians2(): 0.387
Random test #6
rowMedians()/apply(): 0.257
rowMedians()/rowMedians2(): 0.346
Random test #7
rowMedians()/apply(): 0.377
Random test #8
rowMedians()/apply(): 0
rowMedians()/rowMedians2(): NaN
Random test #9
rowMedians()/apply(): 0.167
Random test #10
rowMedians()/apply(): 0.066
Random test #11
rowMedians()/apply(): 0.154
Random test #12
rowMedians()/apply(): 0.119
Random test #13
rowMedians()/apply(): 0.227
rowMedians()/rowMedians2(): 0.625
Random test #14
rowMedians()/apply(): 0.247
Random test #15
rowMedians()/apply(): 0.31
Random test #16
rowMedians()/apply(): 0.138
rowMedians()/rowMedians2(): 0.291
Random test #17
rowMedians()/apply(): 0.27
Random test #18
rowMedians()/apply(): 0.111
rowMedians()/rowMedians2(): 0.571
Random test #19
rowMedians()/apply(): 0.385
rowMedians()/rowMedians2(): 0.754
Random test #20
rowMedians()/apply(): 0.303
rowMedians()/rowMedians2(): 0.378
> 
> 
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Benchmarking
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> cat("Benchmarking:\n")
Benchmarking:
> 
> # Simulate data in a matrix of any shape
> nrow <- 1000
> ncol <- 1000
> x <- rnorm(nrow*ncol)
> dim(x) <- c(nrow, ncol)
> 
> gc()
          used (Mb) gc trigger (Mb) max used (Mb)
Ncells  581034 31.1    1327274 70.9  1044722 55.8
Vcells 2061511 15.8   12255594 93.6 12287603 93.8
> t0 <- system.time({
+   for (rr in 1:20)
+     y0 <- apply(x, MARGIN=1, FUN=median, na.rm=FALSE)
+ })
> gc()
          used (Mb) gc trigger (Mb) max used (Mb)
Ncells  581050 31.1    1327274 70.9  1211070 64.7
Vcells 2062528 15.8   12255594 93.6 12287603 93.8
> t1 <- system.time({
+   for (rr in 1:20)
+     y1 <- rowMedians(x, na.rm=FALSE)
+ })
> gc()
          used (Mb) gc trigger (Mb) max used (Mb)
Ncells  581050 31.1    1327274 70.9  1211070 64.7
Vcells 2061653 15.8   12255594 93.6 12287603 93.8
> stopifnot(all.equal(y0,y1))
> cat(sprintf("rowMedians()/apply(): %.3g\n", (t1/t0)[3]))
rowMedians()/apply(): 0.178
> 
> 
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> # Consistency checks
> # - - - - - - - - - - - - - - - - - - - - - - - - - - - -
> cat("Consistency checks without NAs:\n")
Consistency checks without NAs:
> for (kk in 1:20) {
+   cat("Random test #", kk, "\n", sep="")
+ 
+   # Simulate data in a matrix of any shape
+   nrow <- sample(1000, size=1)
+   ncol <- sample(1000, size=1)
+   x <- rnorm(nrow*ncol)
+   dim(x) <- c(nrow, ncol)
+ 
+   t0 <- system.time({
+     y0 <- apply(x, MARGIN=1, FUN=median, na.rm=FALSE)
+   })
+   t1 <- system.time({
+     y1 <- rowMedians(x, na.rm=FALSE)
+   })
+   stopifnot(all.equal(y0,y1))
+ } # for (kk in ...)
Random test #1
Random test #2
Random test #3
Random test #4
Random test #5
Random test #6
Random test #7
Random test #8
Random test #9
Random test #10
Random test #11
Random test #12
Random test #13
Random test #14
Random test #15
Random test #16
Random test #17
Random test #18
Random test #19
Random test #20
> 
> 
> 
> cat("Consistency checks with NAs:\n")
Consistency checks with NAs:
> for (kk in 1:20) {
+   cat("Random test #", kk, "\n", sep="")
+ 
+   # Simulate data in a matrix of any shape
+   nrow <- sample(1000, size=1)
+   ncol <- sample(1000, size=1)
+   x <- rnorm(nrow*ncol)
+   dim(x) <- c(nrow, ncol)
+ 
+   # Add NAs
+   nna <- sample(nrow*ncol-1, size=1)
+   x[sample(length(x), size=nna)] <- NA
+ 
+   t0 <- system.time({
+     y0 <- apply(x, MARGIN=1, FUN=median, na.rm=TRUE)
+     y0[is.na(y0)] <- NA
+   })
+   t1 <- system.time({
+     y1 <- rowMedians(x, na.rm=TRUE)
+   })
+   stopifnot(all.equal(y0,y1))
+ } # for (kk in ...)
Random test #1
Random test #2
Random test #3
Random test #4
Random test #5
Random test #6
Random test #7
Random test #8
Random test #9
Random test #10
Random test #11
Random test #12
Random test #13
Random test #14
Random test #15
Random test #16
Random test #17
Random test #18
Random test #19
Random test #20
> 
> proc.time()
   user  system elapsed 
 21.819   0.208  22.329