在比较多组个体或观测时,关注的焦点经常是各组的描述性统计信息,而不是样本整体的描
述性统计信息。同样地,在R中完成这个任务有若干种方法。我们将以获取变速箱类型各水平的
描述性统计量开始。
vars<- c("mpg","hp","wt")
> aggregate(mtcars[vars],by=list(am=mtcars$am),median)
am mpg hp wt
1 0 17.3 175 3.52
2 1 22.8 109 2.32
注意list(am=mtcars$am)的使用。如果使用的是list(mtcars$am),则am列将被标注为