Question: Preprocessing Agilent Two-Color
gravatar for babak1146
3.1 years ago by
babak114620 wrote:

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 :

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)

is it true?

gene R genome • 1.2k views
ADD COMMENTlink modified 3.1 years ago by Ram32k • written 3.1 years ago by babak114620

Please use the formatting bar (especially the code option) to present your post better. I've done it for you this time. Formatting bar

ADD REPLYlink written 3.1 years ago by Ram32k

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.

ADD REPLYlink written 3.1 years ago by Kevin Blighe70k
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: 978 users visited in the last hour