[BioC] Problem with intraIDMapper

Marc Carlson mcarlson at fhcrc.org
Mon Apr 23 19:42:33 CEST 2012


Hi Robert,

This method needs to be overhauled and updated now that we have made 
some other improvements to our infrastructure.  A bug has been added to 
my issue tracker for future development and also: thanks for reporting 
this.  It's also good to know that people are using this function.  I 
was actually wondering if maybe I should deprecate it, but now I think 
that maybe I should update it instead.

In the meantime I think you can get much the same functionality by using 
the new select methods.  There is now even a select method for 
inparanoid in the latest development version of AnnotationDbi.  The 
select methods handle missing IDs much more gracefully than the much 
older idConverter function, and they also allow you to recover multiple 
things at once.

In the case you listed below you can already do this by using select 
like so;

library(org.Mm.eg.db)

keys = c("NR_003620","NM_007376")
cols= c("SYMBOL","GENENAME") ## can add however many things you want here.
keytype="REFSEQ"
select(org.Mm.eg.db,keys=keys, cols=cols, keytype=keytype)



   Marc


On 04/16/2012 10:21 AM, Faryabi, Robert (NIH/NCI) [F] wrote:
> Hello,
>
> Here is a problem with idMapper wrappers in AnnotationDbi.
> If you try to map in id that is not in the data base, the function breaks instead of returning a flag that allows handling the situation.
>
>> idConverter("NR_003620", srcSpecies="MUSMU", destSpecies="MUSMU", srcIDType="REFSEQ", destIDType="SYMBOL", keepMultGeneMatches=FALSE, keepMultProtMatches=FALSE, keepMultDestIDMatches = TRUE)
> Error in names(ans)<- unlist(make.name.tree(x, recursive, what.names),  :
>    attempt to set an attribute on NULL
>
>
>> sessionInfo()
> R version 2.14.1 (2011-12-22)
> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
>
> locale:
> [1] C
>
> attached base packages:
> [1] grid      stats     graphics  grDevices utils     datasets  methods
> [8] base
>
> other attached packages:
>   [1] ggplot2_0.9.0        RSvgDevice_0.6.4.1   Vennerable_2.1
>   [4] xtable_1.7-0         reshape_0.8.4        plyr_1.7.1
>   [7] lattice_0.20-0       RBGL_1.30.1          graph_1.32.0
> [10] venneuler_1.1-0      rJava_0.9-3          pheatmap_0.6.0
> [13] heatmap.plus_1.3     RColorBrewer_1.0-5   Heatplus_2.1.0
> [16] gplots_2.10.1        KernSmooth_2.23-7    caTools_1.12
> [19] bitops_1.0-4.1       gdata_2.8.2          gtools_2.6.2
> [22] org.Mm.eg.db_2.6.4   RSQLite_0.11.1       DBI_0.2-5
> [25] AnnotationDbi_1.16.9 Biobase_2.14.0       rj_1.0.3-7
>
> loaded via a namespace (and not attached):
>   [1] IRanges_1.12.4   MASS_7.3-17      colorspace_1.1-1 dichromat_1.2-4
>   [5] digest_0.5.1     memoise_0.1      munsell_0.3      proto_0.3-9.2
>   [9] reshape2_1.2.1   rj.gd_1.0.3-3    scales_0.2.0     stringr_0.6
> [13] tools_2.14.1
>
> Thanks,
> Robert
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list