Question: Gene Set Enrichment Analysis with topGO
gravatar for jacorvar
4.0 years ago by
European Union
jacorvar40 wrote:

Dear BioC community,

I have a boolean vector called pval where 1 means differentially expressed and 0 non-DEwhose element names are the Entrez IDs.

In order to make a GSEA according to GO MF, I do the next:

​GOobj <- new("topGOdata", description = "Simple session",
    ontology = 'MF', allGenes = pval,
    geneSel = function(x) return(x), nodeSize = 10,
    annot =, mapping = annotation(transcript), ID = 'entrez')
allGO = usedGO(object = GOobj)
resultFisher <- runTest(GOobj, algorithm = "classic", statistic = "fisher")
allRes <- GenTable(GOobj, classicFisher = resultFisher,
                   orderBy = "classicFisher", ranksOf = "classicFisher", topNodes = length(allGO))

My main problem is that the lowest p-values always correspond to very general terms, such as binding, and therefore they comprise lots of other (probably more interesting) terms. Is there any way to avoid this with the topGO package or is there another package that solves this issue?

topgo gene set analysis • 909 views
ADD COMMENTlink modified 4.0 years ago by James W. MacDonald51k • written 4.0 years ago by jacorvar40
Answer: Gene Set Enrichment Analysis with topGO
gravatar for James W. MacDonald
4.0 years ago by
United States
James W. MacDonald51k wrote:

The topGO package is based on methods that Adrian Alexa developed to do exactly what you are talking about. Paradoxically it seems that most of the examples in the vignette and the help pages describe the classic Fisher's exact instead. Anyway, you probably want either the elim, weight or weight01 methods (the latter being a combination of the two former methods). The weight01 is the default method, so if you just did

resultFisher <- runTest(GOobj, statistic = "fisher")

you would probably get what you wanted all along. If you care to know about the different methods, see the package vignette, which has references.

ADD COMMENTlink written 4.0 years ago by James W. MacDonald51k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 316 users visited in the last hour