[BioC] split arrays
scholz at Ag.arizona.edu
Thu Sep 29 08:21:08 CEST 2005
Thanks for this response. Unfortunately, there are no oligos in common between
the two arrays. If anyone else has a response to my question (below), I'd like
to hear it.
I am not familiar with the maize arrays, but I am using the following
procedure for Affymetrix moe430 split arrays, which have ~160 probesets
in common between A and B:
1) background-correct each chip separately at probe-level
2) get a measure of expression at probeset-level
3) plot the common probesets against each other for each pair of each
chips. If you observe the same thing as me, you will see that the trend
is linear but with intercept != 0 and slope != 1.
4) scale the B chip with those estimated intercept and slope
Steps 1 and 2 are easily done with rma( , normalize=F).
Wolfgang Huber and I are currently writing a little package which does
steps 3 and 4 automatically.
I'm not sure whether this procedure could make sense or be adapted
somehow to your maize arrays (do they have enough probes in common?),
but anyway, some food for thoughts...
> Recently you advised someone with a split set of maize arrays
> that they could do their analysis by reading all the A slides
> into an RGList and normalizing, then doing the same with the
> B slides, and then combining the two datasets via
> rbind() of the two MAList objects. I have a similar (the
> same?) set of arrays and some of the users of these arrays
> have noted that the A and B slides perform differently, i.e.
> more background on the B slide, for whatever reason. Though
> I'm not actually convinced this is true, it makes me wonder
> whether the two datasets should be combined at all since
> there may be a "between array set"
> source of variation. Am I right to segregate these sets or is
> there some overwhelming benefit to combining them? I'm no
> statistician and would appreciate your take.
College of Agriculture and Life Sciences Web Mail.
More information about the Bioconductor