table of means/medians across bins used for a histogram
Or perhaps a bit simpler: plot(age ~ ave(clock, bin), DF)
On 4/30/06, Gabor Grothendieck <ggrothendieck at gmail.com> wrote:
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:
DF # test data
age rate bin 1 0.002 10.0 A 2 0.045 0.1 B 3 0.130 15.0 A 4 0.150 34.0 D
with(DF, plot(ave(rate, bin), age))
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 at 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 at 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