Calculate gene set scores
calculate_gs_scores.Rd
Calculate gene set scores using the approach described in Tirosh2016
Usage
calculate_gs_scores(
scMuffinList = NULL,
gs_list = NULL,
mc.cores = 1,
nbins = 25,
nmark_min = 5,
ncells_min = 10,
k = 100,
kmin = 50,
score_type = c("relative", "mean"),
verbose = FALSE,
na.rm = TRUE,
overwrite = FALSE
)
Arguments
- scMuffinList
scMuffinList object
- gs_list
list of gene sets
- mc.cores
number of cores
- nbins
number of bins to split the distribution of average gene expression
- nmark_min
number of minimum markers that are required for the succesful calculation of a gene set score
- ncells_min
number of minimum cells in which a gene set has to be succesfully calculated
- k
number of permutations
- kmin
minimum number of permutations; due to missing values it is hard to ensure that a gene set score can be compared to k permutations in every cell
- score_type
type of score. if "relative", than the score is the difference between the observed gene set average expression and that of a k permutations; if "mean" the score is equal to the observed gene set average expression
- verbose
verbosity
- na.rm
whether to use NA or not
- overwrite
whether to update or not gene_set_scoring and gene_set_scoring_full elements of scMuffinList.