###################################################
### chunk number 1:  eval=FALSE
###################################################
## library(ArrayExpress)
## getAE(c('E-GEOD-11121','E-GEOD-12093'), extract=FALSE)


###################################################
### chunk number 2:  eval=FALSE
###################################################
## library(arrayQualityMetrics)
## celDir <- 'PATH'
## celf <- list.celfiles(celDir, full.names=T)
## ab <- read.affybatch(celf)
## arrayQualityMetrics(ab, outdir='QA')


###################################################
### chunk number 3:  eval=FALSE
###################################################
## library(vsn)
## fit <- vsn2(ab, subsample=100000)
## x <- predict(fit, newdata=ab)


###################################################
### chunk number 4:  eval=FALSE
###################################################
## set.seed(1234)
## fit_ref <- vsn2(ab[,1:20])
## fit_add1 <- vsn2(ab[,21:40], reference=fit_ref)
## fit_add2 <- vsn2(ab[,41:60], reference=fit_ref)
## set.seed(1234)
## fit_comp <- vsn2(ab)