This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Plotting Gene Ontology analysis results and questions about Fold Enrichment and other parameters

Hello everyone, I've performed a Gene Ontology analysis on a list of DEGs, and I obtained, as a result, a list of GO terms. For each term, it's shown how many times it's found in up and down-regulated genes. Two p-values, representing statistical significant for such findings, are given too.

I've used the function goana from limma in R Bioconductor.

Now I'm looking for a tool to create Gene Ontology plots.

To date, I've only created with Excel the classic 'bar plots', like these:

enter image description here

However, I've noticed there are more complex types of plots, like the ones where GO terms are represented by circles, whose size is proportional to their up/down regulation. Here are a couple of examples.

enter image description here

enter image description here

These plots refer to parameters not appearing in the results of my previous analysis, including a Fold enrichment, FDR and a z-score.

Is there any way to manually calculate Fold enrichment? Any suggestion in general about how to plot GO results starting from just the GO terms up and down counts and their p-values?

Thanks in advance for your help!

ontology enrichment gene go

1 answer

In case anyone is coming to this thread here is a great package to generate this kind of data https://wencke.github.io/. These kinds of plots require transforming your data in combination with expression data to generate Z score and such.

Log in to answer this question.