GO Enrichment Analysis in R - Count all the parents of each Term
0
0
Entering edit mode
2.7 years ago
rrapopor ▴ 40

Hi,

I'm trying to implement a new GO tool in R. After I have a list of genes, I'm using the following code to retrieve the relevant annotation for each gene:

library(AnnotationDbi)

AnnotationDbi::select(org.Hs.eg.db, keys = genes.list, 
                                     columns = c("ENTREZID", "GO"), keytype = id.type)

Once I compared the results to a published GO tool, it seems that the count for each gene is separate of its parents - meaning that each node in the org.Hs.eg.db database has a unique set of genes, which is not part of its parents genes. If I understand correctly, a gene in one of the nodes should also appear in all of his ancestors.

The problem is when I'm trying to add all the ancestors of a term for each gene, in some cases where the same node has more than 1 parent (example below), I'm counting the ancestor more than once. GOrilla GO terms DAG

Is there any tool / database in r where each gene has it's relevant GO terms all the way up in the tree, and the problem I mentioned above is fixed?

Thank you!

Annotations Analysis Enrichment R GO • 490 views
ADD COMMENT

Login before adding your answer.

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