This function will compress wiggle data by averaging values (including position) across non-overlapping windows. It will take any number of wiggle files, but can only handle one chromosome or chromosomal region per function call. This function assumes that all data is from the chromosome and ignores all chromosomal information. This function ignores all positions without reads mapping which may lead to slight variations in position depending on the wiggle data being used. Note: calls function compress. Written by Tovah Markowitz.
wiggle_compress(inputWiggles, window)
inputWiggles | Will accept two types of wiggle information:
|
---|---|
window | Size of window to be averaged across. Suggested size for chromosome-wide plots for presenting and publication: 200 bp. No default. |
A R (dplyr) data frame with averaged position and value for each wiggle. Column names are defined by input strings/variables.
# NOT RUN { wiggle_compress( red1[[6]], 200 ) wiggle_compress( c( "red1[[6]]", "hop1[[6]]", "smc4[[6]]" ), 200 ) # }