Summary stats in table

Date: Sun, 23 Oct 2011 14:39:08 -0400
From: murdoch.duncan at gmail.com
To: R-help at r-project.org
Subject: [R] Summary stats in table

Suppose I have data like this:

A <- sample(letters[1:3], 1000, replace=TRUE)
B <- sample(LETTERS[1:2], 1000, replace=TRUE)
x <- rnorm(1000)

I can get a table of means via

tapply(x, list(A, B), mean)

and I can add the marginal means to this using cbind/rbind:

main <- tapply(x, list(A,B), mean)
Amargin <- tapply(x, list(A), mean)
Bmargin <- tapply(x, list(B), mean)

rbind(cbind(main, all=Amargin),all=c(Bmargin, mean(x)))

But this is tedious. Has some package got some code that makes this easier?

Duncan Murdoch