From: Gabor Grothendieck <ggrothendieck_at_gmail.com>

Date: Mon 01 May 2006 - 13:15:29 EST

*> with(DF, plot(ave(rate, bin), age))
*

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Mon May 01 14:29:55 2006

Date: Mon 01 May 2006 - 13:15:29 EST

My understanding is that you want to replace each rate with its average over the associated bin and then plot age against that. In that case try this:

1 0.002 10.0 A 2 0.045 0.1 B 3 0.130 15.0 A 4 0.150 34.0 D

Assuming they

are stored in vectors

the columns are age, rate, bin we would have

plot(ave(clock, bin), age)

On 4/30/06, lalitha viswanath <lalithaviswanath@yahoo.com> wrote:

*> Hi
*

> I am trying to get a table of means of parameter 1

*> across BINS of parameter 2.
**>
**> I am working in proteomics and a sample of my data is
**> as follows
**>
**> cluster-age clock-rate(evolutionary rate) scopclass
**> 0.002 10 A
**> 0.045 0.1 B
**> 0.13 15 A
**> 0.15 34 D
**> ....
**> ....
**> ....
**> ....
**>
**> Scop class has only 9 distinct categories (A-I)
**> Whereas cluster-age and clock-rate are discrete
**> variables greater than 0.
**>
**> I am trying to do two things with this kind of data,
**> out of which I managed to accomplish one thanks to the
**> documentation and pre-existing queries on the mailing
**> lists.
**> 1. Plot a histogram of the age distribution with scop
**> class category superimposed on each bin. I managed to
**> do this with barplot2.
**> 2. Now I am trying to plot a scatter plot of the age
**> v/s the clock-rate. However to eliminate possible
**> sampling errors, we are trying to get an average of
**> the clock-rate for each of the bins used above.
**> i.e. before plotting a x-y plot, i wish to compute
**> average clock-rate in each of the bins for the age and
**> then plot a x-y plot of the age v/s clock rate.
**>
**> Can anyone point me to appropriate functions for the
**> same?
**> I am trying to work with prop.table, cut, break, etc.
**> But I am not heading anywhere.
**>
**> Thanks
**> Lalitha
**>
**> ______________________________________________
**> R-help@stat.math.ethz.ch mailing list
**> https://stat.ethz.ch/mailman/listinfo/r-help
**> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
**>
*

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Mon May 01 14:29:55 2006

Archive maintained by Robert King, hosted by
the discipline of
statistics at the
University of Newcastle,
Australia.

Archive generated by hypermail 2.1.8, at Mon 01 May 2006 - 16:09:54 EST.

*
Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help.
Please read the posting
guide before posting to the list.
*