[BioC] IRanges::Rle::[ arguably should return NA when i is out of range
Cook, Malcolm
MEC at stowers.org
Mon Mar 10 19:52:22 CET 2014
I think the semantics of indexing into an Rles should be the same as indexing into the corresponding vector.
But it is not in the case when indices are out of bounds.
Example:
library(IRanges)
> c(1:10)[10:13]
[1] 10 NA NA NA
> Rle(c(1:10))[10:13]
Error in normalizeSingleBracketSubscript(i, x) :
subscript contains NAs or out of bounds indices
I've made similar arguments before which were taken up, as in https://stat.ethz.ch/pipermail/bioconductor/2013-September/054820.html
Is there a good reason NOT to change these semantics?
In the meantime, any suggested clever workarounds to get such extraction semantics?
Thanks!
~ Malcolm Cook
Computational Biology / Shilatifard Lab - Stowers Institute for Medical Research - Kansas City
More information about the Bioconductor
mailing list