The purpose of this pipeline is to streamline the process for analyzing RNA-seq data with potential batch effects.
The pipeline includes:
- quantile normalization
- log-transformation of counts
- ComBat (location) batch correction
- voom calculation of weights.
The functions in this package can be grouped into two main categories:
- The functions used for assessing batch effects.
a.
makeSVD
b.pcRes
c.plotPC
- The functions for removing batch effect and computing weights for limma.
a.
qNorm
b.log2CPM
c.voomMod
d.combatMod
e.batchSEQ
batchSEQ
is the pipeline function. It combines qNorm
, log2CPM
,
voomMod
, and combatMod
into one step.
To begin, install Bioconductor along with a few dependencies that cbcbSEQ uses:
source("http://bioconductor.org/biocLite.R")
biocLite(c('limma', 'preprocessCore', 'sva', 'pasilla'))
Next, use devtools to install the latest version of cbcbSEQ from Github:
require(devtools)
install_github("cbcbSEQ", user="kokrah")
If all went well you should now be able to load cbcbSEQ:
require(cbcbSEQ)
vignette('cbcbSEQIntro', package='cbcbSEQ')