Hi!
I am doing analysis on my RNA-Seq data. When I reached to the GOSeq I encountered a problem as there is no support for the genome and gene references I need. The error message:
> pwf=nullp(genes,"mm10","refGene")
Error in getlength(names(DEgenes), genome, id) :
Length information for genome mm10 and gene ID refGene is not in the geneLenDataBase database. You will have to specify bias.data manually.
Is it even possible to manually add length data or how should I proceed? I could easily get the lengths from the reference files, but how can I import the gene lengths to GOSeq? Or maybe my only option is to wait for the upgrade of geneLenDataBase?
Thanks in advance!
I am doing analysis on my RNA-Seq data. When I reached to the GOSeq I encountered a problem as there is no support for the genome and gene references I need. The error message:
> pwf=nullp(genes,"mm10","refGene")
Error in getlength(names(DEgenes), genome, id) :
Length information for genome mm10 and gene ID refGene is not in the geneLenDataBase database. You will have to specify bias.data manually.
Is it even possible to manually add length data or how should I proceed? I could easily get the lengths from the reference files, but how can I import the gene lengths to GOSeq? Or maybe my only option is to wait for the upgrade of geneLenDataBase?
Thanks in advance!
Comment