Dear List,
I find this behavior of the coverage function quite inconvenient:
> myRanges = IRanges(start=c(10,20,30,100), width=5)
> myCov = coverage(myRanges, start=50, end=80)
> myCov
31-integer "XRleInteger" instance
[1] Error in x...@values[group, drop = TRUE] : subscript out of bounds
When I apply the coverage function, I do not care whether there are is
at least one sequence within the specified region. I think, it would be
more sensible to return an XRleInteger object with 31 "0" in the
situation above.
Best wishes,
Hans-Ulrich
_______________________________________________
Bioc-sig-sequencing mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/bioc-sig-sequencing
