###################################################
### chunk number 1: 
###################################################
library("yaqcaffy")
library("affydata")
data(Dilution)
## probe intensities modification
tmp <- exprs(Dilution)
tmp[,1] <- tmp[,1]*2
exprs(Dilution) <- tmp


###################################################
### chunk number 2: 
###################################################
yqc <-yaqc(Dilution)
show(yqc)


###################################################
### chunk number 3: 
###################################################
plot(yqc)


###################################################
### chunk number 4: 
###################################################
getOutliers(yqc,"sfs")


###################################################
### chunk number 5: 
###################################################
yqc2<-yaqc(Dilution[,2:3])
arrays(yqc) 
arrays(yqc2)
yqc3<-merge(yqc,yqc2)
arrays(yqc3)


###################################################
### chunk number 6: 
###################################################
library(MAQCsubsetAFX)
data(refB)
d<-refB[,1]
sampleNames(d)


###################################################
### chunk number 7:  eval=FALSE
###################################################
##  reprodPlot(d,"refA",normalize="rma")


###################################################
### chunk number 8: 
###################################################
reprodPlot(d,"test",normalize="none")


###################################################
### chunk number 9: 
###################################################
sessionInfo()