Thursday, April 26, 2012

Convert orthologous homologous genes

http://www.ncbi.nlm.nih.gov/homologene



#http://bogdan.org.ua/2010/10/27/batch-retrieve-entrezgene-homologs-using-ncbi-homologene-and-r.html
#source('http://bioconductor.org/biocLite.R')
#biocLite('annotationTools')
library(annotationTools)
homologene <- read.delim('../inputs/homologene-build-51.1.data', header=FALSE)

# use BioMart
# human (species ID 9606)
# mouse (species ID 10090)
mouse.entrez <- read.csv(sep='\t', 'mouse.entrez.ids.csv')

# mouse to human
mouse.ortholog <- getHOMOLOG(unlist(mouse.entrez$EntrezGene.ID), 9606, homologene)

unlist(mouse.ortholog)

No comments: