[BioC] read counts in sliding windows

Jason Lu jasonlu68 at gmail.com
Wed Aug 3 04:46:39 CEST 2011


Hi list,

I would like to seek your advice with this.
Here I have a GRangesList and a GappedAlignments objects (an alignment
to the genome). For each GRanges (with length 1) in the GRangesList, I
would like to do a sliding window on its ranges and count reads from
each window. So far I wasn't able to find an efficient way to do this.
My problem is similar to this one:
http://permalink.gmane.org/gmane.science.biology.informatics.conductor/34431.
But looping over the GRangesList (and using countOverlaps inside the
loop) takes significant time.

Thanks for any suggestions.

Best,

Jason



More information about the Bioconductor mailing list