在比较多组个体或观测时，关注的焦点经常是各组的描述性统计信息，而不是样本整体的描

述性统计信息。同样地，在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列将被标注为