getMarkerGenes.rnaseq: Marker Gene Detection

Description Usage Arguments Details Value Author(s) Examples

View source: R/getMarkerGenes.rnaseq.R

Description

Function to detect marker genes using normalized RNA-seq data

Usage

1
getMarkerGenes.rnaseq(data.mat, class.vec=colnames(data.mat), samples2compare="all", annotate=FALSE, gene.ids.type="ensembl", score.cutoff=1)

Arguments

data.mat

RNA-seq gene expression matrix with genes corresponding to rows and samples corresponding to columns.

class.vec

A character vector containing the classes of samples (columns) of data.mat in the same order as provided in the matrix.

samples2compare

A character vector with the sample names to be compared (e.g. c("liver", "lung", "brain")). By default all samples in the reference matrix are used.

annotate

A boolean value. If TRUE the gene symbol and the entrez gene id are shown.

gene.ids.type

Type of the used gene identifiers, the following gene identifiers are supported: ensembl, refseq and ucsc gene ids. default is ensembl.

score.cutoff

A value in the interval [0,1] to filter the markers according to the specificity score. The default value is 1 (no filtering).

Details

For each marker in the output list, the gene id and the corresponding score are shown. If annotate is TRUE, the gene symbol and the entrez gene id are shown. The score is used to rank the markers according to their specificity. A lower value means a higher specificity.

Value

A list with marker genes associated with each sample type.

Author(s)

Khadija El Amrani <a.khadija@gmx.de>

Examples

1
2
3
4
5
data(ref.mat)
res.list <- getMarkerGenes.rnaseq(ref.mat, class.vec = colnames(ref.mat), samples2compare="all", annotate=TRUE, gene.ids.type="ensembl", score.cutoff=1)
names(res.list)
## show the first 20 markers of liver
res.list[["liver_markers"]][1:20]

MGFR documentation built on Nov. 8, 2020, 11:11 p.m.