In this document, we show that the Wilcoxon-Mann-Whitney test is comparable or superior to alternative methods.
Two alternative methods could be compared with the Wilcoxon-Mann-Whitney (WMW) test proposed by BioQC: the Kolmogorov-Smirnov (KS) test, and the Student’s t-test, or more particularly, the Welch’s test which does not assume equal sample number or equal variance, which is appropriate in the setting of gene expression studies.
Based on these considerations, BioQC implements a computationally efficient version of the WMW test. In order not to confuse end-users, no alternative methods are implemented.
Nevertheless, in order to demonstrate the power of WMW test in comparison with the KS-test or the t-test, we performed the sensitivity benchmark described in the simulation studies, for the two alternative tests respectively.
Figure 1: Sensitivity benchmark. Expression levels of genes in the ovary signature are dedicately sampled randomly from normal distributions with different mean values. The lines show the enrichment score for the Wilcoxon-Mann-Whitney test, the t-test and the Kolmogorov-Smirnov test respectively. In the right panel, outliers were added by adding a random value to 1% of the simulated genes.
As expected, the results suggest, that both the KS-test and the WMW-test are robust to noise, while the performance of the t-test drops significantly on noisy data. Additionally, the WMW-test appears to be superior to the KS-test for low expression differences.
Since the KS-test is so slow, we did not replicate the sensitivity benchmark from the simulation studies using the enrichment score rank. While it takes BioQC about 4 seconds on a single thread to test all 155 signatures, it already takes the KS-test about 2 seconds to test a single signature.
## test replications elapsed relative
## 2 runKS() 5 11.088 1.000
## 1 runWMW() 5 17.830 1.608## R version 4.4.0 beta (2024-04-15 r86425)
## Platform: x86_64-pc-linux-gnu
## Running under: Ubuntu 22.04.4 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.19-bioc/R/lib/libRblas.so
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: America/New_York
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] ggplot2_3.5.1 plyr_1.8.9 reshape2_1.4.4
## [4] hgu133plus2.db_3.13.0 rbenchmark_1.0.0 gplots_3.1.3.1
## [7] gridExtra_2.3 latticeExtra_0.6-30 lattice_0.22-6
## [10] org.Hs.eg.db_3.19.1 AnnotationDbi_1.66.0 IRanges_2.38.0
## [13] S4Vectors_0.42.0 BioQC_1.32.0 Biobase_2.64.0
## [16] BiocGenerics_0.50.0 testthat_3.2.1.1 knitr_1.46
##
## loaded via a namespace (and not attached):
## [1] tidyselect_1.2.1 farver_2.1.1 dplyr_1.1.4
## [4] blob_1.2.4 Biostrings_2.72.0 bitops_1.0-7
## [7] fastmap_1.1.1 digest_0.6.35 lifecycle_1.0.4
## [10] statmod_1.5.0 KEGGREST_1.44.0 RSQLite_2.3.6
## [13] magrittr_2.0.3 compiler_4.4.0 rlang_1.1.3
## [16] sass_0.4.9 tools_4.4.0 utf8_1.2.4
## [19] yaml_2.3.8 labeling_0.4.3 interp_1.1-6
## [22] bit_4.0.5 RColorBrewer_1.1-3 pkgload_1.3.4
## [25] KernSmooth_2.23-22 withr_3.0.0 desc_1.4.3
## [28] grid_4.4.0 fansi_1.0.6 caTools_1.18.2
## [31] colorspace_2.1-0 edgeR_4.2.0 scales_1.3.0
## [34] gtools_3.9.5 cli_3.6.2 rmarkdown_2.26
## [37] crayon_1.5.2 generics_0.1.3 httr_1.4.7
## [40] DBI_1.2.2 cachem_1.0.8 stringr_1.5.1
## [43] zlibbioc_1.50.0 XVector_0.44.0 vctrs_0.6.5
## [46] jsonlite_1.8.8 bit64_4.0.5 jpeg_0.1-10
## [49] locfit_1.5-9.9 limma_3.60.0 jquerylib_0.1.4
## [52] glue_1.7.0 stringi_1.8.3 gtable_0.3.5
## [55] GenomeInfoDb_1.40.0 deldir_2.0-4 UCSC.utils_1.0.0
## [58] munsell_0.5.1 tibble_3.2.1 pillar_1.9.0
## [61] htmltools_0.5.8.1 brio_1.1.5 GenomeInfoDbData_1.2.12
## [64] R6_2.5.1 rprojroot_2.0.4 evaluate_0.23
## [67] highr_0.10 png_0.1-8 memoise_2.0.1
## [70] bslib_0.7.0 Rcpp_1.0.12 xfun_0.43
## [73] pkgconfig_2.0.3