I have hundreds of arrays for which I would like to obtain NUSE and RLE plots. I would like to make one call to affyPLM, but then plot a subset (boxplot), perhaps a few random samples that cover the entire set. Otherwise, there will be too many samples to show effectively in one boxplot. What is a way to go about this? Example below taken from affyPLM vignette.
require(affydata) library("affyPLM") data(Dilution) Dilution = updateObject(Dilution) Pset <- fitPLM(Dilution) boxplot(Pset)