plotErrorRate {transcriptR} | R Documentation |
A simple helper function that plot results of
estimateGapDistance
function call.
plotErrorRate(object, color = c("#1B9E77", "#D95F02", "#7570B3"), xlab = "Gap distance (kb)", ylab = "Error rate", ...) ## S4 method for signature 'TranscriptionDataSet' plotErrorRate(object, color = c("#1B9E77", "#D95F02", "#7570B3"), xlab = "Gap distance (kb)", ylab = "Error rate", ...)
object |
A |
color |
|
xlab |
|
ylab |
|
... |
Further arguments passed to plot. |
The tested gap distances are plotted on the x-axis and corresponding
error rates on the y-axis. Three curved lines depict the two error
types calculated by estimateGapDistance
and the sum of
both errors. The vertical dashed line depicts the gap distance with
the smallest sum of two errors.
plot
Armen R. Karapetyan
### Load TranscriptionDataSet object data(tds) ### Load reference annotations (knownGene from UCSC) data(annot) ### Estimate gap distance minimazing error rate ### Define the range of gap distances to test gdr <- seq(from = 0, to = 10000, by = 1000) estimateGapDistance(object = tds, annot = annot, coverage.cutoff = 5, filter.annot = FALSE, gap.dist.range = gdr) plotErrorRate(object = tds, lwd = 2)