[BioC] followup on reading RNA digestion plots

Swati Ranade ssr3072 at yahoo.com
Tue Mar 29 02:11:07 CEST 2005


Thanks Naomi. Will try your suggestion, but I get the
point, the trends of the curves are more important
than the slopes alone.
--- Naomi Altman <naomi at stat.psu.edu> wrote:
> Looking at the plot, I would not take the slopes too
> seriously!
> 
> I only eliminate an array if it appears to be
> dramatically different from 
> the others.  Your arrays vary (which is expected)
> but I do not see any 
> outlying array.
> 
> I suppose you could do a formal test, but I wouldn't
> bother.  To do the 
> test, you would probably want to use a
> semiparametric model of the form
> 
> curve(i) = m(i) + s(x)
> 
> where s(x) is the mean curve and m(i) is a random
> array effect.  Once you 
> fit the model, you could test the goodness of fit of
> each curve to the 
> model, adjust for multiple comparisons, and reject
> any significant 
> curves.  Sorry, I do not have a reference for this.
> 
> --Naomi
> 
> At 08:35 AM 3/28/2005, Swati Ranade wrote:
> 
> >   Hello,
> >
> >  A while ago I had posted a question about
> amplified>
> >RNA experiment and  guidelines one should follow
> >during data analysis and Dr. Gentleman had
> suggested
> >that RNA digestion plot will be a good indicator
> for
> >deciding which samples should be used in data
> analysis
> >as it is important to comapre likes. For my
> >experiments,I think prima facie the plots look
> quite a
> >bit similar but since I am new to this area
> >I cannot estimate what should be range of
> difference
> >in slopes that one should tolerate while choosing
> the
> >comparable profiles. To give an idea of what I am
> >seeing I am attaching a postscript image and giving
> >the slopes at the bottom of this email. I will be
> >extremely obliged if somebody could provide me
> >further guidance so that I know I am on the right
> >track. Please let me know if I need to provide more
> >than the following information?
> >
> >  Thanks,
> >
> >  Swati
> >
> >Here is the data: (I think sample # 8 seems to have
> >a slope that is different fro the rest!! Should it
> be
> >eliminated or one should observe the trend of the
> >graph?)
> >$slope
> >  [1] 1.6329065 1.5580935 1.5640235 0.9613637
> >  1.2156617 1.2004132 1.4831510
> >  [8] 0.3772531 1.6212330 1.2952673 1.2148458
> >  1.6020575 1.8564214 1.0944649
> >
> >
> >
> >
> >__________________________________________________
> >
> >
> >
> >_______________________________________________
> >Bioconductor mailing list
> >Bioconductor at stat.math.ethz.ch
> >https://stat.ethz.ch/mailman/listinfo/bioconductor
> 
> Naomi S. Altman                               
> 814-865-3791 (voice)
> Associate Professor
> Bioinformatics Consulting Center
> Dept. of Statistics                             
> 814-863-7114 (fax)
> Penn State University                        
> 814-865-1348 (Statistics)
> University Park, PA 16802-2111
> 
>



More information about the Bioconductor mailing list