Error topGO printGenes
0
0
Entering edit mode
@irvingjgalo23-14252
Last seen 6.8 years ago

Hi, I want to extract the enriched genes from the significative GO terms with the printGenes function.

goID <- "GO:0005730"

gt <- printGenes(GOdata, whichTerms = goID,

                 chip = "org.At.tair.db", numChar = 40)

But when I run the script I get the following error:

Error in data.frame(LL.id = as.integer(unlist(mget(all.genes, envir = LL.lib,  : 
  arguments imply differing number of rows: 255, 444

According to that node I have 243 of 255 genes annotated

I hope you can help me

Regards

 

topGo • 1.5k views
ADD COMMENT
0
Entering edit mode

Could you provide how did you create your GOdata? This way the example is reproducible and I can check how to solve it. Thanks

ADD REPLY
0
Entering edit mode

Thanks for your time, that's how I built my Godata

GOdata  <- new ("topGOdata", ontology = "CC",
               allGenes = all_pvlue , geneSel = topDiffGenes,
               description = "Arabidopsis_Genes", nodeSize = 10,
               annot = annFUN.org, mapping = "org.At.tair.db")

 

 

Regards

ADD REPLY
0
Entering edit mode

You don't pass a function to geneSel, do you?

Could you paste the output of print(GOdata)?

ADD REPLY
0
Entering edit mode

Hi, im having exactly the same problem. did you manage to solve it ?

ADD REPLY

Login before adding your answer.

Traffic: 1040 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6