reconsi package: vignette

Introduction

The aim of this pacakage is to improve simultaneous inference for correlated hypotheses using collapsed null distributions. These collapsed null distributions are estimated in an empirical Bayes framework through resampling. Wilcoxon rank sum test and two sample t-test are natively implemented, but any other test can be used.

Installation

library(BiocManager)
BiocManager::install("reconsi")
library(devtools)
install_github("CenterForStatistics-UGent/reconsi")
suppressPackageStartupMessages(library(reconsi))
cat("reconsi package version", as.character(packageVersion("reconsi")), "\n")
## reconsi package version 1.10.0

General use

We illustrate the general use of the package on a synthetic dataset. The default Wilcoxon rank-sum test is used.

The method provides an estimate of the proportion of true null hypothesis, which is close to the true 90%.

## [1] 0.8823949

The result of the procedure can be represented graphically as follows:

## Warning: Removed 2 rows containing missing values (geom_bar).
## Warning: It is deprecated to specify `guide = FALSE` to remove a guide. Please use `guide = "none"` instead.
## It is deprecated to specify `guide = FALSE` to remove a guide. Please use `guide = "none"` instead.

The approximate correlation matrix of test statistic in the univariate \(h\) distribution can be visualized as follows:

On the other hand, one can also view the full variance-covariance matrix of the test statistics as found by the resamples. Here we demonstrate how it can recover a block correlation structure in the data.

It is also possible to provide a custom test function, which must accept at least a ‘y’ response variable and a ‘x’ grouping factor. Additionally, quantile, distribution and density functions should be supplied for transformation through quantiles to z-values.

This framework also accepts more than 2 groups, and additional covariates through the “argList” argument.

If the null distribution of the test statistic is not known, it is also possbile to execute the procedure on the scale of the original test statistics, rather than z-values by setting zValues = FALSE. This may be numerically less stable.

Alternatively, the same resampling instances may be used to determine the marginal null distributions as to estimate the collapsed null distribution, by setting the “resamZvals” flag to true.

When no grouping variable is available, one can perform a bootstrap as resampling procedure. This is achieved by simply not supplying a grouping variable “x”. Here we perform a one sample Wilcoxon rank sum test for equality of the means to 0.

Case study

We illustrate the package using an application from microbiology. The species composition of a community of microorganisms can be determined through sequencing. However, this only yields compositional information, and knowledge of the population size can be acquired by cell counting through flow cytometry. Next, the obtained species compositions can multiplied by the total population size to yield approximate absolute cell counts per species. Evidently, this introduces strong correlation between the tests due to the common factor. In other words: random noise in the estimation of the total cell counts will affect all hypotheses. Therefore, we employ resampling to estimate the collapsed null distribution, that will account for this dependence.

The dataset used is taken from Vandeputte et al., 2017 (see manuscript), a study on gut microbiome in healthy and Crohn’s disease patients. The test looks for differences in absolute abundance between healthy and diseased patients. It relies on the phyloseq package, which is the preferred way to interact with our machinery for microbiome data.

The estimated tail-area false discovery rates can then simply be extracted as

##           0%          25%          50%          75%         100% 
## 1.179539e-16 4.536636e-03 3.724064e-01 1.000000e+00 1.000000e+00

Session info

Finally all info on R and package version is shown

sessionInfo()
## R version 4.2.1 (2022-06-23)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 20.04.5 LTS
## 
## Matrix products: default
## BLAS:   /home/biocbuild/bbs-3.16-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.16-bioc/R/lib/libRlapack.so
## 
## locale:
##  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
##  [3] LC_TIME=en_GB              LC_COLLATE=C              
##  [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       
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
## [1] reconsi_1.10.0
## 
## loaded via a namespace (and not attached):
##  [1] Biobase_2.58.0         sass_0.4.2             jsonlite_1.8.3        
##  [4] splines_4.2.1          foreach_1.5.2          bslib_0.4.0           
##  [7] assertthat_0.2.1       highr_0.9              stats4_4.2.1          
## [10] phyloseq_1.42.0        GenomeInfoDbData_1.2.9 yaml_2.3.6            
## [13] pillar_1.8.1           lattice_0.20-45        glue_1.6.2            
## [16] digest_0.6.30          XVector_0.38.0         colorspace_2.0-3      
## [19] htmltools_0.5.3        Matrix_1.5-1           plyr_1.8.7            
## [22] pkgconfig_2.0.3        zlibbioc_1.44.0        mvtnorm_1.1-3         
## [25] scales_1.2.1           pracma_2.4.2           tibble_3.1.8          
## [28] mgcv_1.8-41            generics_0.1.3         farver_2.1.1          
## [31] IRanges_2.32.0         ggplot2_3.3.6          cachem_1.0.6          
## [34] BiocGenerics_0.44.0    cli_3.4.1              survival_3.4-0        
## [37] magrittr_2.0.3         crayon_1.5.2           mclust_6.0.0          
## [40] evaluate_0.17          ks_1.13.5              fansi_1.0.3           
## [43] nlme_3.1-160           MASS_7.3-58.1          vegan_2.6-4           
## [46] tools_4.2.1            data.table_1.14.4      lifecycle_1.0.3       
## [49] matrixStats_0.62.0     stringr_1.4.1          Rhdf5lib_1.20.0       
## [52] S4Vectors_0.36.0       munsell_0.5.0          cluster_2.1.4         
## [55] Biostrings_2.66.0      ade4_1.7-19            compiler_4.2.1        
## [58] jquerylib_0.1.4        GenomeInfoDb_1.34.0    rlang_1.0.6           
## [61] rhdf5_2.42.0           grid_4.2.1             RCurl_1.98-1.9        
## [64] iterators_1.0.14       rhdf5filters_1.10.0    biomformat_1.26.0     
## [67] igraph_1.3.5           labeling_0.4.2         bitops_1.0-7          
## [70] rmarkdown_2.17         gtable_0.3.1           codetools_0.2-18      
## [73] multtest_2.54.0        DBI_1.1.3              reshape2_1.4.4        
## [76] R6_2.5.1               knitr_1.40             dplyr_1.0.10          
## [79] fastmap_1.1.0          utf8_1.2.2             KernSmooth_2.23-20    
## [82] permute_0.9-7          ape_5.6-2              stringi_1.7.8         
## [85] parallel_4.2.1         Rcpp_1.0.9             vctrs_0.5.0           
## [88] tidyselect_1.2.0       xfun_0.34