[BioC] Affy: Present calls in an eset

Rafael A. Irizarry ririzarr at jhsph.edu
Wed Oct 8 11:53:41 MEST 2003


we will *try* to have mascall fro the next release.

On Wed, 8 Oct 2003 Arne.Muller at aventis.com wrote:

> Hi,
> 
> I get your point with interpreting absent/present calls. Technically it's a
> nice feature, becasue one can just discard the majority of the genes on the
> chip for further analysis. In fact I think absent/present calls make sense in
> terms of biology, since just a fraction of the genes are realy expressed at a
> time. How to express this numerially is a different story (and I guess a
> difficult one).
> 
> Anyway, with MAS the calls are calculated anyway, can't they? So, I'd be nice
> (at least for "completness") to add a "mascall" method to the exprSet objects
> generated by affy. What do you think?
> 
> By the way, if you ignore the call, do you set an arbitrary intensity cutoff
> later in your analysis, or do just reley on the statistics (anova  p-value or
> whatever)?
> 
> 	regards,
> 
> 	Arne
> 
> 
> > -----Original Message-----
> > From: A.J. Rossini [mailto:rossini at blindglobe.net]
> > Sent: 08 October 2003 15:33
> > To: Muller, Arne PH/FR
> > Cc: bioconductor at stat.math.ethz.ch
> > Subject: Re: [BioC] Affy: Present calls in an eset
> > 
> > 
> > <Arne.Muller at aventis.com> writes:
> > 
> > > Hello,
> > >
> > > I'm quite new to Bioconductor/affy, and I was wondering if 
> > there's a simple
> > > way to include the absent/present call for a gene in the 
> > outputfile generated
> > > with write.exprs(eset, file='boo') in theaffy package.
> > >
> > > the eset was generated with
> > >
> > > eset <- expresso(cel, bgcorrect.method = 'rma', 
> > normalize.method =   
> > >              'qspline', pmcorrect.method = 'pmonly', 
> > >              summary.method='liwong')
> > >
> > > For further analyses I'd like to exclude genes that are 
> > absent in all chips.
> > 
> > That's tough.  It isn't clear what a sensible definition of absent
> > is.  Or present.  
> > 
> > Do you mean "expressed" ?  "Differentially expressed" ?  "sort of
> > differentially expressed but not too weakly expressed?".  For any of
> > these, you'll need a precise definition (there isn't any in
> > Bioconductor), and you can compute your own.
> > 
> > (I know that MAS will make these calls; I'm only familiar with Rosetta
> > Resolver's variant, and they don't really make sense to me -- to be
> > precise, I know numerically how they are derived, but fail to why they
> > realistically connect biologically or technologically without a great
> > deal of assumptions and a wild imagination).
> > 
> > best,
> > -tony
> > 
> > -- 
> > rossini at u.washington.edu            
> > http://www.analytics.washington.edu/ 
> > Biomedical and Health Informatics   University of Washington
> > Biostatistics, SCHARP/HVTN          Fred Hutchinson Cancer 
> > Research Center
> > UW (Tu/Th/F): 206-616-7630 FAX=206-543-3461 | Voicemail is unreliable
> > FHCRC  (M/W): 206-667-7025 FAX=206-667-4812 | use Email
> > 
> > CONFIDENTIALITY NOTICE: This e-mail message and any attachments may be
> > confidential and privileged. If you received this message in error,
> > please destroy it and notify the sender. Thank you.
> >
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://www.stat.math.ethz.ch/mailman/listinfo/bioconductor
>



More information about the Bioconductor mailing list