123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354 |
- function [Mean, STD, STE] = Jin_MeanSTE_Line(x,y,varargin)
- y = GetExcludeNan(y);
- % markersize=25;
- markersize=5;
- linewidth=1;
- % Example)
- % Width=0.3;
- % Color.color=53; Color.alpha=0.4;
- % h = JMeanSTE_Line(1,a.data1,Width,Color)
- if nargin == 2
- Width = 0.25;
- Color.color=4; Color.alpha=1;
- end
- if nargin == 3
- Color = varargin{1};
- Width = 0.25;
- 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;
- % p1=plot(x,Mean,'.'); p1.MarkerSize=markersize; p1.Color=c(Color.color,:);
- % Std
- if Mean >= 0
- l1=line([x x],[(Mean-STE) (Mean+STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:); l1.LineWidth=linewidth;
- l1=line([(x-Width/3) (x+Width/3)],[(Mean+STE) (Mean+STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:); l1.LineWidth=linewidth;
- l1=line([(x-Width/3) (x+Width/3)],[(Mean-STE) (Mean-STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:); l1.LineWidth=linewidth;
- elseif Mean < 0
- l1=line([x x],[Mean+STE (Mean-STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:);
- l1=line([(x-Width/3) (x+Width/3)],[(Mean+STE) (Mean+STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:); l1.LineWidth=linewidth;
- l1=line([(x-Width/3) (x+Width/3)],[(Mean-STE) (Mean-STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:);
- end
- g=gca;
- g.YColor=[0 0 0];
- g.XColor=[0 0 0];
- g.ZColor=[0 0 0];
- g.LineWidth=0.6;
|