plotExpression: Plot the expression signal of individual genes

In Biopeak: Identification of Impulse-Like Gene Expression Changes in Short Genomic Series Data

Description

This function plots the expression signal of a defined gene and marks the main peak location with a dashed line.

Usage

plotExpression(exprmat, gene, series, peakdet)

Arguments

exprmat	A numeric matrix with expression series data with variables as rownames.
gene	A character string (not case-sensitive) defining the gene to be plotted.
series	A numeric vector defining the experimental series (e.g. time-points of sample acquisition).
peakdet	A list returned by the peakDetection function.

Value

This function does not return any value but generates a plot.

Author(s)

David Lauenstein

Examples

# Example based on the heat-shock dataset
data(heat)
heat = as.matrix(heat)
# Define series
series <- c(37,40,41,42,43)
# Run the peak detection algorithm
peakdet <- peakDetection(heat, series, type ='rnaseq', actstrength = 1.5,
prominence = 1.3, minexpr = 5000)
# Plot the expression signal of the gene CXCL5
plotExpression(heat, 'CXCL5', series, peakdet)

Biopeak documentation built on Aug. 21, 2019, 5:10 p.m.