方法一：使用aggregate()分组获取描述性统计量
1 aggregate(mtcars[vars],by=list(am=mtcars$am),mean) 2 aggregate(mtcars[vars],by=list(mtcars$am),mean) 3 aggregate(mtcars[vars],by=list(am=mtcars$am),sd)
若有多个分组变量，使用by=list(name1=groupvar1, name2=groupvar2, ... , groupvarN)
这样只能返回单个统计量，若要返回多个统计量，可使用by（）函数
原文：http://www.cnblogs.com/wujiadong2014/p/4925720.html