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.
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!