[BioC] RNA degradation tends & options for analysis

Ben Bolstad bmb at bmbolstad.com
Thu Feb 23 19:25:32 CET 2006


An alternative way to the pairs() method mentioned below is 

MAplot(myAffyBatch,pairs=TRUE)

which will give pairwise MA plots (which tend to be a little more useful
for looking at things in an Y vs X type manner). However, once you get
beyond a relatively small number of arrays this plot becomes unwieldy.
In that case you might consider using it this way:

MAplot(myAffyBatch)

in this case each array will be MA plotted against a "reference" array
created using all the chips in the dataset.

Ben


On Thu, 2006-02-23 at 11:25 -0500, Naomi Altman wrote:
> One thing I find very useful is to look at
> 
> pairs(pm(myAffydata))
> 
> (This takes a while, since you are plotting lots of probes, and I 
> usually put in lower.panel=NULL to get only the upper triangle of plots.)
> 
> If the arrays are comparable, then most of the data should cluster 
> pretty tightly on the diagonal.
> 
> Incidentally, if some ambitious person would write a pairs routine 
> for hexbin, that would be both faster and more informative.
> 
> --Naomi
>



More information about the Bioconductor mailing list