Calculate the average expression levels of each program (cluster) on single cell level, subtracted by the aggregated expression of control gene sets. All analyzed genes are binned based on averaged expression, and the control genes are randomly selected from each bin.
AddModuleScore(object, genes.list = NULL, genes.pool = NULL, n.bin = 25,
seed.use = 1, ctrl.size = 100, use.k = FALSE, enrich.name = "Cluster")
Seurat object
Gene expression programs in list
Number of bins of aggregate expression levels for all analyzed genes
Random seed for sampling
Number of control genes selected from the same bin per analyzed gene
Use gene clusters returned from DoKMeans()
Name for the expression programs
Returns a Seurat object with module scores added to object@meta.data
Tirosh et al, Science (2016)