[BioC] How to save jpegs of all MAplots of an affybatch

James W. MacDonald jmacdon at med.umich.edu
Wed Nov 23 18:13:57 CET 2005


James W. MacDonald wrote:

> If you really don't want these names, you could try your second attempt, 
> but add these two lines after the call to jpeg()
> 
> medchip <- rowMedians(log2(intensity(mysamples)[pms,]))
> MAplot(mysamples[,i], ref=medchip)
> 
Correction - this won't work - a closer reading of the code (or Ben's 
email) would have alerted me to the fact that the ref argument applies 
to the index of a reference chip rather than a set of reference data.

Best,

Jim



-- 
James W. MacDonald
Affymetrix and cDNA Microarray Core
University of Michigan Cancer Center
1500 E. Medical Center Drive
7410 CCGC
Ann Arbor MI 48109
734-647-5623



More information about the Bioconductor mailing list