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:
Post a Comment