【MATLAB】一组含几个bar的柱状图的绘制以及errorbar的添加

在数据统计分析时，往往会用到多个数据分组进行比较，这里介绍一种数据可视化方法。
& {  n. I0 j6 m1 R实例：
% Z4 I- U* G! q  f% 数据
volume_mean=[0.73,0.45;
+ Z5 Y( W3 ~% ^& c2 e                        0.42,0.43;
, d* }$ f( g4 j  S# J                        0.70,0.42];                        
( p; x" A; U; M0 Gvolume_std=[0.65,0.17;
                     0.35,0.14;
                     0.44,0.13];
%绘图                     
close all;figure;
& t& u* h. F* v* F: sh=bar(volume_mean);
set(h,'BarWidth',0.9);        % 柱状图的粗细
  `4 U  O  _" b8 h2 _hold on;
set(h(1),'facecolor',[139 35 35]./255) % 第一列数据视图颜色
set(h(2),'facecolor','k')        % 第二列数据视图颜色


ngroups = size(volume_mean,1);
* [8 R8 y' f0 x2 f0 q5 O. y# hnbars = size(volume_mean,2);
9 z. e( g- M9 l6 O+ k. o  ^groupwidth =min(0.8, nbars/(nbars+1.5));
. t( S$ Q  u2 Z# ]

) c5 Z' m+ F( A+ M% errorbar如果用不同颜色，可以利用colormap的颜色进行循环标记，这个例子没有用到colormap
8 m* u1 _2 N! Q6 P& j8 }% W%colormap(flipud([0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1]));   % blue / red
% color=[0 100/255 0; 220/255 20/255 60/255; 1 215/255 0; 0 0 1];
3 p* L, \& x) Z/ D0 Chold on;
for i = 1:nbars
    x = (1:ngroups) - groupwidth/2 + (2*i-1) * groupwidth / (2*nbars);
    errorbar(x,volume_mean(:,i),volume_std(:,i),'o','color',[.5 .5 .5],'linewidth',2);
end

  t( p# e' n% n, Q' ~
% j7 _; E, c! o  Bset(gca,'XTickLabel',{'2014','2015','2016'},'fontsize',14,'linewidth',2)
( \. b! ~9 y; P- G* v6 K/ wylim([0 1.5])
2 S2 _1 e+ a' o. O/ Fset(gca,'ytick',0:0.5:1.5)
& }1 x" `( ^7 @. m# `/ Y9 I/ c) g$ Vxylabel(gca,' ','Volume[Sv]')
legend('data1','data2','location','NorthEast')
# U( G- r. K7 o+ M3 f

( ?: }5 |3 \% i' ]# i! X以上实例可以参考使用。

+ b  K' |8 J6 _; g
errorbar的局部调整：
1.头部宽度调整
: r6 E+ Z. n* S% G4 c% Create errorbar
X = 0:pi/10:pi;
4 {" }( T2 a" X7 Q$ H; MY = sin(X) + 1;
! L+ W6 @$ m+ ]E = std(Y) * ones(size(X));
ha = errorbar(X, Y, E);
% Width of the top and bottom lines of errorbar
) a1 s- g: b9 n0 e3 Z! S% Fxlength = 0.2;
% Make horizontal lines with 'line'
; G' d0 K8 \% x+ P9 J+ k* qfor k = 1:length(X)
x = [X(k) - xlength, X(k) + xlength];
, d7 B, Q- @: N( C) C  L y_h = [Y(k) + E(k), Y(k) + E(k)];
line(x, y_h);
y_b = [Y(k) - E(k), Y(k) - E(k)];
line(x, y_b);
1 g6 {& ?0 m: R/ y/ f- m0 m; U& {end
# n' Q% o& L0 e参考：www.52ocean.cn
) x- f$ n: u( _; P; p$ w' F

$ i( g$ u0 S; c$ ^9 C6 t7 E! d                    

                                        转载本文请联系原作者获取授权，同时请注明本文来自叶瑞杰科学网博客。