[BioC] GenomicFeatures::makeTranscriptDbFromBiomart - Error in callSuper(...) : could not find function "initRefFields"
Cook, Malcolm
MEC at stowers.org
Thu Feb 2 22:43:07 CET 2012
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.conductor/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/20111108/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
>
More information about the Bioconductor
mailing list