From: Duncan Murdoch <murdoch_at_stats.uwo.ca>

Date: Sun 09 Jul 2006 - 06:00:44 EST

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 Sun Jul 09 06:04:10 2006

Date: Sun 09 Jul 2006 - 06:00:44 EST

On 7/8/2006 3:44 PM, justin rapp wrote:

> I apologize for my constant questions but I am new to R and trying to

*> gain an appreciation for its capabilities. The following task is easy
**> in Excel and I was hoping somebody could give me a quick explanation
**> for how it can be acheived in R so I can avoid having to switch
**> between the two applications.
**>
**> How do I find the Summary Statistics in one Vector of the dataframe by
**> levels in another of the vectors.
**>
**> For example, I have the following headings for my data.frame.
**> Conference
**> Year Drafted
**> Height
**> Weight
**> Ratio
**>
**> I would like to see compute the mean Height, Weight, and Ratio as well
**> as their variances for each of the years under Year
**> Drafted(1980-2000). What is the most efficient way of doing this?
*

I think the quickest is

by(mydf, mydf$Year, summary)

but this won't give you the variance. You'll need your own little function to calculate mean and variance, e.g.

mysummary <- function(df) apply(df, 2,

function(x) c(mean=mean(x), variance=var(x)))

by(mydf, mydf$Year, mysummary)

If you don't like the format of the output, you can play around with the mysummary function. It will be applied to each subset of the data.frame, and the results will be put together into a list with one entry per level of mydf$Year.

Duncan

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 Sun Jul 09 06:04:10 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 Sun 09 Jul 2006 - 08:16:32 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.
*