Issues with .extractEnsemblReleaseFromDbVersion() and get_organism_from_Ensembl_Mart_dataset()
1
0
Entering edit mode
@mattchambers42-10186
Last seen 6.8 years ago

I found a couple issues testing this with older Ensembl archives:

1. .extractEnsemblReleaseFromDbVersion() should support Ensembl 54 (May 2009), which doesn't have "genes" in its version string. Something like this should work:

.extractEnsemblReleaseFromDbVersion <- function (db_version) 
{
    db_version <- tolower(db_version)
    sub("^ensembl(?: genes)? ([0-9]+).*$", "\\1", db_version, perl=TRUE)
}

2. get_organism_from_Ensembl_Mart_dataset() makes the reasonable but flawed assumption that Ensembl will only add, not remove, organisms from its current_mysql directory. Agambiae was removed from the main ensembl site when it was moved to the metazoa site, but the older assembly is still on Ensembl 54. In other words, listDatasets(may2009) says agambiae is available, so get_organism_from_Ensembl_Mart_dataset() should just work. :) The fix should simply be to pass along the Ensembl release to .Ensembl_getMySQLCoreDir().

genomicfeatures • 587 views
ADD COMMENT
0
Entering edit mode
@herve-pages-1542
Last seen 8 hours ago
Seattle, WA, United States

Hi Matt,

Thanks for the report. I made those changes in GenomicFeatures 1.28.3 (release) and 1.29.6 (devel).

H.

ADD COMMENT

Login before adding your answer.

Traffic: 504 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6