I am currently working with a list of GO terms linked to gene ID's.
For completeness like so, but then for 100s of genes:
GENE1 GO:xxxx GO:xxxxx GO:xxxxx
Because I want to group my genes into 'biological-logical' groups like: development, digestion, metabolition etc. I was thinking to get the relevant GO term to do so: the GO term just below 'biological process', e.g.: GO:0008152.
Is there a way with the GO.db package to do so? I have been trying the following:
> get("GO:0006508", GOBPANCESTOR)
But than I get all ancestor terms, is there a way to get only the '2nd from the top'?
I guess I am not the first person who want to group genes based on very broad biological processes, I do not want to reinvent the wheel so if there is another -better- way, please let me know!
Thanks in advance!