某一天，我生成了2个箱图boxplot，突然间老师要求合并成1个箱图，唉～～～于是我各种查资料，学着做这种分组plot，结果如上图所示，主要想用多组boxplot展示同组数据不同分类之间的区别，以及不同组之间的区别。
以上是闲聊
以下是代码
定义个性化boxplot函数
这个函数模块是生成多组boxplot的核心模块，功能是生成没有坐标轴的箱图，有它才能顺利生成想要的箱图。
#定义boxplot的函数
myboxplot <- function(x, col = NULL, xlab, pvalue="auto") {
boxplot(x, axes = FALSE, col = col)
axis(1, at = 1:2, labels =FALSE)
text(1:2, y=par()$usr[3]-0.08*(par()$usr[4]-par()$usr[3]),
srt=60, xpd=T, adj=1, labels = xlab)
#srt 字体旋转角度；usr 绘图区域的用户坐标的极值。，
#xpd 如果为FALSE，则所有绘图都剪切到绘图区域，如果为TRUE，则所有绘图都剪切到图形区域，
#如果为NA，则所有绘图都剪切到设备区域
#adj 值为0左对齐文本，0.5（默认）居中文本和1右对齐文本。
这里采用定义函数的原因是，在函数内可以对boxplot进行个性化的定义。例如函数myboxplot代码中就定义boxplot没有坐标轴(axes = FALSE)， 对注释进行了个性化定义，例如图形的画图区域，文本的旋转，右对齐等。
boxplot画面布局
我的数据共有30组，每组分为2类数据。因此需要划分boxplot画图区域为30块，采用layout将画画空间分为30块，并排成一行，这样加上个性化的没有坐标轴的boxplot，最后生成的图片如首页图所示。
#获得30个位点的数据,第一个位点特殊，专门列出来boxplot
jpeg('这是文件名.jpg',width = 1000,height = 600)#储存图片到电脑
layout(t(1:30))#生成30个空白区域
par(oma=c(5, 2, 0.1, 0),mar=c(5,0.6,0.2,0.2), cex=1)#设置空白区域的图形边界
#30个多组肉眼不太容易看到获得颜色
tcol <- rainbow(30)
总体思路是生成30块空白画布，1个画布代表着一个组的数据，30组数据总共需要画30次，需要采取for循环，每次画一块画布。第一个图较为特殊，画完之后需要另外添加y轴。
我的原始数据是2个单boxplot，下面是从2个单boxplot的值列表和标志列表中根据组，分别筛选出2个单boxplot中的该组的数据，并且合并为新值列表temp_fre和新标志列表temp_label，随后进行空白画布的画图。
p = 0
#获得数据，并画图
for(i in 1:length(h3site)){
#获得该位点site的fre和label
tempsitefre <- c()
tempsitelabel <- c()
for(m in 1:length(site_fre)){
if(as.character(site_lab[m]) == as.character(h3site[i])){
tempsitefre <- c(tempsitefre, site_fre[m])
tempsitelabel <- c(tempsitelabel, site_lab[m])
}
}
print(length(tempsitefre))
print(length(tempsitelabel))
#获得该位点allsite的fre和label
tempallsitefre <- c()#储存该位点allsite 2组数据的fre
tempallsitelabel <- c()#储存该位点allsite 2组数据的label
for(m in 1:length(asite_fre)){
if(as.character(asite_lab[m]) == paste('as', as.character(h3site[i]), sep = '')){
tempallsitefre <- c(tempallsitefre, asite_fre[m])
tempallsitelabel <- c(tempallsitelabel, asite_lab[m])
}
}
#为了做出1个位点2组数据，需要合并site和allsite该位点的数据
temp_fre <- c()#储存该位点2组数据的fre
temp_label <- c()#储存该位点2组数据的label
temp_fre <- c(tempsitefre, tempallsitefre)
temp_label <- c(tempsitelabel,tempallsitelabel)
print(length(temp_fre))
print(length(temp_label))
print(head(temp_fre))
print(head(temp_label))
print(tail(temp_fre))
print(tail(temp_label))
#获得该位点总数据以后，进行boxplot，第一个boxplot比较特殊，因此采用if条件语句
if(i == 1){
p = myboxplot(temp_fre ~ temp_label,col = tcol[i], xlab = c(h3site[i], paste('as', as.character(h3site[i]), sep = '')))
p
axis(2, las=1)#画出y轴
}else{
p = myboxplot(temp_fre ~ temp_label,col = tcol[i], xlab = c(h3site[i], paste('as', as.character(h3site[i]), sep = '')))
p
}
}
dev.off()#结束画图
以上是我这期的代码分享，如果你觉得我写的内容对你有帮助的话，求赞❤️，求分享，求收藏，求关注➕，谢谢小伙伴们的支持，这里是明天吃什么饭的今日代码分享，我们下期见。
参考链接
主体思路：https://guangchuangyu.github.io/cn/2017/10/boxplot-multiple-group/
text函数的具体参数：http://www.mamicode.com/info-detail-1776049.html
layout函数用法：https://blog.csdn.net/mr_muli/article/details/80000449
par函数的具体参数：https://www.jianshu.com/p/5d44a08f67a5
oma和mar区别：https://www.cnblogs.com/xudongliang/p/6741843.html