Jin_MeanSTE_Line.m 1.5 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354
  1. function [Mean, STD, STE] = Jin_MeanSTE_Line(x,y,varargin)
  2. y = GetExcludeNan(y);
  3. % markersize=25;
  4. markersize=5;
  5. linewidth=1;
  6. % Example)
  7. % Width=0.3;
  8. % Color.color=53; Color.alpha=0.4;
  9. % h = JMeanSTE_Line(1,a.data1,Width,Color)
  10. if nargin == 2
  11. Width = 0.25;
  12. Color.color=4; Color.alpha=1;
  13. end
  14. if nargin == 3
  15. Color = varargin{1};
  16. Width = 0.25;
  17. end
  18. if nargin == 4
  19. Width = varargin{1};
  20. Color = varargin{2};
  21. end
  22. c=mapcolor(100,0);
  23. Mean=nanmean(y);
  24. STD=nanstd(y);
  25. STE=nanstd(y)/sqrt(length(y));
  26. x1=x-Width;
  27. x2=x+Width;
  28. y1=0;
  29. y2=Mean;
  30. % p1=plot(x,Mean,'.'); p1.MarkerSize=markersize; p1.Color=c(Color.color,:);
  31. % Std
  32. if Mean >= 0
  33. l1=line([x x],[(Mean-STE) (Mean+STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:); l1.LineWidth=linewidth;
  34. l1=line([(x-Width/3) (x+Width/3)],[(Mean+STE) (Mean+STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:); l1.LineWidth=linewidth;
  35. l1=line([(x-Width/3) (x+Width/3)],[(Mean-STE) (Mean-STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:); l1.LineWidth=linewidth;
  36. elseif Mean < 0
  37. l1=line([x x],[Mean+STE (Mean-STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:);
  38. l1=line([(x-Width/3) (x+Width/3)],[(Mean+STE) (Mean+STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:); l1.LineWidth=linewidth;
  39. l1=line([(x-Width/3) (x+Width/3)],[(Mean-STE) (Mean-STE)],'LineWidth',1,'LineStyle','-'); l1.Color=c(Color.color,:);
  40. end
  41. g=gca;
  42. g.YColor=[0 0 0];
  43. g.XColor=[0 0 0];
  44. g.ZColor=[0 0 0];
  45. g.LineWidth=0.6;