一、行列求和、平均值和频度

rowSums, colSums, rowMeans, colMeans 可以简单理解为按行或列求和或求均值，table把数字或字符当成因子统计频度，都相当简单：

1. > a <- array(rep(1:3, each=3), dim=c(3,3))
2. > a
3. [,1] [,2] [,3]
4. [1,]    1    2    3
5. [2,]    1    2    3
6. [3,]    1    2    3
7. > rowSums(a)
8. [1] 6 6 6
9. > colSums(a)
10. [1] 3 6 9
11. > table(a)
12. a
13. 1 2 3
14. 3 3 3

对于多维数组，rowSums, colSums, rowMeans, colMeans的使用稍为复杂点。它们的参数为：

1. colSums (x, na.rm = FALSE, dims = 1)
2. rowSums (x, na.rm = FALSE, dims = 1)
3. colMeans(x, na.rm = FALSE, dims = 1)
4. rowMeans(x, na.rm = FALSE, dims = 1)

其中dims为整数，表示哪个或哪些维数被看做行或列，对于row统计函数，dims+1及以后的维度被看做行，对于col函数，dims及以前的维度（1:dims）被看做列：

1. > b <- array(rep(1:3, each=9), dim=c(3,3,3))
2. > b
3. , , 1
4. [,1] [,2] [,3]
5. [1,]    1    1    1
6. [2,]    1    1    1
7. [3,]    1    1    1
8. , , 2
9. [,1] [,2] [,3]
10. [1,]    2    2    2
11. [2,]    2    2    2
12. [3,]    2    2    2
13. , , 3
14. [,1] [,2] [,3]
15. [1,]    3    3    3
16. [2,]    3    3    3
17. [3,]    3    3    3
18. > rowSums(b)
19. [1] 18 18 18
20. > rowSums(b,dims=1)
21. [1] 18 18 18
22. > rowSums(b,dims=2)
23. [,1] [,2] [,3]
24. [1,]    6    6    6
25. [2,]    6    6    6
26. [3,]    6    6    6
27. > colSums(b)
28. [,1] [,2] [,3]
29. [1,]    3    6    9
30. [2,]    3    6    9
31. [3,]    3    6    9
32. > colSums(b,dims=2)
33. [1]  9 18 27

table可以统计数字出现的频率，也可以统计其他可以被看做因子的数据类型：

1. > table(b)
2. b
3. 1 2 3
4. 9 9 9
5. > c <- sample(letters[1:5], 10, replace=TRUE)
6. > c
7. [1] "a" "c" "b" "d" "a" "e" "d" "e" "c" "a"
8. > table(c)
9. c
10. a b c d e
11. 3 1 2 2 2

如果参数不只一个，它们的长度应该一样，结果是不同因子组合的频度表：

1. > a <- rep(letters[1:3], each=4)
2. > b <- sample(LETTERS[1:3],12,replace=T)
3. > table(a,b)
4. b
5. a   A B C
6. a 0 3 1
7. b 3 0 1
8. c 1 1 2

二、apply系列函数：

如果我们关心的不仅仅是求和、平均值和频度这些指标的计算，可以用apply系列函数来处理，这些函数包括apply、lapply、sapply、vapply、tapply和mapply。这些函数的使用可以从目标数据类型和返回值类型两个方面进行了解。

1、apply函数：

这个函数的使用格式为：apply(X, MARGIN, FUN, ...)。它应用的数据类型是数组或矩阵，返回值类型由FUN函数结果的长度确定。

X参数为数组或矩阵；MARGIN为要应用计算函数的边/维，MARGIN=1为第一维（行），2为第二维（列），...；FUN为要应用的计算函数，后面可以加FUN的有名参数。比如，要按行或列计算数组a的标准差就可以这样：

1. > apply(a, MARGIN=1, FUN=sd)
2. [1] 1 1 1
3. > apply(a, MARGIN=2, FUN=sd)
4. [1] 0 0 0

MARGIN的长度可以不是1（多维应用），如果长度等于X的维数，应用到FUN函数的数据就只有一个值，结果没什么意义，甚至函数会获得无效值：

1. > apply(b, MARGIN=3, FUN=sum)
2. [1]  9 18 27
3. > apply(b, MARGIN=1:2, FUN=sum)
4. [,1] [,2] [,3]
5. [1,]    6    6    6
6. [2,]    6    6    6
7. [3,]    6    6    6
8. > apply(a, MARGIN=1:2, FUN=sd)
9. [,1] [,2] [,3]
10. [1,]   NA   NA   NA
11. [2,]   NA   NA   NA
12. [3,]   NA   NA   NA

上面我们使用的sd、sum或mean函数的返回值的向量长度都是1（每一次单独计算），apply函数结果的维数与MARGIN的向量长度相同；如果FUN函数返回值的长度不是1而是每次都为n，apply函数的结果是维度为c(n, dim(X)[MARGIN])：

1. > a
2. [,1] [,2] [,3]
3. [1,]    1    2    3
4. [2,]    1    2    3
5. [3,]    1    2    3
6. > apply(a, MARGIN=1, FUN=quantile, probs=seq(0,1, 0.25))
7. [,1] [,2] [,3]
8. 0%    1.0  1.0  1.0
9. 25%   1.5  1.5  1.5
10. 50%   2.0  2.0  2.0
11. 75%   2.5  2.5  2.5
12. 100%  3.0  3.0  3.0
13. > apply(a, MARGIN=2, FUN=quantile, probs=seq(0,1, 0.25))
14. [,1] [,2] [,3]
15. 0%      1    2    3
16. 25%     1    2    3
17. 50%     1    2    3
18. 75%     1    2    3
19. 100%    1    2    3

如果FUN函数返回值的长度不一样，情况就复杂了，apply函数的结果会是列表。

2、lapply、sapply和vapply函数：

这几个函数是一套，前两个参数都为X和FUN，其他参数在R的函数帮助文档里有相信介绍。它们应用的数据类型都是列表，对每一个列表元素应用FUN函数，但返回值类型不大一样。lappy是最基本的原型函数，sapply和vapply都是lapply的改进版。

2.1  lapply返回的结果为列表，长度与X相同

1. > scores <- list(YuWen=c(80,88,94,70), ShuXue=c(99,87,100,68,77))
2. > lapply(scores, mean)
3. $YuWen
4. [1] 83
5. $ShuXue
6. [1] 86.2
7. > lapply(scores, quantile, probs=c(0.5,0.7,0.9))
8. $YuWen
9. 50%  70%  90%
10. 84.0 88.6 92.2
11. $ShuXue
12. 50%  70%  90%
13. 87.0 96.6 99.6

2.2  sapply返回的结果比较“友好”，如果结果很整齐，就会得到向量或矩阵或数组

sapply是simplify了的lapply，所谓的simplify，是指对结果的数据结构进行了simplify，方便后续处理。

1. > sapply(scores, mean)
2. YuWen ShuXue
3. 83.0   86.2
4. > sapply(scores, quantile, probs=c(0.5,0.7,0.9))
5. YuWen ShuXue
6. 50%  84.0   87.0
7. 70%  88.6   96.6
8. 90%  92.2   99.6

2.3  vapply函数：对返回结果（value）进行类型检查的sapply

虽然sapply的返回值比lapply好多了，但可预测性还是不好，如果是大规模的数据处理，后续的类型判断工作会很麻烦而且很费时。vapply增加的FUN.VALUE参数可以直接对返回值类型进行检查，这样的好处是不仅运算速度快，而且程序运算更安全（因为结果可控）。下面代码的rt.value变量设置返回值长度和类型，如果FUN函数获得的结果和rt.value设置的不一致（长度和类型）都会出错：

1. > probs <- c(1:3/4)
2. > rt.value <- c(0,0,0)  #设置返回值为3个数字
3. > vapply(scores, quantile, FUN.VALUE=rt.value, probsprobs=probs)
4. YuWen ShuXue
5. 25%  77.5     77
6. 50%  84.0     87
7. 75%  89.5     99
8. > probs <- c(1:4/4)
9. > vapply(scores, quantile, FUN.VALUE=rt.value, probsprobs=probs)

错误于vapply(scores, quantile, FUN.VALUE = rt.value, probs = probs) :

值的长度必需为3，

但FUN(X[[1]])结果的长度却是4

1. > rt.value <- c(0,0,0,0)   #返回值类型为4个数字
2. > vapply(scores, quantile, FUN.VALUE=rt.value, probsprobs=probs)
3. YuWen ShuXue
4. 25%   77.5     77
5. 50%   84.0     87
6. 75%   89.5     99
7. 100%  94.0    100
8. > rt.value <- c(0,0,0,'')  #设置返回值为3个数字和1个字符串
9. > vapply(scores, quantile, FUN.VALUE=rt.value, probsprobs=probs)

错误于vapply(scores, quantile, FUN.VALUE = rt.value, probs = probs) :

值的种类必需是'character'，

但FUN(X[[1]])结果的种类却是'double'

FUN.VALUE为必需参数。

3、 mapply函数：

R的在线文档说mapply是sapply的多变量版本（multivariate sapply），但它的参数顺序和sapply却不一样：

mapply(FUN, ..., MoreArgs = NULL, SIMPLIFY = TRUE, USE.NAMES = TRUE)

mapply应用的数据类型为向量或列表，FUN函数对每个数据元素应用FUN函数；如果参数长度为1，得到的结果和sapply是一样的；但如果参数长度不是1，FUN函数将按向量顺序和循环规则（短向量重复）逐个取参数应用到对应数据元素：

1. > sapply(X=1:4, FUN=rep, times=4)
2. [,1] [,2] [,3] [,4]
3. [1,]    1    2    3    4
4. [2,]    1    2    3    4
5. [3,]    1    2    3    4
6. [4,]    1    2    3    4
7. > mapply(rep, x = 1:4, times=4)
8. [,1] [,2] [,3] [,4]
9. [1,]    1    2    3    4
10. [2,]    1    2    3    4
11. [3,]    1    2    3    4
12. [4,]    1    2    3    4
13. > mapply(rep, x = 1:4, times=1:4)
14. [[1]]
15. [1] 1
16. [[2]]
17. [1] 2 2
18. [[3]]
19. [1] 3 3 3
20. [[4]]
21. [1] 4 4 4 4
22. > mapply(rep, x = 1:4, times=1:2)
23. [[1]]
24. [1] 1
25. [[2]]
26. [1] 2 2
27. [[3]]
28. [1] 3
29. [[4]]
30. [1] 4 4

4、tapply 和 by 函数：

tapply函数可以看做是table函数的扩展：table函数按因子组合计算频度，而tapply可以按因子组合应用各种函数。使用格式为：tapply(X, INDEX, FUN = NULL, ..., simplify = TRUE)

X为要应用函数的数据，通常为向量；INDEX为因子，和table函数一样，它的长度必需和X相同。

1. > (x <- 1:10)
2. [1]  1  2  3  4  5  6  7  8  9 10
3. > (f <- gl(2,5, labels=c("CK", "T")))
4. [1] CK CK CK CK CK T  T  T  T  T
5. Levels: CK T
6. > tapply(x, f, length)  #FUN函数是length，得到的结果和table类似
7. CK  T
8. 5  5
9. > table(f)
10. f
11. CK  T
12. 5  5
13. > tapply(x, f, sum)
14. CK  T
15. 15 40

by函数是tapply函数针对数据框类型数据的应用，但结果不怎么友好，你可以用下面语句看看情况：

1. with(mtcars, by(mtcars, cyl, summary))

三、aggregate函数

这个函数的功能比较强大，它首先将数据进行分组（按行），然后对每一组数据进行函数统计，最后把结果组合成一个比较nice的表格返回。根据数据对象不同它有三种用法，分别应用于数据框（data.frame）、公式（formula）和时间序列（ts）：

1. aggregate(x, by, FUN, ..., simplify = TRUE)
2. aggregate(formula, data, FUN, ..., subset, nana.action = na.omit)
3. aggregate(x, nfrequency = 1, FUN = sum, ndeltat = 1, ts.eps = getOption("ts.eps"), ...)

我们通过 mtcars 数据集的操作对这个函数进行简单了解。mtcars 是不同类型汽车道路测试的数据框类型数据：

1. > str(mtcars)
2. 'data.frame':   32 obs. of  11 variables:
3. $ mpg : num  21 21 22.8 21.4 18.7 18.1 14.3 24.4 22.8 19.2 ...
4. $ cyl : num  6 6 4 6 8 6 8 4 4 6 ...
5. $ disp: num  160 160 108 258 360 ...
6. $ hp  : num  110 110 93 110 175 105 245 62 95 123 ...
7. $ drat: num  3.9 3.9 3.85 3.08 3.15 2.76 3.21 3.69 3.92 3.92 ...
8. $ wt  : num  2.62 2.88 2.32 3.21 3.44 ...
9. $ qsec: num  16.5 17 18.6 19.4 17 ...
10. $ vs  : num  0 0 1 1 0 1 0 1 1 1 ...
11. $ am  : num  1 1 1 0 0 0 0 0 0 0 ...
12. $ gear: num  4 4 4 3 3 3 3 4 4 4 ...
13. $ carb: num  4 4 1 1 2 1 4 2 2 4 ...

先用attach函数把mtcars的列变量名称加入到变量搜索范围内，然后使用aggregate函数按cyl（汽缸数）进行分类计算平均值：

1. > attach(mtcars)
2. > aggregate(mtcars, by=list(cyl), FUN=mean)
3. Group.1      mpg cyl     disp        hp     drat       wt     qsec        vs        am     gear     carb
4. 1       4 26.66364   4 105.1364  82.63636 4.070909 2.285727 19.13727 0.9090909 0.7272727 4.090909 1.545455
5. 2       6 19.74286   6 183.3143 122.28571 3.585714 3.117143 17.97714 0.5714286 0.4285714 3.857143 3.428571
6. 3       8 15.10000   8 353.1000 209.21429 3.229286 3.999214 16.77214 0.0000000 0.1428571 3.285714 3.500000

by参数也可以包含多个类型的因子，得到的就是每个不同因子组合的统计结果：

1. > aggregate(mtcars, by=list(cyl, gear), FUN=mean)
2. Group.1 Group.2    mpg cyl     disp       hp     drat       wt    qsec  vs   am gear     carb
3. 1       4       3 21.500   4 120.1000  97.0000 3.700000 2.465000 20.0100 1.0 0.00    3 1.000000
4. 2       6       3 19.750   6 241.5000 107.5000 2.920000 3.337500 19.8300 1.0 0.00    3 1.000000
5. 3       8       3 15.050   8 357.6167 194.1667 3.120833 4.104083 17.1425 0.0 0.00    3 3.083333
6. 4       4       4 26.925   4 102.6250  76.0000 4.110000 2.378125 19.6125 1.0 0.75    4 1.500000
7. 5       6       4 19.750   6 163.8000 116.5000 3.910000 3.093750 17.6700 0.5 0.50    4 4.000000
8. 6       4       5 28.200   4 107.7000 102.0000 4.100000 1.826500 16.8000 0.5 1.00    5 2.000000
9. 7       6       5 19.700   6 145.0000 175.0000 3.620000 2.770000 15.5000 0.0 1.00    5 6.000000
10. 8       8       5 15.400   8 326.0000 299.5000 3.880000 3.370000 14.5500 0.0 1.00    5 6.000000

公式（formula）是一种特殊的R数据对象，在aggregate函数中使用公式参数可以对数据框的部分指标进行统计：

1. > aggregate(cbind(mpg,hp) ~ cyl+gear, FUN=mean)
2. cyl gear    mpg       hp
3. 1   4    3 21.500  97.0000
4. 2   6    3 19.750 107.5000
5. 3   8    3 15.050 194.1667
6. 4   4    4 26.925  76.0000
7. 5   6    4 19.750 116.5000
8. 6   4    5 28.200 102.0000
9. 7   6    5 19.700 175.0000
10. 8   8    5 15.400 299.5000

上面的公式 cbind(mpg,hp) ~ cyl+gear 表示使用 cyl 和 gear 的因子组合对 cbind(mpg,hp) 数据进行操作。

aggregate在时间序列数据上的应用请参考R的函数说明文档。

R语言进阶之3：数据汇总/透视/提炼相关推荐

1. 【R语言进阶】对数据组排序

   R语言中的排序命令 jiaxue_jiyin_jiance<-c(8900,4500,1500,8900,8900,3500,18000,18000) 用sort()对佳学基因检测价格进行排序 ...

2. 价值1143元的《R语言统计分析微生物组数据(Statistical Analysis of Microbiome Data with R)》系列图书

   文章目录 <R语言统计分析微生物组数据> 本书简介 作者简介 章节简介 猜你喜欢 写在后面 <R语言统计分析微生物组数据> Statistical Analysis of Mi ...

3. R语言data.table导入数据实战:data.table中编写函数并使用SD数据对象

   R语言data.table导入数据实战:data.table中编写函数并使用SD数据对象 目录 R语言data.table导入数据实战:data.table中编写函数并使用SD数据对象 #data.t ...

4. R语言data.table导入数据实战：把data.frame数据转化为data.table数据

   R语言data.table导入数据实战:把data.frame数据转化为data.table数据 目录 R语言data.table导入数据实战:把data.frame数据转化为data.table数据 ...

5. R语言data.table导入数据实战:data.table使用by函数进行数据分组（aggregate）

   R语言data.table导入数据实战:data.table使用by函数进行数据分组(aggregate) 目录 R语言data.table导入数据实战:data.table使用by函数进行数据分组( ...

6. R语言单变量分析实战：汇总统计（Summary Statistics）、频率表（Frequency Table）、图表（charts: boxplot、histogram、density）

   R语言单变量分析实战:汇总统计(Summary Statistics).频率表(Frequency Table).图表(charts: boxplot.histogram.density) 目录

7. R语言dplyr包通过数据列的索引重命名数据列实战（Rename Column by Index Position）

   R语言dplyr包通过数据列的索引重命名数据列实战(Rename Column by Index Position) 目录 R语言dplyr包通过数据列的索引重命名数据列实战(Rename Colum ...

8. R语言rev函数对数据对象（向量、dataframe通过行或者列）反序实战

   R语言rev函数对数据对象(向量.dataframe通过行或者列)反序实战 目录 R语言rev函数对数据对象(向量.dataframe通过行或者列)反序实战

9. R语言ggplot2可视化时间序列数据并添加稳定趋势线识别数据的稳定趋势

   R语言ggplot2可视化时间序列数据并添加稳定趋势线识别数据的稳定趋势 目录 R语言ggplot2可视化时间序列数据并添

最新文章

1. 视频录制，压缩实现源码
2. html不支持css3,ie不支持css3吗？
3. 解决android 异常导致应用程序停止运行的错误
4. Logstash Introduction
5. python popen sqlplus_Python基于Select模型实现Popen输出
6. P2824-[HEOI2016/TJOI2016]排序【线段树,二分】
7. Flex Gumbo中如何自定义HSlider数据Tip样式的例子
8. Python 密文转换为明文
9. Vant Weapp小程序蹲坑之事件名称问题及解决方案
10. Mac 查看jdk版本
11. 数据预处理——以GSE3494为例
12. web前端期末大作业 html+css+javascript防天天生鲜官网网页设计实例 企业网站制作
13. 解决vscode没有文件图标（亲测有效）
14. 前端js正则表达式2
15. vsftp虚拟账户登录失败331 Please specify the password.
16. Python爬虫实战(一) QQ音乐评论爬取及可视化分析
17. 如何写好科研论文2021期末考试答案
18. ipam：allocate ip 172.16.3.7 to node x.x.x.112 allocator failed, provided IP is already allocated
19. RGB与HSV颜色空间转换
20. c oracle udt,c#调用Oracle UDT自定义结构数组类型源代码

热门文章

1. 芝加哥 计算机 录取,花样Offer来袭~藤门收获芝大、哥大、杜克、宾大、康奈尔、帝国理工等多枚录取！...
2. flash 脚本 2
3. 共享软件开发者、独立软件开发者
4. (课程类 Course)
5. VMware Workstation download
6. UVA - 12627 Erratic Expansion(分治)
7. 熟悉超声心动图对心脏功能的评估及超声心动图的正常值
8. DE、SaDE、JADE、SHADE、L-SHADE算法整理
9. 解析大型.NET ERP系统：十三种界面设计模式
10. 什么是前端脚手架？脚手架原理？