I am a newbie at gene ontology analysis in R. I have a list of Uniprot ID from different organisms. These ids were found as a result of a blast search. I want to find GO term hierarchy of these proteins for Biological process category. Then I want to sum up all categories from leaves to root with number of proteins in corresponding category. In the end I want to have something like:
GO:0008150 (7 proteins)
-> Child 1 (7 proteins)
-> Child 1.1 (5 proteins) -> Child 1.2 (2 proteins)
I found that GO.db library returns ancestors, children. However, they are not ordered by hierarchy. Uniprot IDs can be queried via UniProt.ws package. But I haven't tried it yet.
If you can help me, I would be appreciated.