[BioC] How to do RMA without summary to probeset level?

[Benilton Carvalho]

now that i think a bit more about this... if you want the control probes too, what i suggest is to use the following (yeah, yeah, it's not exposed):

idx = oligo:::getFidProbeset(data)
theProbes = exprs(data[idx,])

b

On Dec 26, 2009, at 4:07 PM, Benilton Carvalho wrote:

> pm(data)
> 
> b
> 
> On Dec 26, 2009, at 2:21 PM, Peng Yu wrote:
> 
>> I use the following code to do RMA. I'm wondering how get the probe
>> level values before the summary to the probeset level values.
>> 
>> library(oligo)
>> data<-read.celfiles(list.celfiles(data_dir, full.names=TRUE))
>> eset<-rma(data)
>> write.exprs(eset, file='output_file_name', sep="\t")
>> 
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor