BioMart Attributes/Filters missing
1
0
Entering edit mode
terka • 0
@terka-9304
Last seen 5.8 years ago
European Union

Hi, 

I would like to use getBM function to get ENSEMBL gene identifiers to ProbeNames from Agilent Whole Mouse Genome Microarray Kit (design ID 014868), 4x44k. However, in the available filters, I can only find Agilent WholeGenome 4x44k v1 probe ID(s) and Agilent WholeGenome 4x44k v2 probe ID(s). Do you have any suggestions about how can I find the IDs to the specific microarray that I used?

 

Thanks a lot in advance, 

 

Teresa

microarray agilent biomart annotation • 866 views
ADD COMMENT
1
Entering edit mode
@james-w-macdonald-5106
Last seen 9 hours ago
United States

There are at least two other ways you could approach this. First, if you are using limma to do the analysis (and you probably should be doing that), then when you read the data into an RGList, the 'genes' list item will contain accession numbers that you can use to map to Ensembl using the org.Mm.eg.db package. Something like

library(limma)
library(org.Mm.eg.db)
dat <- read.maimages(<filenames go here>, source = "agilent.median")
gns <- dat$genes

## there should be a column called 'Systematic name' or something like that, which contains the accession numbers. You can then map using mapIds():

ensgene <- mapIds(org.Mm.eg.db, as.character(gns[,<correct column number goes here>]), "ENSEMBL","ACCNUM")

An alternative would be to get the data from GEO. You could use GEOquery to get those data, or just download by hand and parse out the column containing the Ensembl IDs.

ADD COMMENT
0
Entering edit mode

Great, this helps a lot! Using mapIds did it really well. Thank you, James!

ADD REPLY

Login before adding your answer.

Traffic: 241 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6