Sample mapping code

sample mapping code

library(biomaRt)
ensembl=useMart("ensembl")
#show possible datasets
listDatasets(ensembl)

#your results are listed here
dataset <- read.csv('Hits_PvsM_HighinMimic_filt.csv')
#you have several different ensembl results
unique(substr(dataset$Gene, 1,7))

#get a mapping file to get externam gene names from stickleback
stickleback_table <- getBM(attributes=c('ensembl_gene_id','external_gene_name','entrezgene'),
                      filters='ensembl_gene_id',
                      values=dataset$Gene,
                      mart=useDataset("gaculeatus_gene_ensembl",mart=ensembl))
#repeat this with other species in your results set by changing the species in the use dataset.  
#To find the abbreviation for that species use listDatasets(ensembl)
mapped_table <- rbind(stickleback_table) #add other tables

write.csv(mapped_table, file="Mapped Genes.csv")