DEsingle: DEsingle: Detecting differentially expressed genes from...

Description Usage Arguments Value Author(s) See Also Examples

View source: R/DEsingle.R

Description

This function is used to detect differentially expressed genes between two specified groups of cells in a raw read counts matrix of single-cell RNA-seq (scRNA-seq) data. It takes a non-negative integer matrix of scRNA-seq raw read counts or a SingleCellExperiment object as input. So users should map the reads (obtained from sequencing libraries of the samples) to the corresponding genome and count the reads mapped to each gene according to the gene annotation to get the raw read counts matrix in advance.

Usage

1
DEsingle(counts, group, parallel = FALSE, BPPARAM = bpparam())

Arguments

counts

A non-negative integer matrix of scRNA-seq raw read counts or a SingleCellExperiment object which contains the read counts matrix. The rows of the matrix are genes and columns are samples/cells.

group

A vector of factor which specifies the two groups to be compared, corresponding to the columns in the counts matrix.

parallel

If FALSE (default), no parallel computation is used; if TRUE, parallel computation using BiocParallel, with argument BPPARAM.

BPPARAM

An optional parameter object passed internally to bplapply when parallel=TRUE. If not specified, bpparam() (default) will be used.

Value

A data frame containing the differential expression (DE) analysis results, rows are genes and columns contain the following items:

Author(s)

Zhun Miao.

See Also

DEtype, for the classification of differentially expressed genes found by DEsingle.

TestData, a test dataset for DEsingle.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
# Load test data for DEsingle
data(TestData)

# Specifying the two groups to be compared
# The sample number in group 1 and group 2 is 50 and 100 respectively
group <- factor(c(rep(1,50), rep(2,100)))

# Detecting the differentially expressed genes
results <- DEsingle(counts = counts, group = group)

# Dividing the differentially expressed genes into 3 categories
results.classified <- DEtype(results = results, threshold = 0.05)

miaozhun/DEsingle documentation built on June 21, 2018, 3:51 a.m.