You can convert the resulting output to a matrix, and then link them:
a <- 1:50 b <- 3:53 c <- rnorm(500) cbind(as.matrix(summary(a)), as.matrix(summary(b)), as.matrix(summary(c)))
Alternatively, you can combine them into a list and use the apply (or plyr ) function:
library(plyr) ldply(list(a, b, c), summary)
Shane
source share