Question: How to use detect call in GCRMA (or RMA) analysis?
0
gravatar for Wang mingyi
12.9 years ago by
Wang mingyi10
Wang mingyi10 wrote:
Dear all, I have a question for Affymetrix data analysis by GCRMA (or RMA). From the GCRMA (RMA) results, we cannot get detect call information which can be returned from MAS5. I think this information are useful when we want to perform some downstream t-tests. For example, three replicates (after drug) comparing with other three replicates (before drug), if given one gene, most detect values are "A" in all these 6 replicates, we will think the t-test doesnot make sense. So how to deal with this situation? Can we incorporate with detect call information from MAS5? Thank you in advance! _________________________________________________________________ ?????????????????????????????? MSN Hotmail?? http://www.hotmail.com
gcrma • 585 views
ADD COMMENTlink modified 12.9 years ago by Jianping Jin890 • written 12.9 years ago by Wang mingyi10
Answer: How to use detect call in GCRMA (or RMA) analysis?
0
gravatar for Kasper Daniel Hansen
12.9 years ago by
United States
Kasper Daniel Hansen6.4k wrote:
On Dec 6, 2006, at 2:47 PM, Wang mingyi wrote: > Dear all, > > I have a question for Affymetrix data analysis by GCRMA (or RMA). > From the GCRMA (RMA) results, we cannot get detect call information > which can be returned from MAS5. I think this information are > useful when we want to perform some downstream t-tests. For > example, three replicates (after drug) comparing with other three > replicates (before drug), if given one gene, most detect values are > "A" in all these 6 replicates, we will think the t-test doesnot > make sense. So how to deal with this situation? Can we incorporate > with detect call information from MAS5? > Thank you in advance! GCRMA or RMA does not make presence/absence calls like MAS5 does. One way to deal with this is to do the relevant routine and then subsequently remove the probesets which are declared absent using MAS5. You could also remove them before doing RMA/GCRMA but that is a bit harder. Kasper > _________________________________________________________________ > ??????????????? MSN Hotmail? http:// > www.hotmail.com > > _______________________________________________ > 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 COMMENTlink written 12.9 years ago by Kasper Daniel Hansen6.4k
Answer: How to use detect call in GCRMA (or RMA) analysis?
0
gravatar for Jianping Jin
12.9 years ago by
Jianping Jin890
Jianping Jin890 wrote:
You may want to try "mas5calls", a function in the affy package and then filter out the probesets with most absent calls. For example: PAcalls <- mas5calls(affybatch.obj) HTH, JP- --On Wednesday, December 06, 2006 10:47 PM +0000 Wang mingyi <mingyiwang at="" hotmail.com=""> wrote: > Dear all, > > I have a question for Affymetrix data analysis by GCRMA (or RMA). From > the GCRMA (RMA) results, we cannot get detect call information which can > be returned from MAS5. I think this information are useful when we want > to perform some downstream t-tests. For example, three replicates (after > drug) comparing with other three replicates (before drug), if given one > gene, most detect values are "A" in all these 6 replicates, we will think > the t-test doesnot make sense. So how to deal with this situation? Can we > incorporate with detect call information from MAS5? > Thank you in advance! > > _________________________________________________________________ > ??????????????? MSN Hotmail? > http://www.hotmail.com > ################################## Jianping Jin Ph.D. Bioinformatics scientist Center for Bioinformatics Room 3133 Bioinformatics building CB# 7104 University of Chapel Hill Chapel Hill, NC 27599 Phone: (919)843-6105 FAX: (919)843-3103 E-Mail: jjin at email.unc.edu
ADD COMMENTlink written 12.9 years ago by Jianping Jin890
If you have a PM-only array, the mas5calls function will not do. I've seen a couple of ad hoc strategies for calling P/A only from PM probes, but never accompanied with comparison to MAS5 or some other quality metric. Has anyone dealt with this problem? -Amit On Dec 7, 2006, at 8:21 AM, Jianping Jin wrote: > You may want to try "mas5calls", a function in the affy package and > then > filter out the probesets with most absent calls. > For example: > > PAcalls <- mas5calls(affybatch.obj) > > HTH, > > JP- > > --On Wednesday, December 06, 2006 10:47 PM +0000 Wang mingyi > <mingyiwang at="" hotmail.com=""> wrote: > >> Dear all, >> >> I have a question for Affymetrix data analysis by GCRMA (or RMA). >> From >> the GCRMA (RMA) results, we cannot get detect call information >> which can >> be returned from MAS5. I think this information are useful when we >> want >> to perform some downstream t-tests. For example, three replicates >> (after >> drug) comparing with other three replicates (before drug), if >> given one >> gene, most detect values are "A" in all these 6 replicates, we >> will think >> the t-test doesnot make sense. So how to deal with this situation? >> Can we >> incorporate with detect call information from MAS5? >> Thank you in advance! >> >> _________________________________________________________________ >> ??????????????? MSN Hotmail? >> http://www.hotmail.com >> > > > > ################################## > Jianping Jin Ph.D. > Bioinformatics scientist > Center for Bioinformatics > Room 3133 Bioinformatics building > CB# 7104 > University of Chapel Hill > Chapel Hill, NC 27599 > Phone: (919)843-6105 > FAX: (919)843-3103 > E-Mail: jjin at email.unc.edu > > _______________________________________________ > 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 REPLYlink written 12.9 years ago by Amit Bahl20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 121 users visited in the last hour