= R Annotation Packages from Bioconductor = === Example from Zebrafish === Install package {{{ source("http://bioconductor.org/biocLite.R") biocLite("org.Dr.eg.db") # load package library(org.Dr.eg.db) }}} Now let's list the contents to see what's available {{{ > ls("package:org.Dr.eg.db") [1] "org.Dr.eg" "org.Dr.eg.db" [3] "org.Dr.eg_dbconn" "org.Dr.eg_dbfile" [5] "org.Dr.eg_dbInfo" "org.Dr.eg_dbschema" [7] "org.Dr.egACCNUM" "org.Dr.egACCNUM2EG" [9] "org.Dr.egALIAS2EG" "org.Dr.egCHR" [11] "org.Dr.egCHRLENGTHS" "org.Dr.egCHRLOC" [13] "org.Dr.egCHRLOCEND" "org.Dr.egENSEMBL" [15] "org.Dr.egENSEMBL2EG" "org.Dr.egENSEMBLPROT" [17] "org.Dr.egENSEMBLPROT2EG" "org.Dr.egENSEMBLTRANS" [19] "org.Dr.egENSEMBLTRANS2EG" "org.Dr.egENZYME" [21] "org.Dr.egENZYME2EG" "org.Dr.egGENENAME" [23] "org.Dr.egGO" "org.Dr.egGO2ALLEGS" [25] "org.Dr.egGO2EG" "org.Dr.egMAPCOUNTS" [27] "org.Dr.egORGANISM" "org.Dr.egPATH" [29] "org.Dr.egPATH2EG" "org.Dr.egPFAM" [31] "org.Dr.egPMID" "org.Dr.egPMID2EG" }}} I want to convert a between ensembl gene ids and Entrez gene ids {{{ # Get the entrez gene identifiers that are mapped to an ACCNUM mapped_genes <- mappedkeys(org.Dr.egACCNUM) # Convert to a list eg2acc <- as.list(org.Dr.egACCNUM[mapped_genes]) }}}