Entering edit mode
4.7 years ago
trueflyingsheep
•
0
Is there an easy API to merge individual databases and the GO term database? The following seems clumsy.
library(AnnotationDbi)
library(org.Mm.eg.db)
library(GO.db)
library(dplyr)
select(org.Mm.eg.db, gene_symbols, 'GO', 'SYMBOL') %>%
left_join(select(GO.db, .$GO, 'TERM', 'GOID'), c(GO = 'GOID'))