Question: Summarize by gene for agilent microarray data - one color
gravatar for salamandra
5 months ago by
salamandra350 wrote:

I want to normalize agilent data (one color array).

Followed limma package user guide, but in the end got some genes repeated for different probes. How can I get data summarized by gene instead of probe? Is it ok, to just average probes of same gene?

Code so far:


# read SDRF file - equivalent to "targets" file in limma:
SDRF <- read.delim(SDRF.path,check.names=FALSE,stringsAsFactors=FALSE)

# read the gene expression intensity data:
x <- read.maimages(paste(dataDir, studyId, Data.Name, SDRF[,"Array Data File"], sep='/'), source="agilent", green.only=TRUE, other.columns="gIsWellAboveBG")

# gene annotation:
x$genes$EntrezID <- mapIds(RnAgilentDesign028282.db, x$genes$ProbeName, keytype="PROBEID", column="ENTREZID")
x$genes$Symbol <- mapIds(RnAgilentDesign028282.db, x$genes$ProbeName, keytype="PROBEID", column="SYMBOL")

# normexp background correction and quantile normalization:
y <- backgroundCorrect(x, method="normexp")
y <- normalizeBetweenArrays(y, method="quantile")

## filter probes:
Control <- y$genes$ControlType==1 
NoEntrez <-$genes$EntrezID) 
IsExpr <- rowSums(y$other$gIsWellAboveBG > 0) >= 4 
yfilt <- y[!Control & !NoEntrez & IsExpr, ]

# remove annotation columns not needed:
yfilt$genes <- yfilt$genes[,c("ProbeName","Symbol","EntrezID")]

# convert info to dataframe:
normData =
limma microarray R • 216 views
ADD COMMENTlink modified 5 months ago by Kevin Blighe69k • written 5 months ago by salamandra350
gravatar for Kevin Blighe
5 months ago by
Kevin Blighe69k
Republic of Ireland
Kevin Blighe69k wrote:

You can use avereps(), as I do here: A: How to process (seems) Agilent microarrry data?


ADD COMMENTlink written 5 months ago by Kevin Blighe69k
Please log in to add an answer.


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