在之前的文章中，分享了Matlab带缺口的分组箱线图的绘制模板：
进一步，再来分享一下带缺口的分组填充箱线图的绘制模板。
先来看一下成品效果：
特别提示：本期内容『数据+代码』已上传资源群中，加群的朋友请自行下载。有需要的朋友可以关注同名公号【阿昆的科研日常】，后台回复关键词【绘图桶】查看加入方式。
模板中最关键的部分内容：
1. 数据准备
此部分主要是读取原始数据。
% 读取数据load data.mat
2. 颜色定义
作图不配色就好比做菜不放盐，总让人感觉少些味道。
但颜色搭配比较考验个人审美，需要多加尝试。
这里直接使用TheColor配色工具中的XKCD配色库：
%% 颜色定义C = TheColor('xkcd',[260,426]);CA = [repmat(C(1,:),size(group1,2),1);repmat(C(2,:),size(group2,2),1)];
3. 带缺口的分组填充箱线图绘制
通过调用‘boxplot’和‘patch’命令，绘制初始带缺口的分组填充箱线图。
b1 = boxplot(group1,'positions',p1,'Colors','k','Notch','on','Widths',0.7,'Symbol','o');b2 = boxplot(group2,'positions',p2,'Colors','k','Notch','on','Widths',0.7,'Symbol','o');boxobj = findobj(gca, 'Tag', 'Box');for i = 1:length(boxobj)X = get(boxobj(i),'XData');Y = get(boxobj(i),'YData');p(i) = patch(X,Y,CA(i,:),'EdgeColor','k','FaceAlpha',0.7,'LineWidth',1.5);enddelete(b1,b2)box_1 = boxplot(group1,'positions',p1,'Colors','k','Notch','on',...'Widths',0.7,'Symbol','o','OutlierSize',5,'Jitter',0.2);box_2 = boxplot(group2,'positions',p2,'Colors','k','Notch','on',...'Widths',0.7,'Symbol','o','OutlierSize',5,'Jitter',0.2);hTitle = title('Miles per Gallon by Vehicle Origin');hXLabel = xlabel('Samples');hYLabel = ylabel('Miles per Gallon (MPG)');
4. 细节优化
为了插图的美观，对带缺口的分组填充箱线图线宽线型属性进行调整：
% 线宽set([box_1,box_2],'LineWidth',1.5)set([box_1(1:2,:),box_2(1:2,:)],'LineStyle','-')
然后，对坐标轴细节等进行美化：
% 坐标轴美化set(gca, 'Box', 'off', ... % 边框'LineWidth', 1,... % 线宽'XGrid', 'off', 'YGrid', 'on', ... % 网格'TickDir', 'out', 'TickLength', [.015 .015], ... % 刻度'XMinorTick', 'off', 'YMinorTick', 'off', ... % 小刻度'XColor', [.1 .1 .1], 'YColor', [.1 .1 .1]) % 坐标轴颜色set(gca, 'XTick', (p1+p2)/2,...'XTickLabel', {'Sample1','Sample2','Sample3','Sample4'},...'Xlim',[-0.5 11.5],...'Ylim',[-4 4]);% legendhLegend = legend([p(5),p(1)],...'Group1','Group2');% 字体和字号set(gca, 'FontName', 'Arial', 'FontSize', 9)set([hLegend,hXLabel,hYLabel], 'FontSize', 11, 'FontName', 'Arial')set(hTitle, 'FontSize', 12, 'FontWeight' , 'bold')% 背景颜色set(gcf,'Color',[1 1 1])% 添加上、右框线xc = get(gca,'XColor');yc = get(gca,'YColor');unit = get(gca,'units');ax = axes( 'Units', unit,...'Position',get(gca,'Position'),...'XAxisLocation','top',...'YAxisLocation','right',...'Color','none',...'XColor',xc,...'YColor',yc);set(ax, 'linewidth',1,...'XTick', [],...'YTick', []);
设置完毕后，以期刊所需分辨率、格式输出图片。
%% 图片输出figW = figureWidth;figH = figureHeight;set(figureHandle,'PaperUnits',figureUnits);set(figureHandle,'PaperPosition',[0 0 figW figH]);fileout = 'test';print(figureHandle,[fileout,'.png'],'-r300','-dpng');
以上。
