[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