[BioC] FW: GenomicFeatures::makeTranscriptDbFromBiomart - Error in callSuper(...) : could not find function "initRefFields"
Cook, Malcolm
MEC at stowers.org
Fri Feb 3 21:08:37 CET 2012
Hi Marc,
> A host argument has indeed been added to the devel branch.
Great
> But more
> pertinent to your current issue is the fact that I am able to easily
> load your database using loadFeatures() or loadDB(). So I don't see
> anything wrong with the DB on this end. We tested it in both Linux and
> Mac and are unable to recreate the problem.
Hmmm...
> This makes me think that there is something wrong with your current
> setup. Your version of R could be a little more current, I don't really
> anticipate that this would cause a problem like this, but it might be
> worth updating to 2.14.1 anyhow. And the only package I have that is
> more current than yours is XML, which should not matter here. Are you
> sure that you have read permissions for your copy of this DB?
Yes. I am sure.
FWIW: google finds a few others similarly afflicted , but only in very the recent past:
5 days ago - http://lists.hexdump.org/pipermail/twitter-users-hexdump.org/2012-January/000027.html
5 days ago - same problem in ggplot- https://groups.google.com/forum/?fromgroups#!topic/ggplot2/FciMhMfDZp4
But wait. SUCCESS:
I just tried rebuilding from source, and, after the following I was able to load the database
biocLite(c('RCurl','rtracklayer','GenomicFeatures','AnnotationDbi','Biobase','GenomicRanges','IRanges'), type='source')
HOORAY!
I guess we're not going to find the root cause this time.
Like last time: http://permalink.gmane.org/gmane.science.biology.informatics.conductor/38404
That's OK.
If you have any suggestions for me to try something, fire away.
Other than that....
Thanks for your help sleuthing with me. !!!
Cheers,
Malcolm
>
>
> Marc
>
>
> On 02/03/2012 07:43 AM, Cook, Malcolm wrote:
> > Hi Marc,
> >
> >> I cannot say why your older annotation database is not working, but I
> >> suspect that it is just out of phase with recent changes introduced when
> >> you updated the code. But the error you reported from when you tried
> to
> >> generate a new TxDb package is an unhappy message indeed. It means
> that
> >> the data for cds ranges and total cds length (both from biomaRt) no
> >> longer agree with each other. In other words, the data from the current
> >> drosophila ranges in biomaRt seems to disagree with itself, and so the
> >> code is refusing to make a package out of this data as a result.
> >>
> >> So to try and get the 1st issue fixed I need more information. I
> >> suspect that what has happened to your old DB is that is is missing some
> >> metadata that is required by the newer versions of GenomicFeatures. Is
> >> there a place where you could put your DB so that I could DL it and have
> >> a look at it?
> > Sure - great - thanks - here it is:
> >
> > ftp://ftp.stowers-
> institute.org/pub/MEC/dmelanogaster_gene_ensembl.sqlite
> >
> >> To get the 2nd issue fixed probably involves talking to ensembl about
> >> their CDS data for fly to see if we can resolve the discrepancy.
> > I would be happy to take this to them.
> >
> > Can you recommend a best way to get a more diagnostic trace from the
> attempt at txdb creation so we can correctly report to ensembl team the
> errant transcript(s) ?
> >
> > Finally, do you have any information on the progress of adding a 'host'
> option to makeTranscriptDbFromBiomart ?
> > If I install Bioc-dev should I expect to find this implemented?
> > And if I do, will a txdb so-created be backward compatible with current
> released GenomicFeatures (i.e. should I be able to load it)?
> >
> > Thanks so much for you help,
> >
> > Cheers,
> >
> > Malcolm
> >
> >
> >
> >>> txdb<-makeTranscriptDbFromBiomart(host='
> >> jun2011.archive.ensembl.org/biomart/martservice',biomart="ensembl",
> >> dataset="dmelanogaster_gene_ensembl",circ_seqs=NULL)
> >>> but I found that the `host` option is not (yet?) supported (in released).
> >>>
> >>> So, most urgently to me, is to find a workaround to my primary problem,
> >> namely, of getting past the missing 'initRefFields' when calling
> `loadFeatures`
> >> .
> >>> Is there a fix/patch/workaround that can address this?
> >>>
> >>> Thanks,
> >>>
> >>> Malcolm Cook
> >>>
> >>>
> >>>
> >>>> sessionInfo()
> >>> R version 2.14.0 (2011-10-31)
> >>> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
> >>>
> >>> locale:
> >>> [1] C
> >>>
> >>> attached base packages:
> >>> [1] stats graphics grDevices utils datasets methods base
> >>>
> >>> other attached packages:
> >>> [1] GenomicFeatures_1.6.7 AnnotationDbi_1.16.11 Biobase_2.14.0
> >>> [4] GenomicRanges_1.6.6 IRanges_1.12.5
> >>>
> >>> loaded via a namespace (and not attached):
> >>> [1] BSgenome_1.22.0 Biostrings_2.22.0 DBI_0.2-5 RCurl_1.9-5
> >>> [5] RSQLite_0.11.1 XML_3.9-2 biomaRt_2.10.0
> rtracklayer_1.14.4
> >>> [9] tools_2.14.0 zlibbioc_1.0.0
> >>> _______________________________________________
> >>> 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.conduct
> >> or
> >>
> >> _______________________________________________
> >> 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