Hi there

I would like to generate a histogramof my data after normalision using oligo
rma.  See code below.

library(oligo)
celfiles <- list.celfiles("cel_files",full.names = T)
twfs <- read.celfiles(celfiles)
results <- rma(twfs)
hist(results)

The hist implementation in this instance takes the log2 of the expression
set elements before making the plot.  This is a problem since the data has
already been logged, so my x-axis range is now [0,4] instead of [0,16].

Is there a way to square each element of "results" before calling hist?

Regards

Andrew

	[[alternative HTML version deleted]]

