I don't know the answer to your question, but please be informed that on the AnnotationHub such OrgDb for tomato is available. Maybe useful to know. I don't work with tomato, so I don't know how extensive the annotation information really is (AFAIK these are only based on the info available at the NCBI), and whether it is comparable to you OrgDb.
> library(AnnotationHub)
> hub = AnnotationHub()
snapshotDate(): 2016-01-25
> query(hub, c("Solanum"))
AnnotationHub with 4 records
# snapshotDate(): 2016-01-25
# $dataprovider: Inparanoid8, ftp://ftp.ncbi.nlm.nih.gov/gene/DATA/
# $species: Solanum lycopersicum, Solanum tuberosum
# $rdataclass: Inparanoid8Db, OrgDb
# additional mcols(): taxonomyid, genome, description, tags, sourceurl,
# sourcetype
# retrieve records with, e.g., 'object[["AH10593"]]'
title
AH10593 | hom.Solanum_lycopersicum.inp8.sqlite
AH10606 | hom.Solanum_tuberosum.inp8.sqlite
AH48047 | org.Solanum_lycopersicum.eg.sqlite
AH48064 | org.Solanum_tuberosum.eg.sqlite
>
> org.Sl.eg.db <- hub[["AH48047"]]
>
> columns(org.Sl.eg.db) #which annotation data can be retrieved?
> keytypes(org.Sl.eg.db) #which identifiers can be queried with?
What information (bimaps?) is missing from the OrgDb in AnnotationHub? How are you creating the OrgDb package (using code from AnnotationForge)? Providing a code example of the information you're trying to extract and how you are creating the new package will make it easier for others to help.
Also, are you aware of the select() interface for interacting with the OrgDbs? If not, See ?AnnotationDbi::select.
Thank you for your kindness about select() interface.
It seems there is no AnnotationDb object in my package. (but there are AnnDbBimaps)
> require(org.Slycopersicum.eg.db)
> columns(org.Slycopersicum.eg.db)
Error in columns(org.Slycopersicum.eg.db) :
error in evaluating the argument 'x' in selecting a method for function 'columns': Error: object 'org.Slycopersicum.eg.db' not found
> x <- org.Slycopersicum.egITAG
> mapped_genes <- mappedkeys(x)
> mapped_genes
Thank you for your kindness about select() interface.
It seems there is no AnnotationDb object in my package. (but there are AnnDbBimaps)
> require(org.Slycopersicum.eg.db)
> columns(org.Slycopersicum.eg.db)
Error in columns(org.Slycopersicum.eg.db) :
error in evaluating the argument 'x' in selecting a method for function 'columns': Error: object 'org.Slycopersicum.eg.db' not found
> x <- org.Slycopersicum.egITAG
> mapped_genes <- mappedkeys(x)
> mapped_genes
Thank you for your kindness about select() interface.
It seems there is no AnnotationDb object in my package. (but there are AnnDbBimaps)
> require(org.Slycopersicum.eg.db)
> columns(org.Slycopersicum.eg.db)
Error in columns(org.Slycopersicum.eg.db) :
error in evaluating the argument 'x' in selecting a method for function 'columns': Error: object 'org.Slycopersicum.eg.db' not found
> x <- org.Slycopersicum.egITAG
> mapped_genes <- mappedkeys(x)
> mapped_genes
I'm fairly new to the annotations so I had to seek counsel from Marc on this one (thanks Marc!).
The fact that your original OrgDb did not have an AnnotationDb object indicates you had an old version of AnnotationForge. Current release is 1.12.2 and devel is 1.13.8. It looks like you updated at some point and the second package did contain an AnnotationDb. When you post a question on the support site please include the output of sessionInfo() so we can see what versions you're using.
The bimaps have been replaced with the select() interface. The makeOrgPackage() was not intended to support bimaps but instead it creates an OrgDb from any set of IDs you choose. (Old OrgDb packages with bimaps were created with hard-coded with templates etc.) The important things to check in an OrgDb made with makeOrgPackage() is that an AnnotationDbi object is present and the select() interface returns the fields originally defined. MAPCOUNTS is not important or relevant here.
What information (bimaps?) is missing from the OrgDb in AnnotationHub? How are you creating the OrgDb package (using code from AnnotationForge)? Providing a code example of the information you're trying to extract and how you are creating the new package will make it easier for others to help.
Also, are you aware of the select() interface for interacting with the OrgDbs? If not, See ?AnnotationDbi::select.
Valerie
The missing info is like this.
I created the sqlite file with https://github.com/kozo2/org.Slycopersicum.eg.db-gen
, and create the OrgDb with https://github.com/kozo2/org.Slycopersicum.eg.db
I'm using AnnotationForge in https://github.com/kozo2/org.Slycopersicum.eg.db-gen/blob/master/tomatoMakeOrgPkg.R
This is just for the sqlite not for the package itself.
A code example is what I already mentioned.
I just want to check the MAPCOUNTS and know why this line does not work.
https://github.com/kozo2/org.Slycopersicum.eg.db/blob/master/R/objects.R#L119
Thank you for your kindness about select() interface.
It seems there is no AnnotationDb object in my package. (but there are AnnDbBimaps)
mapped_genes output is too long so I ommited.
The missing info is like this.
I created the sqlite file with https://github.com/kozo2/org.Slycopersicum.eg.db-gen
, and create the OrgDb with https://github.com/kozo2/org.Slycopersicum.eg.db
I'm using AnnotationForge in https://github.com/kozo2/org.Slycopersicum.eg.db-gen/blob/master/tomatoMakeOrgPkg.R
This is just for the sqlite not for the package itself.
A code example is what I already mentioned.
I just want to check the MAPCOUNTS and know why this line does not work.
https://github.com/kozo2/org.Slycopersicum.eg.db/blob/master/R/objects.R#L119
Thank you for your kindness about select() interface.
It seems there is no AnnotationDb object in my package. (but there are AnnDbBimaps)
mapped_genes output is too long so I ommited.
The missing info is like this.
I created the sqlite file with https://github.com/kozo2/org.Slycopersicum.eg.db-gen
, and create the OrgDb with https://github.com/kozo2/org.Slycopersicum.eg.db
I'm using AnnotationForge in https://github.com/kozo2/org.Slycopersicum.eg.db-gen/blob/master/tomatoMakeOrgPkg.R
This is just for the sqlite not for the package itself.
A code example is what I already mentioned.
I just want to check the MAPCOUNTS and know why this line does not work.
https://github.com/kozo2/org.Slycopersicum.eg.db/blob/master/R/objects.R#L119
Thank you for your kindness about select() interface.
It seems there is no AnnotationDb object in my package. (but there are AnnDbBimaps)
mapped_genes output is too long so I ommited.