Get GO Term Size from GO ID
1
0
Entering edit mode
21 months ago
cthangav ▴ 110

I have a matrix of GO Term ID's and their GO term names. I'd like to order them by Term size, from largest to smallest number of genes.

So far I have been pulling all the genes from org.Mm.eg.db and counting the number of rows, but this is a slow process and difficult to (l)apply to many GO terms.

nrow(AnnotationDbi::select(org.Mm.eg.db, keytype="GOALL", keys=c("GO:0072089"), columns="SYMBOL"))

Is there a simple way to get the GO size of a list of Mm.db GO Terms so I can arrange them by size?

GO-term R • 860 views
ADD COMMENT
2
Entering edit mode
21 months ago

Try something along the line of the following:

library(GO.db)
library("org.Mm.eg.db")
annotations <- mapIds(org.Mm.eg.db, keys(org.Mm.eg.db, "GO"), "ENSEMBL", "GO", multiVals = "list")

This gives list of all GO term IDs with a vector of all Ensembl genes associated with each term. columns(org.Mm.eg.db) will give you a list of available identifiers you could use instead of Ensembl genes.

ADD COMMENT

Login before adding your answer.

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