Jin_Lee_2021/Matlab_function/Jin_Errorbar.m

function [h, Mean, STD, STE] = Jin_Errorbar(x,y,varargin)


% Example)
% Width=0.3;
% Color.color=53; Color.alpha=0.4;
% h = JErrorbar(1,a.data1,Width,Color)

if nargin == 2
    Width = 0.35;
    Color.color=4; Color.alpha=1;
end

if nargin == 3
    Color.color=4; Color.alpha=1;
    Width = varargin{1};
end

if nargin == 4
    Width = varargin{1};
    Color = varargin{2};
end


c=mapcolor(100,0);
Mean=nanmean(y);
STD=nanstd(y);
STE=nanstd(y)/sqrt(length(y));

x1=x-Width;
x2=x+Width;
y1=0;
y2=Mean;
% Box
% h=fill([x1 x1 x2 x2],[y2 y1 y1 y2],[0.6 0.6 0.6]);
h=fill([x1 x1 x2 x2],[y2 y1 y1 y2],c(Color.color,:));
h.FaceAlpha=Color.alpha;
h.Marker='none';
h.EdgeColor='none';
h.EdgeAlpha=0;
% Line
l1=line([x1 x1],[y2 y1],'LineWidth',0.6,'LineStyle','-'); l1.Color=[0 0 0];%l1.Color=c(Color.color,:);
% l1=line([x1 x2],[y1 y1],'LineWidth',0.6,'LineStyle','-'); l1.Color=c(Color.color,:);
l1=line([x1 x2],[y2 y2],'LineWidth',0.6,'LineStyle','-'); l1.Color=[0 0 0];%l1.Color=c(Color.color,:);
l1=line([x2 x2],[y2 y1],'LineWidth',0.6,'LineStyle','-'); l1.Color=[0 0 0];%l1.Color=c(Color.color,:);

% Std
if Mean >= 0
    l1=line([x x],[Mean (Mean+STE)],'LineWidth',0.6,'LineStyle','-'); l1.Color=[0 0 0]; %l1.Color=c(Color.color,:);
    l1=line([(x-Width/3) (x+Width/3)],[(Mean+STE) (Mean+STE)],'LineWidth',0.6,'LineStyle','-'); l1.Color=[0 0 0]; %l1.Color=c(Color.color,:);
elseif Mean < 0
    l1=line([x x],[Mean (Mean-STE)],'LineWidth',0.6,'LineStyle','-'); l1.Color=[0 0 0]; %l1.Color=c(Color.color,:);
    l1=line([(x-Width/3) (x+Width/3)],[(Mean-STE) (Mean-STE)],'LineWidth',0.6,'LineStyle','-'); l1.Color=[0 0 0]; %l1.Color=c(Color.color,:);
end


% for i=1:length(y)
%     p1=plot(GetJitter(x,0.1),y(i),'.'); p1.Color=[0.8 0.2 0.2];
% end
% g=gca;
% g.YColor=[0 0 0];
% g.LineWidth=0.6;