[BioC] Proper use of annotation package

Naxerova, Kamila naxerova at fas.harvard.edu
Mon Mar 18 21:46:10 CET 2013

Hi all,

I have a little problem using an annotation package that I know how to work around, but I am wondering how to do it more elegantly and efficiently.

I am analyzing a bunch of Mouse Gene 2.0 ST arrays. I built my own annotation package (with much help from all of you!). I am using Limma and want to look up annotation for diff exp genes provided by topTable(). So it's really a standard situation. On the Bioconductor website, this sequence of commands is suggested (http://www.bioconductor.org/help/workflows/annotation-data/)

tbl <- topTable(efit, coef=2)
ids <- tbl[["ID"]]
entrez <- hgu95av2ENTREZID[ids]

Looks beautiful! But when I try to do the same thing, I get:

ids <- tbl[["ID"]]

Error in .checkKeys(value, Lkeys(x), x at ifnotfound) : 
  value for "17549282" not found

This error is evidently generated because some of the ids don't map to any accession numbers. I can work around this by filtering my ids first, but am I doing it wrong? Of course lots of probe ids on the array are not going to map to any accession numbers or symbols or names -- why can't they just come back with NA instead of an error message and abortion of the whole process?


More information about the Bioconductor mailing list