[BioC] FW: GenomicFeatures::makeTranscriptDbFromBiomart - Error in callSuper(...) : could not find function "initRefFields"

Cook, Malcolm MEC at stowers.org
Fri Feb 3 16:36:48 CET 2012


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.
> 
> 
> 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.
> 
> 
>    Marc
> 
> 
> 
> On 02/02/2012 01:58 PM, Cook, Malcolm wrote:
> > Hi Herve and Marc and Bioc fellow-travellers,
> >
> > I am unhappy to report that....
> >
> > It now happens that:
> >
> > 	>  txdb<-loadFeatures('dmelanogaster_gene_ensembl.sqlite')
> > 	Error in callSuper(...) : could not find function "initRefFields"
> >
> > ... just after I updated installed Bioc packages today, using
> >
> > 	source("http://bioconductor.org/biocLite.R")
> > 	old.packages(repos=biocinstallRepos())
> >
> > At first glance it seems to be the same as reported earlier here:
> >
> >
> 	http://permalink.gmane.org/gmane.science.biology.informatics.con
> ductor/38404
> >
> > However, trying the same workaround, namely, re-updating from source
> >
> > 	>  biocLite(c('GenomicFeatures'), type='source')
> >
> > and starting a new session, fails to fix... the error still occurs.
> >
> > As a potential workaround, I tried re-building the txdb from ensembl, but
> this is now failing (possibly due to new schema/data at Ensembl?), yielding
> the following error:
> >
> > 	txdb<-makeTranscriptDbFromBiomart(biomart="ensembl",
> dataset="dmelanogaster_gene_ensembl",circ_seqs=NULL)
> > 	Download and preprocess the 'transcripts' data frame ... OK
> > 	Download and preprocess the 'chrominfo' data frame ... OK
> > 	Download and preprocess the 'splicings' data frame ... Error in
> .extractCdsRangesFromBiomartTable(bm_table) :
> >    	BioMart data anomaly: for some transcripts, the cds cumulative
> length inferred from the exon and UTR info doesn't match the "cds_length"
> attribute from BioMart
> >
> > So, I looked for the hoped-for further workaround, of being able to specify
> an (archive) 'host' to makeTranscriptDbFromBiomart, (as
> discussed:https://stat.ethz.ch/pipermail/bioconductor/attachments/201111
> 08/a5ece78c/attachment.pl)_ which I would expect to call as
> >
> > 	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