### R code from vignette source 'vignettes/RPA/inst/doc/RPA.Rnw'

###################################################
### code chunk number 1: rpa
###################################################
require(affy)
require(affydata)
data(Dilution)


###################################################
### code chunk number 2: rpa (eval = FALSE)
###################################################
## eset <- rpa(Dilution)


###################################################
### code chunk number 3: RPA.pointestimate
###################################################
require(RPA)
sets <- geneNames(Dilution)[1:2]
rpa.results <- RPA.pointestimate(Dilution, sets)


###################################################
### code chunk number 4: visu (eval = FALSE)
###################################################
## plot(rpa.results, set = "1000_at", plots = "all")
## #plot(rpa.results, set = set, plots = "data")
## #plot(rpa.results, set = set, plots = "toydata.comparison")


###################################################
### code chunk number 5: barplots
###################################################
dat <- plot(rpa.results, set = "1000_at", plots = "all")


###################################################
### code chunk number 6: noise
###################################################
noise <- get.probe.noise.estimates(rpa.results)


###################################################
### code chunk number 7: setpriors
###################################################
priors <- list(alpha = NULL, beta = NULL, d = NULL)
set <- "1000_at"
priors$alpha <- 2
priors$beta[[set]] <- rep(1, 16)
probe.index <- 5
priors$beta[[set]][[probe.index]] <- 3
rpa.results <- RPA.pointestimate(Dilution, sets, priors = priors)


###################################################
### code chunk number 8: iter (eval = FALSE)
###################################################
## res <- rpa.fit(S)


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