# NOT RUN {
# 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)
# }
Run the code above in your browser using DataLab