Batch Correction Tab
Single Cell Toolkit
David Jenkins
The Section of Computational Biomedicine, Boston University School of Medicine, Boston, MA;Program in Bioinformatics, Boston University, Boston, MAdfj@bu.edu
Mohammed Muzamil Khan
The Section of Computational Biomedicine, Boston University School of Medicine, Boston, MA;Program in Bioinformatics, Boston University, Boston, MAmmkhan@bu.edu
W. Evan Johnson
The Section of Computational Biomedicine, Boston University School of Medicine, Boston, MA;Program in Bioinformatics, Boston University, Boston, MAwej@bu.edu Source:
vignettes/v06-tab04_Batch-Correction.Rmd
v06-tab04_Batch-Correction.Rmd
Introduction
Batch correction is available as a subtab in the QC & Filtering tab. Because of the complexities of the library preparation and the low starting material in scRNA-Seq experiments, non-biological variation (batch effects) are present and can be a major source of variation present in single cell experiments (Hicks et al. 2017). ComBat is a widely used method for adjusting for batch effects in microarray and RNA-Seq data (Johnson, Li, and Rabinovic 2007). If users identify variation associated with a technical effect, ComBat can be run within the SCTK to remove this variation before further downstream analysis. Users can choose an annotation present in the annotation data frame and add additional covariates to the ComBat model before performing batch correction. After batch correction, the ComBat results are stored as an additional assay in the SCExperiment object, which can then be used in the other analysis tabs within the SCTK.
Plot Batch Effect
Analysis on the batch correction tab is performed on the assay selected in the “Select Assay” field. To visualize the batch effect present in the data, select an annotation column from the annotation data frame in the “Select Batch Annotation” drop down and an experimental condition annotation in the “Select Condition Annotation” drop down. A set of boxplots will appear that show the percent variation explained by condition+batch, condition alone, and batch alone.
Run Batch Correction
Select a batch correction method from the “Select Method” drop-down. Currently, only ComBat is supported, but additional methods will be added in later versions of the toolkit. After batch correction, the corrected data will be saved as an additional assay in the SCtkExperiment object. Choose a name for this assay in the “Assay Name to Use” field.
ComBat
To run ComBat batch correction, select a batch annotation, add any additional covariates to the model, and adjust any of the ComBat parameters available. For details about the available options for ComBat analysis, see the ComBat documentation.
Session info
## R version 4.1.2 (2021-11-01)
## Platform: x86_64-apple-darwin17.0 (64-bit)
## Running under: macOS Big Sur 10.16
##
## Matrix products: default
## BLAS: /Library/Frameworks/R.framework/Versions/4.1/Resources/lib/libRblas.0.dylib
## LAPACK: /Library/Frameworks/R.framework/Versions/4.1/Resources/lib/libRlapack.dylib
##
## locale:
## [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] BiocStyle_2.22.0
##
## loaded via a namespace (and not attached):
## [1] rstudioapi_0.13 knitr_1.38 magrittr_2.0.3
## [4] R6_2.5.1 ragg_1.2.2 rlang_1.0.2
## [7] fastmap_1.1.0 stringr_1.4.0 tools_4.1.2
## [10] xfun_0.30 cli_3.2.0 jquerylib_0.1.4
## [13] systemfonts_1.0.4 htmltools_0.5.2 yaml_2.3.5
## [16] digest_0.6.29 rprojroot_2.0.3 pkgdown_2.0.2
## [19] bookdown_0.25 textshaping_0.3.6 BiocManager_1.30.16
## [22] purrr_0.3.4 sass_0.4.1 fs_1.5.2
## [25] memoise_2.0.1 cachem_1.0.6 evaluate_0.15
## [28] rmarkdown_2.14 stringi_1.7.6 compiler_4.1.2
## [31] bslib_0.3.1 desc_1.4.1 jsonlite_1.8.0