Finds markers (differentially expressed genes) for each of the identity classes in a dataset
FindAllMarkers(object, genes.use = NULL, logfc.threshold = 0.25,
test.use = "wilcox", min.pct = 0.1, min.diff.pct = -Inf,
print.bar = TRUE, only.pos = FALSE, max.cells.per.ident = Inf,
return.thresh = 0.01, do.print = FALSE, random.seed = 1,
min.cells.gene = 3, min.cells.group = 3, latent.vars = NULL,
assay.type = "RNA", ...)Seurat object
Genes to test. Default is to use all genes
Limit testing to genes which show, on average, at least X-fold difference (log-scale) between the two groups of cells. Default is 0.25 Increasing logfc.threshold speeds up the function, but can miss weaker signals.
Denotes which test to use. Available options are:
"wilcox" : Wilcoxon rank sum test (default)
"bimod" : Likelihood-ratio test for single cell gene expression, (McDavid et al., Bioinformatics, 2013)
"roc" : Standard AUC classifier
"t" : Student's t-test
"tobit" : Tobit-test for differential gene expression (Trapnell et al., Nature Biotech, 2014)
"poisson" : Likelihood ratio test assuming an underlying poisson distribution. Use only for UMI-based datasets
"negbinom" : Likelihood ratio test assuming an underlying negative binomial distribution. Use only for UMI-based datasets
"MAST : GLM-framework that treates cellular detection rate as a covariate (Finak et al, Genome Biology, 2015)
"DESeq2 : DE based on a model using the negative binomial distribution (Love et al, Genome Biology, 2014)
only test genes that are detected in a minimum fraction of min.pct cells in either of the two populations. Meant to speed up the function by not testing genes that are very infrequently expressed. Default is 0.1
only test genes that show a minimum difference in the fraction of detection between the two groups. Set to -Inf by default
Print a progress bar once expression testing begins (uses pbapply to do this)
Only return positive markers (FALSE by default)
Down sample each identity class to a max number. Default is no downsampling.
Only return markers that have a p-value < return.thresh, or a power > return.thresh (if the test is ROC)
FALSE by default. If TRUE, outputs updates on progress.
Random seed for downsampling
Minimum number of cells expressing the gene in at least one of the two groups, currently only used for poisson and negative binomial tests
Minimum number of cells in one of the groups
Remove the effects of these variables, used only when test.use is one of
'negbinom', 'poisson', or 'MAST'
Type of assay to perform DE for (default is RNA)
Additional parameters to pass to specific DE functions
Matrix containing a ranked list of putative markers, and associated statistics (p-values, ROC score, etc.)
# NOT RUN {
all_markers <- FindAllMarkers(object = pbmc_small)
head(x = all_markers)
# }
Run the code above in your browser using DataCamp Workspace