方法一：使用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（）函数