[Bioc-devel] GenomicRanges: Seqinfo indexing with non-existing key
Julian Gehring
julian.gehring at embl.de
Thu Feb 27 13:49:55 CET 2014
Hi,
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
library(GenomicRanges)
library(BSgenome.Hsapiens.UCSC.hg19)
si = seqinfo(BSgenome.Hsapiens.UCSC.hg19)
ind2 = "999"
ind2 %in% seqnames(si) ## FALSE, no seqlevel with this name
si[ind2] ## works, but should not
#+END_SRC
Shouldn't this throw an error?
Best wishes
Julian
More information about the Bioc-devel
mailing list