分组计算描述性统计量来看一个栗子 myvars-c("mpg","hp","wt")aggregate(mtcars[myvars],by=list(am=mtcars$am),mean) #均值 am mpg hp wt1 0 17. 160.2632 3. 1 24. 126.8462 2.aggregate(mtcars[myvars],by=list(am=mtcars$am),sd) #标准差 am mpg hp wt1 0 3. 53. 0. 1 6. 84.06232 0.#上述的描述了以am分组标签,对数据框的mpg,hp及wt进行了均值mea和标准差sd的描述,但是很遗憾,aggregate函数每次只能单一调用和返回相关的统计量,所有#我们需要借助于by函数来解决这个问题,格式为:#by(data,indices,f n),其中data为数据框或者矩阵,indices是一个因子或者因子组成的列表,f n是任意函数.再来看一个栗子 dstats-f nction(x) sapply(x,mystats) #定义函数,调用函数 myvars-c("mpg","hp","wt") by(mtcars[myvars],mtcars$am,dstats) #by函数结果如下图: by函数统计结果 引用phych包中的 describeby分组计算概述统计量library(hmisc) myvars - c("mpg", "hp", "wt") describeby(mtcars[myvars], list(am=mtcars$am)) #中间必须有空格结果如下: 分组描述结果 好了,我的小伙伴们,今天就先到这儿吧,下期见!