topGO: how do I retrieve my genes contributing to the GO term enrichment
0
0
Entering edit mode
Priya • 0
@e7e6cc87
Last seen 19 months ago
Australia

I am not able to retrieve my genes contributing to the GO term being enriched.

Code should be placed in three backticks as shown below

#this is my code

allRes <- GenTable(hs_topGOdata, classicFisher = resultFisher,
 classicKS = resultFisher,
 orderBy = "classicFisher", ranksOf = "classicKS", topNodes = 20)

allRes
# retrieve genes2GO list from the "expanded" annotation in GOdata

allGO = genesInTerm(hs_topGOdata)

allGO["GO:0051241"]
SAM_ANOTATION <-  lapply(allGO,function(x) x[x %in% allGenes.hs] )

SAM_ANOTATION[["GO:0051241"]]

Everything runs fine and I get all the annotated genes which in my case for that term is 20 through allGO, however, I am unable to run the SAM_ANNOTATION function of lapply(), I am unsure if allGenes.hs is right. The allGene.hs is my named vector containing all of my genes and the p-values.

Please help

gene enrichment topGO GO retrieve • 631 views
ADD COMMENT

Login before adding your answer.

Traffic: 477 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