方法一：使用aggregate()分组获取描述性统计量
aggregate(mtcars[vars],by=list(am=mtcars$am),mean)
aggregate(mtcars[vars],by=list(mtcars$am),mean)
aggregate(mtcars[vars],by=list(am=mtcars$am),sd)
若有多个分组变量，使用by=list(name1=groupvar1, name2=groupvar2, ... , groupvarN)
这样只能返回单个统计量，若要返回多个统计量，可使用by（）函数