Search
Question: Help with annotating eset (Affy Arrays)
0
gravatar for rumroses12
2.0 years ago by
rumroses120
United States
rumroses120 wrote:

I'm trying to add annotation to an RMA normalized spreadsheet of microarrray expression data (probe ID in the first column, sample names in the first row). The arrays are from Affymetrix.

Combing through the forums it looks like I need to use the AnnotationDbi package. I've looked at the package documentation, but can't figure out how to use it. I'm new to R and programming.

So far, I have succeeded in transforming my data into an Expressionset object:

 

source("http://bioconductor.org/biocLite.R")
biocLite(c("Biobase"))
library("Biobase")
exprsFile <- "C:/filepath/file.txt"
exprs <- as.matrix(read.table(exprsFile, header = TRUE, sep="\t", row.names=1, as.is=TRUE))
minimalSet <- ExpressionSet(assayData=exprs)

 

Any thoughts on how I can add annotation to this using AnnotationDbi or some other package?

ADD COMMENTlink modified 2.0 years ago • written 2.0 years ago by rumroses120
0
gravatar for James W. MacDonald
2.0 years ago by
United States
James W. MacDonald45k wrote:

You have to know the type of array, and from that divine the annotation package name. Usually the annotation package name is based on taking the array name, going to all lower case, and stripping off all the non-letters. So as an example, the HG-U133_plus2 array's annotation package is hgu133plus2.db. It is a bit more complex for the more recent arrays, as they can be summarized at different levels. As an example, the HuGene 1.0 ST array has two annotation packages, the hugene10sttranscriptcluster.db and hugene10sttprobeset.db package, reflecting the summarization at the transcript and probeset levels, respectively.

Once you know the annotation package, you simply install, and use select() to get what you want. For most arrays this will result in a one-to-many mapping, as lots of annotations are not unique. You will have to decide how you want to deal with that complexity. One possibility is to add the annotations one at a time, using mapIds(), which by default will just take the first of any multi-hits.

ADD COMMENTlink written 2.0 years ago by James W. MacDonald45k
0
gravatar for rumroses12
2.0 years ago by
rumroses120
United States
rumroses120 wrote:

Thanks- I appreciate the help!

ADD COMMENTlink written 2.0 years ago by rumroses120
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 2.2.0
Traffic: 285 users visited in the last hour