Question: Custom cdf files and genes
0
gravatar for 12021560-002
2.8 years ago by
12021560-0020 wrote:

Hello,

I am working on microarray data, as I am new to it that's why I need your help. I want to ask that if I download custom cdf files from brain array then how would I be able to extract genes specific to my disease of interest i.e t2d. Which packages would I use to do this analysis? Can someone guide me about this?

Your help will be really valued.

Thank You.

microarray cdf r3.2.2 • 436 views
ADD COMMENTlink modified 2.8 years ago by James W. MacDonald51k • written 2.8 years ago by 12021560-0020
Answer: Custom cdf files and genes
0
gravatar for James W. MacDonald
2.8 years ago by
United States
James W. MacDonald51k wrote:

The probeset IDs for the MBNI custom cdfs are the database ID concatenated to _st. So e.g., if you were using the Entrez Gene version, given an Entrez Gene ID of 1223, the probeset ID would be 1223_st.

If we assume you have run rma to get an ExpressionSet (that you called 'eset'), you can then make a featureData slot that includes all the annotation by doing something like (assuming this is human).

library(org.Hs.eg.db)
egids <- gsub("_st", "", featureNames(eset))
adf <- data.frame(PROBEID = featureNames(eset),
                         SYMBOL = mapIds(org.Hs.eg.db, egids, "SYMBOL", "ENTREZID"),
                         ENTREZID = egids,
                         GENENAME = mapIds(org.Hs.eg.db, egids, "GENENAME","ENTREZID"))

featureData(eset) <- AnnotatedDataFrame(data = adf)

Now if you process your data using limma, the annotation data will end up in the topTable output.

ADD COMMENTlink written 2.8 years ago by James W. MacDonald51k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 447 users visited in the last hour