Question: Summarize by gene for agilent microarray data - one color
0
gravatar for salamandra
5 months ago by
salamandra350
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:

library(limma)
library(tidyverse)
library(RnAgilentDesign028282.db)

# 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 <- is.na(y$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 = as.data.frame.EList(yfilt)
limma microarray R • 216 views
ADD COMMENTlink modified 5 months ago by Kevin Blighe69k • written 5 months ago by salamandra350
2
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?

Kevin

ADD COMMENTlink written 5 months ago by Kevin Blighe69k
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.3.0
Traffic: 2085 users visited in the last hour
_