I am trying to use the
getGeneLengthAndGCContent function from the
EDASeq library to retrieve gene lengths for c. elegans. I retrieved Ensemble gene ids from Biomart for c.elegans and I'm using those as input to getGeneLengthAndGCContent (I picked "ensemble_gene_id" from the filters of the worm annotation Mart). Here's an example of how I'm trying to use the function:
getGeneLengthAndGCContent("WBGene00001042", org, mode=c("biomart", "org.db"))
When I try to run this, I get
Error in getGeneLengthAndGCContent("WBGene00001042", org, mode = c("biomart", : Only ENTREZ or ENSEMBL gene IDs are supported.