coverageplot {chipseq} | R Documentation |
A function that plots one or two coverage vectors over a relatively small interval in the genome.
coverageplot(peaks1, peaks2 = NULL, i = 1, xlab = "Position", ylab = "Coverage", opposite = TRUE, ...)
peaks1, peaks2 |
A set of peaks as described by ranges over a coverage vector. |
i |
Which peak to use. |
xlab, ylab |
Axis labels. |
opposite |
Logical specifying whether the two peaks should be plotted on opposite sides (appropriate for positive and negative strand peaks). |
... |
extra arguments. |
Deepayan Sarkar
cov <- Rle(c(1:10, seq(10, 1, -2), seq(1,5,2), 4:1), rep(1:2, 11)) peaks <- slice(cov, 3) peaks.cov <- Views(cov, ranges(peaks)) peaks.cov.rev <- rev(peaks.cov) coverageplot(peaks.cov, peaks.cov.rev, ylab = "Example")