[Bioc-devel] GenomicRanges: Seqinfo indexing with non-existing key
tooyoung at gmail.com
Thu Feb 27 14:00:05 CET 2014
To me it sounds like what I am expecting, from a R base point of view.
On Thu, 2014-02-27 at 13:49 +0100, Julian Gehring wrote:
> In the current bioc-stable and bioc-devel, a 'Seqinfo' object can be
> indexed successfully with any character key, even it is non-existing:
> #+BEGIN_SRC R
> si = seqinfo(BSgenome.Hsapiens.UCSC.hg19)
> ind2 = "999"
> ind2 %in% seqnames(si) ## FALSE, no seqlevel with this name
> si[ind2] ## works, but should not
> Shouldn't this throw an error?
> Best wishes
> Bioc-devel at r-project.org mailing list
More information about the Bioc-devel