How to do RMA without summary to probeset level?
1
0
Entering edit mode
Peng Yu ▴ 940
@peng-yu-3586
Last seen 9.6 years ago
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")
• 1.2k views
ADD COMMENT
0
Entering edit mode
@benilton-carvalho-1375
Last seen 4.0 years ago
Brazil/Campinas/UNICAMP
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
ADD COMMENT
0
Entering edit mode
library(oligo) data<-read.celfiles(list.celfiles(data_dir, full.names=TRUE)) pm(data) Do you mean the above commands? If that is the case, you misunderstood me. I don't need the original probe values. As far as I know, RMA does three things: background correction, quantile normalization, and summary from probes to probesets. I want the probe values after background correction and quantile normalization but before the summary. On Sat, Dec 26, 2009 at 12:07 PM, Benilton Carvalho <bcarvalh at="" jhsph.edu=""> 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 > >
ADD REPLY
0
Entering edit mode
library(oligo) data = read.celfiles(list.celfiles(data_dir, full=T)) idx = oligo:::getFidProbeset(data) theProbes = exprs(data[idx,]) bgData = rma.background.correct(theProbes) normData = normalize.quantiles(log2(bgData)) On Dec 26, 2009, at 5:56 PM, Peng Yu wrote: > library(oligo) > data<-read.celfiles(list.celfiles(data_dir, full.names=TRUE)) > pm(data) > > Do you mean the above commands? If that is the case, you misunderstood > me. I don't need the original probe values. > > As far as I know, RMA does three things: background correction, > quantile normalization, and summary from probes to probesets. I want > the probe values after background correction and quantile > normalization but before the summary. > > On Sat, Dec 26, 2009 at 12:07 PM, Benilton Carvalho <bcarvalh at="" jhsph.edu=""> 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
ADD REPLY
0
Entering edit mode
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
ADD REPLY

Login before adding your answer.

Traffic: 731 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6