Preprocessing Agilent Two-Color
I have an Agilent 2-Color dataset (GSE23611), i am trying to normalize it with Limma package. when I normalize dataset , hclust cannot distinguish between case and control after normalization. my code is :
library(limma)
files_case <- dir(pattern="*\\.txt$")
dat <- read.maimages(files_case,source="agilent")
dat2<-normalizeWithinArrays(dat, method="loess","normexp", offset=50)
dat2 <- normalizeBetweenArrays(dat2, method="Aquantile")
dat2 <- avereps(dat2, ID=dat2$genes$ProbeName)
dat.m<-dat2$M
rownames(dat.m)<-dat2$genes$ProbeName
dat.m=na.omit(dat.m)
dat.dist<-dist(t(dat.m))
plot(hclust(dat.dist))
is it true?
• 2,210 views
•
link
updated
by
Ram
•
written
by
babak1146 •
•
written
by
babak1146 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Agilent array in R
written by Monika •
-
error in processing Agilent two channel microarray raw data
written by firstorthopedicdoctor •
-
How to get gene level of differential expression from 2 color channel array with limma?
written by MatthewP
-
Summarize by gene for agilent microarray data - one color
written by salamandra
-
Agilent single color array
written by Leite
-
Two channel 4x44k array analysis in Limma
written by nemo566 •
-
Problem in Heat map generation of DE genes in Agilent Microarray data analysis
written by Vinay Singh •
-
Heatmap of Agilent two colors Dye Swap experiment
written by antgomo •
-
Creating Heatmaps Of Diffrentially Expressed Genes In Angilent Single-Color Array Analysis
written by Tonig
-
Normalize Red Channel Across Arrays Using Limma
written by David Quigley
Please use the formatting bar (especially the
codeoption) to present your post better. I've done it for you this time.Why were you expecting it to segregate cases and controls? The type of clustering that you are doing is unsupervised, i.e., using all probe expression values; thus, it should never be expected that this will segregate your dataset as you desire.
If you perform differential expression analysis between cases and controls and then perform supervised clustering using the statistically significant probes, then you should see segregation.