Remove centromeric regions from wiggle data

This function allows you to remove the region containing the centromere in each chromosome. The original motivation was to analyse Red1 ChIP-seq data in rec8 mutants (show a strong enrichment around the centromeres).

remove_centromeres(wiggleData, regionSize = 50000)

Arguments

wiggleData	As a list of the 16 chr wiggle data (output of readall_tab). No default.
regionSize	Number indicating the size (in bp) of the region to remove (centered on the centromere of each chromosome). Defaults to 50'000 bp.

Value

A similar list of the 16 chr wiggle data after removal of the rows for the centromeric regions (both positions and signal) in each chromosome data frame.

Examples

# NOT RUN {
remove_centromeres(rec8)

remove_centromeres(rec8, regionSize = 30000)
# }