This function is designed to work with function wiggle_scatter_plot
to
allow for an overall generalized comparison of two data sets.
Great for comparing replicates, checking controls, etc.
The two functions are used as follows:
wiggle_scatter
will extract the necessary information using
wiggle_compress
for all chromosomes
wiggle_scatter_plot
will make a simple .eps file of the results.
Written by Tovah Markowitz
wiggle_scatter(wiggleData1, wiggleData2, window = 5000)
wiggleData1 | A list of 16 chr wiggle data (output of |
---|---|
wiggleData2 | A second set of wiggle data (output of |
window | Parameter of |
Output is a large (dplyr) data frame with
compressed data from all chromosomes. See wiggle_compress
for more
information.
# NOT RUN { a <- wiggle_scatter(red1, hop1, 4000) wiggle_scatter_plot(a) # }