Learn R Programming

celda (version 0.0.0.9000)

diffExpBetweenCellStates: Gene expression markers for cell clusters using MAST

Description

Finds markers (differentially expressed genes) for cell clusters using MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data (Finak et al, Genome Biology, 2015)

Usage

diffExpBetweenCellStates(counts, celda.mod, c1, c2 = NULL, only.pos = FALSE,
  log2fc.threshold = NULL, fdr.threshold = 1)

Arguments

counts

A numeric count matrix.

celda.mod

An object of class celda_C or celda_CG.

c1

Numeric. Cell cluster(s) to define markers for.

c2

Numeric. Second cell cluster(s) for comparison. If NULL (default) - use all other cells for comparison.

only.pos

Logical. Only return markers with positive log2fc (FALSE by default).

log2fc.threshold

Numeric. Only return markers whose the absolute values of log2fc are greater than this threshold (NULL by default).

fdr.threshold

Numeric. Only return markers whose false discovery rates (FDRs) are less than this threshold (1 by default).

Value

Data frame containing a ranked list (based on the absolute value of log2fc) of putative markers, and associated statistics (p-value, log2fc and FDR).