> ddply(ma, .(variable), summarise, mean = mean(value), sd = sd(value), > skewness = skewness(value), median = median(value), > mean.gt.med = mean.gt.med(value))
In principle, you should be able to do: ddply(ma, .(variable), colwise(each(mean, sd, skewness, median, mean.gt.med))) but currently colwise and each don't work together that well. Hadley -- Assistant Professor / Dobelman Family Junior Chair Department of Statistics / Rice University http://had.co.nz/ ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.