This function allows you to plot average chromosome signal for one or two selected samples.
It takes as input the output of chr_coverage
(16x2 R data frame).
chr_coverage_plot(coverageDataA, coverageDataB, genome, meanNorm = FALSE, yMax, onScreen = TRUE, fileName, colorA = "grey50", colorB = "green")
coverageDataA | A 16x2 data frame of coverage: chromosome and average signal. No default. |
---|---|
coverageDataB | Optional 16x2 data frame of coverage: chromosome and average signal. No default. |
genome | A string representing the genome used for mapping. No default. |
meanNorm | Boolean indicating whether input data are normalized to genome-wide
averages ( |
yMax | Optional number to be used as the max Y scale value in the plots. No default. |
onScreen | Boolean indicating plots should be returned to the screen ( |
fileName | A string to name the output .pdf file in case ( |
colorA | Optional R color for sample A. Defaults to |
colorB | Optional R color for sample B. Defaults to |
A dot plot of one or two samples, either on screen or as a .pdf file (in the working directory).
# NOT RUN { chr_coverage_plot(WT, rec8, genome = 'SK1', onScreen = TRUE, colorB = 'red') chr_coverage_plot(WT, dot1, genome = 'S288C', meanNorm = TRUE, onScreen = FALSE, fileName='chr_coverage_WT_and_dot1.pdf') # }