[BioC] Get Gene annotation (GO terms)
Marc Carlson
mcarlson at fhcrc.org
Thu Mar 17 22:16:14 CET 2011
What species? Those look like human gene symbosl, so I will give an
example for that:
library(org.Hs.eg.db)
sym <- c("ACSL1","ADAMTSL1","ATG2B","BACH2")
## Then you can just merge together the GO to entrez gene mappings with
the SYMBOL to entrez gene mappings like this:
merge(toTable(org.Hs.egSYMBOL2EG[sym]), toTable(org.Hs.egGO), by.x=
"gene_id", by.y="gene_id")
Does that make sense?
Marc
On 03/17/2011 02:11 PM, David martin wrote:
> Hi,
> Which package should i use to convert a list of gene names
> (ACSL1,ADAMTSL1,ATG2B,BACH2...) into Go id. The idea behind is to get
> a list of significant GO terms based on the input genes).thanks for
> pointing me to right direction.
>
> thanks,
> david
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list