[BioC] differences between "manual" boxplots and "standard" boxplots of a oligo object

Javier Pérez Florido jpflorido at gmail.com
Wed Apr 20 09:04:56 CEST 2011


Dear list,
I'm checking the quality of a set of raw Human Gene ST 1.0 arrays 
through the use of boxplots and histograms. CEL files are read using 
read.celfiles from oligo package.
Then, boxplots and histograms are plot using "boxplot" and "hist" 
commands respectively, but I observe some differences between the plots 
given by these commands and the plots of "manual" boxplots and 
histograms such these:

pms = log2(pm(OligoRaw))
dens = oligo:::matDensity(pms)
cols = darkColors(ncol(pms))
boxplot(as.data.frame(pms), range=0, col=cols)
matplot(dens[["x"]], dens[["y"]], type="l", lwd=2, 
col=cols,xlab="log-intensity", ylab="density")

where "OligoRaw" is the object that holds raw data from the CEL files.
To be more precise, in boxplots there are differences in the values 
related to Q3+1.5*IQR. In the "manual" boxplot, these values are bigger 
than in the standard boxplot command.
On the other hand, in histograms, the differences are related to the CEL 
files with the biggest density values. In the "manual" histogram 
(matplot), these values are bigger than in the standard hist command.

Why do these differences exist?
Thanks,
Javier



More information about the Bioconductor mailing list