| Back to Build/check report for BioC 3.18: simplified long |
|
This page was generated on 2023-05-10 10:04:37 -0000 (Wed, 10 May 2023).
| Hostname | OS | Arch (*) | R version | Installed pkgs |
|---|---|---|---|---|
| kunpeng1 | Linux (Ubuntu 22.04.1 LTS) | aarch64 | R Under development (unstable) (2023-03-12 r83975) -- "Unsuffered Consequences" | 6211 |
| Click on any hostname to see more info about the system (e.g. compilers) (*) as reported by 'uname -p', except on Windows and Mac OS X | ||||
|
To the developers/maintainers of the preprocessCore package: - Allow up to 24 hours (and sometimes 48 hours) for your latest push to git@git.bioconductor.org:packages/preprocessCore.git to reflect on this report. See Troubleshooting Build Report for more information. - Use the following Renviron settings to reproduce errors and warnings. Note: If "R CMD check" recently failed on the Linux builder over a missing dependency, add the missing dependency to "Suggests" in your DESCRIPTION file. See the Renviron.bioc for details. |
| Package 1525/2194 | Hostname | OS / Arch | INSTALL | BUILD | CHECK | BUILD BIN | ||||||||
| preprocessCore 1.63.1 (landing page) Ben Bolstad
| kunpeng1 | Linux (Ubuntu 22.04.1 LTS) / aarch64 | OK | OK | WARNINGS | |||||||||
| Package: preprocessCore |
| Version: 1.63.1 |
| Command: /home/biocbuild/bbs-3.17-bioc/R/bin/R CMD check --install=check:preprocessCore.install-out.txt --library=/home/biocbuild/bbs-3.17-bioc/R/site-library --timings preprocessCore_1.63.1.tar.gz |
| StartedAt: 2023-05-10 00:51:17 -0000 (Wed, 10 May 2023) |
| EndedAt: 2023-05-10 00:51:34 -0000 (Wed, 10 May 2023) |
| EllapsedTime: 17.0 seconds |
| RetCode: 0 |
| Status: WARNINGS |
| CheckDir: preprocessCore.Rcheck |
| Warnings: 1 |
##############################################################################
##############################################################################
###
### Running command:
###
### /home/biocbuild/bbs-3.17-bioc/R/bin/R CMD check --install=check:preprocessCore.install-out.txt --library=/home/biocbuild/bbs-3.17-bioc/R/site-library --timings preprocessCore_1.63.1.tar.gz
###
##############################################################################
##############################################################################
* using log directory ‘/home/biocbuild/bbs-3.17-bioc/meat/preprocessCore.Rcheck’
* using R Under development (unstable) (2023-03-12 r83975)
* 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.2 LTS
* using session charset: UTF-8
* checking for file ‘preprocessCore/DESCRIPTION’ ... OK
* this is package ‘preprocessCore’ version ‘1.63.1’
* 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 ‘preprocessCore’ can be installed ... OK
* used C compiler: ‘gcc (Ubuntu 11.3.0-1ubuntu1~22.04) 11.3.0’
* checking installed package size ... OK
* checking package directory ... OK
* checking DESCRIPTION meta-information ... OK
* checking top-level files ... NOTE
Found a ‘configure.in’ file: ‘configure.ac’ has long been preferred.
* 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 ... WARNING
checkRd: (7) normalize.quantiles.Rd:47: Invalid email address: bmbolstad.com
checkRd: (7) rma.background.correct.Rd:35: Invalid email address: bmbolstad.com
* 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 line endings in shell scripts ... 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 examples ... OK
* checking for unstated dependencies in ‘tests’ ... OK
* checking tests ...
Running ‘PLMdtest.R’
Running ‘qnormtest.R’
OK
* checking PDF version of manual ... OK
* DONE
Status: 1 WARNING, 2 NOTEs
See
‘/home/biocbuild/bbs-3.17-bioc/meat/preprocessCore.Rcheck/00check.log’
for details.
preprocessCore.Rcheck/00install.out
############################################################################## ############################################################################## ### ### Running command: ### ### /home/biocbuild/bbs-3.17-bioc/R/bin/R CMD INSTALL preprocessCore ### ############################################################################## ############################################################################## * installing to library ‘/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/site-library’ * installing *source* package ‘preprocessCore’ ... ** using staged installation 'config' variable 'CPP' is defunct checking for gcc... gcc checking whether the C compiler works... yes checking for C compiler default output file name... a.out checking for suffix of executables... checking whether we are cross compiling... no checking for suffix of object files... o checking whether we are using the GNU C compiler... yes checking whether gcc accepts -g... yes checking for gcc option to accept ISO C89... none needed checking how to run the C preprocessor... gcc -E checking for library containing pthread_create... none required checking for grep that handles long lines and -e... /usr/bin/grep checking for egrep... /usr/bin/grep -E checking for ANSI C header files... yes checking for sys/types.h... yes checking for sys/stat.h... yes checking for stdlib.h... yes checking for string.h... yes checking for memory.h... yes checking for strings.h... yes checking for inttypes.h... yes checking for stdint.h... yes checking for unistd.h... yes checking for stdlib.h... (cached) yes checking if PTHREAD_STACK_MIN is defined... yes checking if R is using flexiblas... flexiblas not found. preprocessCore threading will not be disabled configure: Enabling threading for preprocessCore configure: creating ./config.status config.status: creating src/Makevars ** libs using C compiler: ‘gcc (Ubuntu 11.3.0-1ubuntu1~22.04) 11.3.0’ gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c R_colSummarize.c -o R_colSummarize.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c R_plmd_interfaces.c -o R_plmd_interfaces.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c R_plmr_interfaces.c -o R_plmr_interfaces.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c R_rlm_interfaces.c -o R_rlm_interfaces.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c R_subColSummarize.c -o R_subColSummarize.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c R_subrcModel_interfaces.c -o R_subrcModel_interfaces.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c avg.c -o avg.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c avg_log.c -o avg_log.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c biweight.c -o biweight.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c init_package.c -o init_package.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c lm.c -o lm.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c log_avg.c -o log_avg.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c log_median.c -o log_median.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c matrix_functions.c -o matrix_functions.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c median.c -o median.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c median_log.c -o median_log.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c medianpolish.c -o medianpolish.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c plmd.c -o plmd.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c plmr.c -o plmr.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c psi_fns.c -o psi_fns.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c qnorm.c -o qnorm.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c rlm.c -o rlm.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c rlm_anova.c -o rlm_anova.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c rlm_se.c -o rlm_se.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c rma_background4.c -o rma_background4.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c rma_common.c -o rma_common.o gcc -I"/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/include" -DNDEBUG -I/usr/local/include -I/usr/local/include -g -O2 -DPACKAGE_NAME=\"\" -DPACKAGE_TARNAME=\"\" -DPACKAGE_VERSION=\"\" -DPACKAGE_STRING=\"\" -DPACKAGE_BUGREPORT=\"\" -DPACKAGE_URL=\"\" -DSTDC_HEADERS=1 -DHAVE_SYS_TYPES_H=1 -DHAVE_SYS_STAT_H=1 -DHAVE_STDLIB_H=1 -DHAVE_STRING_H=1 -DHAVE_MEMORY_H=1 -DHAVE_STRINGS_H=1 -DHAVE_INTTYPES_H=1 -DHAVE_STDINT_H=1 -DHAVE_UNISTD_H=1 -DHAVE_STDLIB_H=1 -DUSE_PTHREADS=1 -fPIC -g -O2 -c weightedkerneldensity.c -o weightedkerneldensity.o gcc -shared -L/usr/local/lib -o preprocessCore.so R_colSummarize.o R_plmd_interfaces.o R_plmr_interfaces.o R_rlm_interfaces.o R_subColSummarize.o R_subrcModel_interfaces.o avg.o avg_log.o biweight.o init_package.o lm.o log_avg.o log_median.o matrix_functions.o median.o median_log.o medianpolish.o plmd.o plmr.o psi_fns.o qnorm.o rlm.o rlm_anova.o rlm_se.o rma_background4.o rma_common.o weightedkerneldensity.o -llapack -L/home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/lib/R/lib -lRblas -lgfortran -lm installing to /home/biocbuild/bbs-3.17-bioc/R-devel_2023-03-12_r83975-bin/site-library/00LOCK-preprocessCore/00new/preprocessCore/libs ** R ** inst ** byte-compile and prepare package for lazy loading ** help *** installing help indices ** building package indices ** 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 (preprocessCore)
preprocessCore.Rcheck/tests/PLMdtest.Rout
R Under development (unstable) (2023-03-12 r83975) -- "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(preprocessCore)
>
>
> values <- rnorm(100)
> group.labels <- sample(0:4,replace=TRUE, 100)
>
> results <- double(10000)
> ngroups <- 2
>
>
> for (i in 1:10000){
+ values <- rnorm(100,sd=1)
+ values <- values/sd(values)
+ group.labels <- sample(0:(ngroups-1),replace=TRUE, 100)
+ blah <- .C("R_split_test",as.double(values), as.integer(100), as.integer(ngroups), as.integer(group.labels),double(1))
+ results[i] <- blah[[5]]
+ }
>
> plot(sort(results),qchisq(0:9999/10000,ngroups-1))
> lm(qchisq(0:9999/10000,ngroups-1) ~ sort(results))
Call:
lm(formula = qchisq(0:9999/10000, ngroups - 1) ~ sort(results))
Coefficients:
(Intercept) sort(results)
0.007376 0.951762
>
>
>
> boxplot(values ~ group.labels,ylim=c(-2,2))
>
>
>
> sc <- median(abs(resid(lm(values ~ 1))))/0.6745
> sum((resid(lm(values ~ 1))/sc)^2)/2
[1] 45.71354
> sum((resid(lm(values ~ as.factor(group.labels)))/sc)^2)/2
[1] 44.49536
>
>
> values <- rnorm(100)
> group.labels <- sample(0:4,replace=TRUE, 100)
> values[group.labels == 1] <- values[group.labels == 1] + 0.4
>
>
> blah <- .C("R_split_test",as.double(values), as.integer(100), as.integer(5), as.integer(group.labels),double(1))
>
> boxplot(values ~ group.labels,ylim=c(-2,2))
>
>
>
> library(preprocessCore)
>
> .C("R_test_get_design_matrix",as.integer(4),as.integer(5))
1.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00
1.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00
1.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
1.00 0.00 0.00 0.00 0.00 -1.00 -1.00 -1.00
0.00 1.00 0.00 0.00 0.00 1.00 0.00 0.00
0.00 1.00 0.00 0.00 0.00 0.00 1.00 0.00
0.00 1.00 0.00 0.00 0.00 0.00 0.00 1.00
0.00 1.00 0.00 0.00 0.00 -1.00 -1.00 -1.00
0.00 0.00 1.00 0.00 0.00 1.00 0.00 0.00
0.00 0.00 1.00 0.00 0.00 0.00 1.00 0.00
0.00 0.00 1.00 0.00 0.00 0.00 0.00 1.00
0.00 0.00 1.00 0.00 0.00 -1.00 -1.00 -1.00
0.00 0.00 0.00 1.00 0.00 1.00 0.00 0.00
0.00 0.00 0.00 1.00 0.00 0.00 1.00 0.00
0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 1.00 0.00 -1.00 -1.00 -1.00
0.00 0.00 0.00 0.00 1.00 1.00 0.00 0.00
0.00 0.00 0.00 0.00 1.00 0.00 1.00 0.00
0.00 0.00 0.00 0.00 1.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 1.00 -1.00 -1.00 -1.00
1.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00
1.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00
1.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
1.00 0.00 0.00 0.00 0.00 -1.00 -1.00 -1.00 -1.00
0.00 1.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00
0.00 1.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00
0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
0.00 1.00 0.00 0.00 0.00 -1.00 -1.00 -1.00 -1.00
0.00 0.00 1.00 0.00 0.00 1.00 0.00 0.00 0.00
0.00 0.00 1.00 0.00 0.00 0.00 0.00 1.00 0.00
0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 1.00
0.00 0.00 1.00 0.00 0.00 -1.00 -1.00 -1.00 -1.00
0.00 0.00 0.00 1.00 0.00 1.00 0.00 0.00 0.00
0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00 0.00
0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 1.00 0.00 -1.00 -1.00 -1.00 -1.00
0.00 0.00 0.00 0.00 1.00 1.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 1.00 0.00 0.00 1.00 0.00
0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 1.00 -1.00 -1.00 -1.00 -1.00
1.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00
1.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00
1.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00
1.00 0.00 0.00 0.00 0.00 -1.00 -1.00 -1.00 -1.00 -1.00
0.00 1.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00
0.00 1.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00
0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00
0.00 1.00 0.00 0.00 0.00 -1.00 -1.00 -1.00 -1.00 -1.00
0.00 0.00 1.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00
0.00 0.00 1.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00
0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00
0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 1.00 0.00 1.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00 0.00 0.00
0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 1.00 0.00
0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 1.00
0.00 0.00 0.00 0.00 1.00 1.00 0.00 0.00 0.00 0.00
0.00 0.00 0.00 0.00 1.00 0.00 0.00 1.00 0.00 0.00
0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 1.00 0.00
0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 1.00
[[1]]
[1] 4
[[2]]
[1] 5
>
>
>
> chips <- as.factor(rep(c(1,2,3,4,5,6),c(5,5,5,5,5,5)))
> probes <- rep(c(1,3,4,5,6),6)
>
> probes[c(1,6,11)] <- 2
> ##probes[24 + c(8,16,24)] <- 10
> probes <- as.factor(probes)
>
>
> model.matrix(~ -1 + probes)%*%contr.sum(6)
[,1] [,2] [,3] [,4] [,5]
1 0 1 0 0 0
2 0 0 1 0 0
3 0 0 0 1 0
4 0 0 0 0 1
5 -1 -1 -1 -1 -1
6 0 1 0 0 0
7 0 0 1 0 0
8 0 0 0 1 0
9 0 0 0 0 1
10 -1 -1 -1 -1 -1
11 0 1 0 0 0
12 0 0 1 0 0
13 0 0 0 1 0
14 0 0 0 0 1
15 -1 -1 -1 -1 -1
16 1 0 0 0 0
17 0 0 1 0 0
18 0 0 0 1 0
19 0 0 0 0 1
20 -1 -1 -1 -1 -1
21 1 0 0 0 0
22 0 0 1 0 0
23 0 0 0 1 0
24 0 0 0 0 1
25 -1 -1 -1 -1 -1
26 1 0 0 0 0
27 0 0 1 0 0
28 0 0 0 1 0
29 0 0 0 0 1
30 -1 -1 -1 -1 -1
>
>
> probes <- rep(c(1,3,4,5,6),6)
>
> probes[c(1,6,11)] <- 2
> probes[c(20,25,30)] <- 7
> probes <- as.factor(probes)
> model.matrix(~ -1 + probes)%*%contr.sum(7)
[,1] [,2] [,3] [,4] [,5] [,6]
1 0 1 0 0 0 0
2 0 0 1 0 0 0
3 0 0 0 1 0 0
4 0 0 0 0 1 0
5 0 0 0 0 0 1
6 0 1 0 0 0 0
7 0 0 1 0 0 0
8 0 0 0 1 0 0
9 0 0 0 0 1 0
10 0 0 0 0 0 1
11 0 1 0 0 0 0
12 0 0 1 0 0 0
13 0 0 0 1 0 0
14 0 0 0 0 1 0
15 0 0 0 0 0 1
16 1 0 0 0 0 0
17 0 0 1 0 0 0
18 0 0 0 1 0 0
19 0 0 0 0 1 0
20 -1 -1 -1 -1 -1 -1
21 1 0 0 0 0 0
22 0 0 1 0 0 0
23 0 0 0 1 0 0
24 0 0 0 0 1 0
25 -1 -1 -1 -1 -1 -1
26 1 0 0 0 0 0
27 0 0 1 0 0 0
28 0 0 0 1 0 0
29 0 0 0 0 1 0
30 -1 -1 -1 -1 -1 -1
>
>
>
>
> probes <- rep(c(1,3,4,5,6),6)
>
> probes[c(1,6,11)] <- 2
> probes[c(5,10,15)] <- 7
> probes <- as.factor(probes)
> model.matrix(~ -1 + probes)%*%contr.sum(7)
[,1] [,2] [,3] [,4] [,5] [,6]
1 0 1 0 0 0 0
2 0 0 1 0 0 0
3 0 0 0 1 0 0
4 0 0 0 0 1 0
5 -1 -1 -1 -1 -1 -1
6 0 1 0 0 0 0
7 0 0 1 0 0 0
8 0 0 0 1 0 0
9 0 0 0 0 1 0
10 -1 -1 -1 -1 -1 -1
11 0 1 0 0 0 0
12 0 0 1 0 0 0
13 0 0 0 1 0 0
14 0 0 0 0 1 0
15 -1 -1 -1 -1 -1 -1
16 1 0 0 0 0 0
17 0 0 1 0 0 0
18 0 0 0 1 0 0
19 0 0 0 0 1 0
20 0 0 0 0 0 1
21 1 0 0 0 0 0
22 0 0 1 0 0 0
23 0 0 0 1 0 0
24 0 0 0 0 1 0
25 0 0 0 0 0 1
26 1 0 0 0 0 0
27 0 0 1 0 0 0
28 0 0 0 1 0 0
29 0 0 0 0 1 0
30 0 0 0 0 0 1
>
>
>
> probes <- rep(c(1,3,4,5,6),6)
>
> probes[c(1,6,11)] <- 2
> probes[1+c(1,6,11)] <- 8
> probes[2+c(1,6,11)] <- 9
> probes[3+c(1,6,11)] <- 10
> probes[c(5,10,15)] <- 7
> probes <- as.factor(probes)
> model.matrix(~ -1 + probes)%*%contr.sum(10)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9]
1 0 1 0 0 0 0 0 0 0
2 0 0 0 0 0 0 0 1 0
3 0 0 0 0 0 0 0 0 1
4 -1 -1 -1 -1 -1 -1 -1 -1 -1
5 0 0 0 0 0 0 1 0 0
6 0 1 0 0 0 0 0 0 0
7 0 0 0 0 0 0 0 1 0
8 0 0 0 0 0 0 0 0 1
9 -1 -1 -1 -1 -1 -1 -1 -1 -1
10 0 0 0 0 0 0 1 0 0
11 0 1 0 0 0 0 0 0 0
12 0 0 0 0 0 0 0 1 0
13 0 0 0 0 0 0 0 0 1
14 -1 -1 -1 -1 -1 -1 -1 -1 -1
15 0 0 0 0 0 0 1 0 0
16 1 0 0 0 0 0 0 0 0
17 0 0 1 0 0 0 0 0 0
18 0 0 0 1 0 0 0 0 0
19 0 0 0 0 1 0 0 0 0
20 0 0 0 0 0 1 0 0 0
21 1 0 0 0 0 0 0 0 0
22 0 0 1 0 0 0 0 0 0
23 0 0 0 1 0 0 0 0 0
24 0 0 0 0 1 0 0 0 0
25 0 0 0 0 0 1 0 0 0
26 1 0 0 0 0 0 0 0 0
27 0 0 1 0 0 0 0 0 0
28 0 0 0 1 0 0 0 0 0
29 0 0 0 0 1 0 0 0 0
30 0 0 0 0 0 1 0 0 0
>
>
>
>
>
>
>
>
>
> true.probes <- c(4,3,2,1,-1,-2,-3,-4)
>
> true.chips <- c(8,9,10,11,12,13)
>
>
> y <- outer(true.probes,true.chips,"+")
>
>
>
> estimate.coefficients <- function(y){
+
+
+ colmean <- apply(y,2,mean)
+
+ y <- sweep(y,2,FUN="-",colmean)
+
+ rowmean <- apply(y,1,mean)
+ y <- sweep(y,1,FUN="-",rowmean)
+
+
+ list(y,colmean,rowmean)
+ }
> estimate.coefficients(y)
[[1]]
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 0 0 0 0 0 0
[2,] 0 0 0 0 0 0
[3,] 0 0 0 0 0 0
[4,] 0 0 0 0 0 0
[5,] 0 0 0 0 0 0
[6,] 0 0 0 0 0 0
[7,] 0 0 0 0 0 0
[8,] 0 0 0 0 0 0
[[2]]
[1] 8 9 10 11 12 13
[[3]]
[1] 4 3 2 1 -1 -2 -3 -4
>
>
>
> y <- outer(true.probes,true.chips,"+")
>
>
> estimate.coefficients(y)
[[1]]
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 0 0 0 0 0 0
[2,] 0 0 0 0 0 0
[3,] 0 0 0 0 0 0
[4,] 0 0 0 0 0 0
[5,] 0 0 0 0 0 0
[6,] 0 0 0 0 0 0
[7,] 0 0 0 0 0 0
[8,] 0 0 0 0 0 0
[[2]]
[1] 8 9 10 11 12 13
[[3]]
[1] 4 3 2 1 -1 -2 -3 -4
>
>
>
>
> y2 <- sweep(y,2,FUN="-",apply(y,2,mean))
>
>
>
> c(3.875, 2.875, 1.875, 0.875,
+ -1.125, -2.125, -3.125, -4, -2.25)
[1] 3.875 2.875 1.875 0.875 -1.125 -2.125 -3.125 -4.000 -2.250
>
>
>
>
> cp <- rep(c(1,2,3,4,5,6),rep(8,6))
> pr <- rep(c(1,2,3,4,5,6,7,8),6)
>
>
> pr[c(32,40,48)] <- 9
>
>
>
>
> true.probes <- c(4,3,2,1,-1,-2,-3,-4)
>
> true.chips <- c(8,9,10,11,12,10)
>
>
> y <- outer(true.probes,true.chips,"+") + rnorm(48,0,0.1)
>
> y[8,4:6] <- c(11,12,10)+2 + rnorm(3,0,0.1)
>
>
> lm(as.vector(y) ~ -1 + as.factor(cp) + C(as.factor(pr),"contr.sum"))
Call:
lm(formula = as.vector(y) ~ -1 + as.factor(cp) + C(as.factor(pr),
"contr.sum"))
Coefficients:
as.factor(cp)1 as.factor(cp)2
8.1946 9.2152
as.factor(cp)3 as.factor(cp)4
10.1711 11.2037
as.factor(cp)5 as.factor(cp)6
12.2991 10.2854
C(as.factor(pr), "contr.sum")1 C(as.factor(pr), "contr.sum")2
3.7743 2.8126
C(as.factor(pr), "contr.sum")3 C(as.factor(pr), "contr.sum")4
1.7686 0.7622
C(as.factor(pr), "contr.sum")5 C(as.factor(pr), "contr.sum")6
-1.2893 -2.1321
C(as.factor(pr), "contr.sum")7 C(as.factor(pr), "contr.sum")8
-3.2770 -4.1353
>
>
> matplot(y,type="l")
> matplot(matrix(fitted( lm(as.vector(y) ~ -1 + as.factor(cp) +
+ C(as.factor(pr),"contr.sum"))),ncol=6),type="l")
>
>
> library(preprocessCore)
> true.probes <- c(4,3,2,1,-1,-2,-3,-4)
>
> true.chips <- c(8,9,10,11,12,10)
>
> y <- outer(true.probes,true.chips,"+") + rnorm(48,0,0.25)
>
> y[8,4:6] <- c(11,12,10)+ 2.5 + rnorm(3,0,0.25)
> y[5,4:6] <- c(11,12,10)+-2.5 + rnorm(3,0,0.25)
>
>
>
> ###.C("plmd_fit_R", as.double(y), as.integer(8), as.integer(6),
> ### as.integer(2), as.integer(c(1,1,1,2,2,2) - 1),
> ### double(6 +2*8),
> ### double(48),
> ### double(48))
>
> ###matplot(matrix(.C("plmd_fit_R", as.double(y), as.integer(8), as.integer(6),
> ### as.integer(2), as.integer(c(1,1,1,2,2,2) - 1),
> ### double(6 +2*8),
> ### double(48),
> ### double(48))[[7]],ncol=6))
> ###
>
>
> ##.Call("R_plmd_model",y,0,1.3345,as.integer(c(1,1,1,2,2,2) - 1),as.integer(2))
> rcModelPLM(y)
$Estimates
[1] 8.3169023 9.2809276 10.4168042 11.2819293 12.1641491 10.3398623
[7] 3.6347409 2.5905522 1.7366863 0.7652716 -2.0548139 -2.3629422
[13] -3.3650431 -0.9444519
$Weights
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1.0000000 1.0000000 1.0000000 1.0000000 1.000000 1.0000000
[2,] 0.8730168 1.0000000 1.0000000 1.0000000 1.000000 1.0000000
[3,] 1.0000000 1.0000000 1.0000000 1.0000000 1.000000 1.0000000
[4,] 1.0000000 1.0000000 0.6963267 1.0000000 1.000000 1.0000000
[5,] 0.6975639 0.6158733 1.0000000 0.4873403 1.000000 0.6075995
[6,] 1.0000000 1.0000000 1.0000000 1.0000000 1.000000 1.0000000
[7,] 1.0000000 1.0000000 1.0000000 1.0000000 1.000000 1.0000000
[8,] 0.1099013 0.1329995 0.1308422 0.1290849 0.122359 0.1199671
$Residuals
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 0.14796819 0.03712246 -0.301050205 -0.05627386 0.199022368 -0.02678895
[2,] -0.46903675 0.10525974 0.333661640 0.34829319 -0.278989932 -0.09874766
[3,] 0.18612963 0.17264411 -0.009738675 -0.21170757 -0.011532707 -0.12579479
[4,] 0.19557522 -0.31820550 0.588089436 -0.26630757 0.006166356 -0.02673089
[5,] 0.58675206 0.66460741 0.323402537 -0.84047397 -0.322787098 -0.67417698
[6,] -0.05306488 -0.08392657 -0.154465963 -0.16482977 0.295285146 0.16100204
[7,] -0.06695385 0.08726733 -0.191844329 0.35095151 -0.296639940 0.11721928
[8,] -3.72583620 -3.07877488 -3.129474470 3.17210451 3.346512226 3.41319130
$StdErrors
[1] 0.2303741 0.2287995 0.2275630 0.2312391 0.2234607 0.2295250 0.2298418
[8] 0.2317326 0.2298418 0.2347044 0.2602041 0.2298418 0.2298418 0.5783378
$Scale
[1] 0.3044699
> rcModelPLMd(y,c(1,1,1,2,2,2))
$Estimates
[1] 8.005369 9.034160 10.082258 10.882428 11.899884 10.041785 3.945547
[8] 2.866283 2.045802 1.060794 -1.232278 -2.346678 -2.053827 -3.080891
[15] -3.891517 2.686766
$Weights
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1.0000000 1.0000000 0.9634219 1.0000000 1.0000000 1
[2,] 0.6166764 1.0000000 0.6808531 0.5660368 0.9197825 1
[3,] 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000 1
[4,] 1.0000000 0.7280995 0.4260700 1.0000000 1.0000000 1
[5,] 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000 1
[6,] 1.0000000 1.0000000 1.0000000 1.0000000 1.0000000 1
[7,] 1.0000000 1.0000000 1.0000000 0.5729866 0.8441314 1
[8,] 0.5718592 1.0000000 1.0000000 1.0000000 1.0000000 1
$Residuals
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 0.14869574 -0.02691561 -0.27731032 0.03242202 0.15248180 -0.03951711
[2,] -0.43323402 0.07629683 0.39247670 0.47206424 -0.29045532 -0.07640064
[3,] 0.18854776 0.11029661 0.01569179 -0.12132110 -0.05638269 -0.13683237
[4,] 0.21158582 -0.36696053 0.62711237 -0.16232864 -0.02509116 -0.02417600
[5,] 0.07574916 0.08883888 -0.16458803 -0.14910779 0.23334264 -0.08423484
[6,] -0.05064675 -0.14627407 -0.12903550 -0.07444331 0.25043516 0.14996446
[7,] -0.03957232 0.04988323 -0.14145046 0.46630137 -0.31652653 0.13114510
[8,] -0.46723777 0.11505792 0.15213630 -0.05961134 -0.02044007 0.08005141
$StdErrors
[1] 0.10057998 0.09251231 0.11350253 0.10277317 0.10277317 0.08542204
[7] 0.10225195 0.15766343 0.09331564 0.11594469 0.13431167 0.13407482
[13] 0.09331564 0.11433878 0.16442599 0.00000000
$WasSplit
[1] 0 0 0 0 1 0 0 1
>
> ###R_plmd_model(SEXP Y, SEXP PsiCode, SEXP PsiK, SEXP Groups, SEXP Ngroups)
>
>
>
>
>
> pr[seq(3,48,8)][1:3] <- 10
>
> y[seq(3,48,8)][1:3] <- c(8,9,10) -3 + rnorm(3,0,0.1)
> lm(as.vector(y) ~ -1 + as.factor(cp) + C(as.factor(pr),"contr.sum"))
Call:
lm(formula = as.vector(y) ~ -1 + as.factor(cp) + C(as.factor(pr),
"contr.sum"))
Coefficients:
as.factor(cp)1 as.factor(cp)2
7.907 8.955
as.factor(cp)3 as.factor(cp)4
10.072 10.713
as.factor(cp)5 as.factor(cp)6
11.670 9.821
C(as.factor(pr), "contr.sum")1 C(as.factor(pr), "contr.sum")2
4.078 3.024
C(as.factor(pr), "contr.sum")3 C(as.factor(pr), "contr.sum")4
2.147 1.239
C(as.factor(pr), "contr.sum")5 C(as.factor(pr), "contr.sum")6
-1.655 -1.919
C(as.factor(pr), "contr.sum")7 C(as.factor(pr), "contr.sum")8
-2.921 -3.896
C(as.factor(pr), "contr.sum")9
2.893
>
>
> proc.time()
user system elapsed
2.222 0.044 2.246
preprocessCore.Rcheck/tests/qnormtest.Rout
R Under development (unstable) (2023-03-12 r83975) -- "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(preprocessCore)
>
> err.tol <- 10^-8
>
> x <- matrix(c(100,15,200,250,110,16.5,220,275,120,18,240,300),ncol=3)
> x
[,1] [,2] [,3]
[1,] 100 110.0 120
[2,] 15 16.5 18
[3,] 200 220.0 240
[4,] 250 275.0 300
> normalize.quantiles(x)
[,1] [,2] [,3]
[1,] 110.0 110.0 110.0
[2,] 16.5 16.5 16.5
[3,] 220.0 220.0 220.0
[4,] 275.0 275.0 275.0
>
> x.norm.truth <- matrix(rep(c(110.0,16.5,220,275.0),3),ncol=3)
>
> if (all(abs(x.norm.truth - normalize.quantiles(x)) < err.tol) != TRUE){
+ stop("Disagreement in normalize.quantiles(x)")
+ }
>
> normalize.quantiles.determine.target(x)
[1] 16.5 110.0 220.0 275.0
>
> x.norm.target.truth <- c(16.5,110.0,220.0,275.0)
>
> if (all(abs(x.norm.target.truth - normalize.quantiles.determine.target(x)) < err.tol) != TRUE){
+ stop("Disagreement in normalize.quantiles.determine.target(x)")
+ }
>
>
> y <- x
> y[2,2] <- NA
> y
[,1] [,2] [,3]
[1,] 100 110 120
[2,] 15 NA 18
[3,] 200 220 240
[4,] 250 275 300
> normalize.quantiles(y)
[,1] [,2] [,3]
[1,] 134.44444 47.66667 134.44444
[2,] 47.66667 NA 47.66667
[3,] 226.11111 180.27778 226.11111
[4,] 275.00000 275.00000 275.00000
>
> y.norm.target.truth <- c(47.6666666666667,134.4444444444444,226.1111111111111,275.0000000000000)
>
> y.norm.truth <- matrix(c(134.4444444444444, 47.6666666666667, 134.4444444444444,
+ 47.6666666666667, NA, 47.6666666666667,
+ 226.1111111111111, 180.2777777777778, 226.1111111111111,
+ 275.0000000000000, 275.0000000000000, 275.0000000000000),byrow=TRUE,ncol=3)
>
>
> if (all(abs(y.norm.truth - normalize.quantiles(y)) < err.tol,na.rm=TRUE) != TRUE){
+ stop("Disagreement in normalize.quantiles(y)")
+ }
>
>
>
> if (all(abs(y.norm.target.truth - normalize.quantiles.determine.target(y)) < err.tol) != TRUE){
+ stop("Disagreement in normalize.quantiles.determine.target(y)")
+ }
>
>
>
> if (all(abs(normalize.quantiles.use.target(y,y.norm.target.truth) - y.norm.truth) < err.tol,na.rm=TRUE) != TRUE){
+ stop("Disagreement in normalize.quantiles.use.target(y)")
+ }
>
>
> x <- matrix(c(100,15,200,250,110,16.5,220,275,120,18,240,300),ncol=3)
> rownames(x) <- letters[1:4]
> colnames(x) <- LETTERS[1:3]
> y <- normalize.quantiles(x, keep.names = TRUE)
> if(!all(colnames(x)==colnames(y))){
+ stop("Disagreement between initial and final column names despite keep.names=TRUE")
+ }
> if(!all(rownames(x)==rownames(y))){
+ stop("Disagreement between initial and final row names despite keep.names=TRUE")
+ }
>
> proc.time()
user system elapsed
0.200 0.036 0.217
preprocessCore.Rcheck/preprocessCore-Ex.timings
| name | user | system | elapsed | |
| colSummarize | 0.002 | 0.000 | 0.004 | |
| normalize.quantiles.in.blocks | 0.043 | 0.001 | 0.042 | |
| rcModelPLMd | 0.018 | 0.000 | 0.019 | |
| rcModelPLMr | 0.037 | 0.004 | 0.040 | |
| rcModels | 0.004 | 0.003 | 0.008 | |
| subColSummarize | 0.005 | 0.001 | 0.005 | |
| subrcModels | 0.004 | 0.007 | 0.011 | |