###################################################
### chunk number 1: 
###################################################
#line 48 "vignettes/affydata/inst/doc/affydata.Rnw"
library(affydata)
data(Dilution)


###################################################
### chunk number 2: 
###################################################
#line 62 "vignettes/affydata/inst/doc/affydata.Rnw"
Dilution


###################################################
### chunk number 3: 
###################################################
#line 70 "vignettes/affydata/inst/doc/affydata.Rnw"
phenoData(Dilution)
pData(Dilution)


###################################################
### chunk number 4: 
###################################################
#line 84 "vignettes/affydata/inst/doc/affydata.Rnw"
pData(Dilution) ##notice the scanner covariate


###################################################
### chunk number 5: 
###################################################
#line 95 "vignettes/affydata/inst/doc/affydata.Rnw"
par(mfrow=c(1,1))
boxplot(Dilution,col=c(2,2,3,3))


###################################################
### chunk number 6: 
###################################################
#line 118 "vignettes/affydata/inst/doc/affydata.Rnw"
options(warn=-1)


###################################################
### chunk number 7: 
###################################################
#line 123 "vignettes/affydata/inst/doc/affydata.Rnw"
gn <- sample(geneNames(Dilution),100) ##pick only a few genes
pms <- pm(Dilution[,3:4], gn)
mva.pairs(pms)


###################################################
### chunk number 8: 
###################################################
#line 135 "vignettes/affydata/inst/doc/affydata.Rnw"
normalized.Dilution <- Biobase::combine(normalize(Dilution[, 1:2]),
                             normalize(Dilution[, 3:4]))


###################################################
### chunk number 9: 
###################################################
#line 146 "vignettes/affydata/inst/doc/affydata.Rnw"
normalize.methods(Dilution)


###################################################
### chunk number 10: 
###################################################
#line 159 "vignettes/affydata/inst/doc/affydata.Rnw"
boxplot(normalized.Dilution, col=c(2,2,3,3), main="Normalized Arrays")


###################################################
### chunk number 11: 
###################################################
#line 170 "vignettes/affydata/inst/doc/affydata.Rnw"
pms <- pm(normalized.Dilution[, 3:4],gn)
mva.pairs(pms)