[BioC] "islands" and number of reads in the island in package chipseq/ShortRead

Nora Rieber n.rieber at dkfz-heidelberg.de
Tue Dec 1 13:17:25 CET 2009


Dear all,

I'm looking for a fast way of computing how many reads map into a
certain location (e.g. between position 300 and 400 on chromosome 14).

Is there any function to compute this in the packages chipseq or
ShortRead? After reading the vignette for chipseq, I thought
viewSums(islands) computes the number of reads, but it seems to compute
the sum of the coverage.

I've written my own code to do this but it is pretty slow and I wonder
if there's a better way.

Thanks,
Nora


More information about the Bioconductor mailing list