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:

  1. wiggle_scatter will extract the necessary information using wiggle_compress for all chromosomes

  2. wiggle_scatter_plot will make a simple .eps file of the results.

Written by Tovah Markowitz

wiggle_scatter(wiggleData1, wiggleData2, window = 5000)

Arguments

wiggleData1

A list of 16 chr wiggle data (output of readall_tab). No default.

wiggleData2

A second set of wiggle data (output of readall_tab). No default.

window

Parameter of wiggle_scatter. Size of window to be compressed. Default is 5000 bp.

Value

Output is a large (dplyr) data frame with compressed data from all chromosomes. See wiggle_compress for more information.

Examples

# NOT RUN {
a <- wiggle_scatter(red1, hop1, 4000)

wiggle_scatter_plot(a)
# }