Parsing and sorting of uncorrected read and sequence information files
Loads WIG files for readcount, GC, and mappability data for non-overlapping windows of fixed length (i.e. bins), and returns a structure ready to used for readcount correction. See Details for specifics about file assumptions.
Pathname to WIG file containing readcounts per bin.
Pathname to WIG file containing GC content per bin.
Pathname to WIG file containing average mappability per bin.
Set to TRUE if messages are desired
The number of lines in the three input files are expected to be identical, although the order and names of chromsomes in the file need not be identical. Chromosome lengths are required to be identical and unique, and if the latter is not true, the order of the chromosomes must then be identical.
At present, these three WIG files are expected to be generated by external programs, namely those from the HMMcopy suite (see See Also), rather than by existing R packages out of space and memory considerations when working with high coverage full genome samples.
RangedData object, where each row
entry represents a bin, with the three values from the input as columns named
reads, gc, and map.
- correctedReadcount Suite
correctReadcount, to correct the readcounts in the resultant
1 2 3 4 5
Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.