Hi all, I'm working in published array expression data sets (affymetrix u133b) . I've tried to do some loess normalization but I found some problems. Here the code:
library(affy) celfiles <- list.files() rawData <- ReadAffy(filenames =celfiles) exprs(rawData) boxplot(rawData) normData <-normalize.AffyBatch.loess(rawData) exprs(normData) boxplot(normData)
So the code compiles properly but when I plot
normData the boxplot doesn't look like normalized.