[Bioc-sig-seq] IRanges - coverage

[Hans-Ulrich Klein]

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 at 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