[BioC] RMA on a batch with NA values
Tarca, Adi
atarca at med.wayne.edu
Fri Aug 8 21:53:06 CEST 2008
Hi all,
I am trying to run RMA using only a limited number of probes. I have the
x and y coordinates of those probes that I want to use.
The approach I have tried was to replace with NA all probe intensities
that I want to discard, but the RGUI gives an error and tries to end the
session. The console is still running but it does not return any
results.
Here is the code I use:
#####
abatch<-ReadAffy(filenames=fnss,celfile.path="./CELL")
indok<-xy2indices(x=okx,y=oky,abatch=abatch)
ints<-intensity(abatch)
ints[-indok,]<-NA
intensity(abatch)<-ints
esetb<-rma(abatch)
Any ideas?
Thanks,
Adi L. Tarca
More information about the Bioconductor
mailing list