[BioC] Loading multiple Affy chips simultaneously

Paul Boutros Paul.Boutros at utoronto.ca
Sun Jan 28 00:57:42 CET 2007


I just realized I forgot to give some basic information on my system:
R 2.4.1
affy 1.12.2
gcrma 2.6.0
Windows XP

Paul

-----Original Message-----
From: Paul Boutros [mailto:Paul.Boutros at utoronto.ca] 
Sent: Saturday, January 27, 2007 1:59 PM
To: 'bioconductor at stat.math.ethz.ch'
Subject: Loading multiple Affy chips simultaneously

Hello,

I'm trying to analyze & merge some rat data that was generated on two
separate Affymetrix chips (RAE230-A and RAE230-2).  One chip (RAE230-2) is a
superset of the other.  Approximately equal numbers of animals were analyzed
with the smaller array (n=9) as the larger (n=10).  

I would like to:
a) pre-process (via GC-RMA) the two datasets separately
b) subset the summarized expression levels of the larger chip
c) merge the subsetted dataset with the smaller one
d) perform statistical analysis on the overall data (e.g. with limma)

My code loads the data like this:
data.batch1 <- ReadAffy(filenames=cel.files.batch1,
phenoData="phenodata_batch1.txt");
data.batch2 <- ReadAffy(filenames=cel.files.batch2,
phenoData="phenodata_batch2.txt");

Then does some diagnostics such as:
par( mfrow = c(1,2) );
hist(data.batch1, main = "Batch 1");
hist(data.batch2, main = "Batch 2");

Which generates the warnings:
Attaching package: 'rat2302cdf'

        The following object(s) are masked from package:rae230acdf :
         i2xy 

        The following object(s) are masked from package:rae230acdf :
         xy2i

I was not sure if this is a problem or not.  I later go on to use the
functions AffyRNAdeg and gcrma, and I am wondering if this masking will
interfere with their correct operation?  In other words, do I need to
process the two batches sequentially, or can I do it in parallel as I've
shown here?

Many thanks for any help,
Paul



More information about the Bioconductor mailing list