1
0

plot3lines_models.m 2.4 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778
  1. function plot3lines_models(time,a,colores,xl,yl,legend1,legend2,titu,simple,xlab,addingcluster,timesignclust,coloresclus,rdline,rdlinecolor,legend3,showleg)
  2. a_plot(20)=plot(time,rdline,'--','LineWidth',2,'Color',rdlinecolor);hold on;
  3. for pl=[1 2];
  4. dia_acc=a{pl}
  5. timetoplot=time(150:end);
  6. mean_da=nanmean(dia_acc(150:end,:),2);
  7. sem_da=nanstd(dia_acc(150:end,:),0,2)/sqrt(sum(~isnan(dia_acc(1,:))));
  8. line([timetoplot(1),timetoplot(end)], [0,0], 'Color', 'k','LineStyle','-','LineWidth',2);hold on;
  9. % line([timetoplot(1),timetoplot(end)], [chancelevel,chancelevel], 'Color', 'k','LineStyle','--');hold on;
  10. a_plot(pl)=plot(timetoplot,mean_da,'Color',colores(pl,:),'LineWidth',3);hold on;
  11. fill([timetoplot,flipdim(timetoplot,2),timetoplot(1)],[mean_da+sem_da;flipdim(mean_da-sem_da,1);mean_da(1)+sem_da(1)],colores(pl,:),'LineStyle','none','FaceAlpha',0.30);
  12. hold on
  13. end
  14. if addingcluster
  15. % Adding lines of cluster perm test
  16. strp=-0.005;
  17. dist=0.005;
  18. for lll=1:1
  19. clustertoplot=nan(size(time))
  20. t1=find(time==timesignclust(lll,1))
  21. t2=find((time==timesignclust(lll,2)))
  22. clustertoplot(t1:t2)=strp-dist*lll
  23. hold on;
  24. a_plot(10+lll)=plot(time,clustertoplot,'o',...
  25. 'MarkerSize',5,...
  26. 'MarkerEdgeColor',coloresclus(lll,:),...
  27. 'MarkerFaceColor',coloresclus(lll,:))
  28. legend('off')
  29. end
  30. end
  31. if simple==0
  32. a_plot(5)=line([0,0],[-10 10], 'Color', 'k','LineStyle',':','LineWidth',2);hold on;
  33. a_plot(6)=line([1000,1000],[-10 10], 'Color', 'k','LineStyle',':','LineWidth',2);hold on;
  34. % a_plot(7)=line([2000,2000],[-10 10], 'Color', 'k','LineStyle','--','LineWidth',2);hold on;
  35. a_plot(8)=fill([2000 2350 2000 2350],[20 20 0 0],'k','LineStyle','none','FaceAlpha',0.2);hold on;
  36. a_plot(9)=fill([0 500 0 500],[20 20 0 0],'k','LineStyle','none','FaceAlpha',0.1);hold on;
  37. a_plot(10)=fill([1000 1500 1000 1500],[20 20 0 0],'k','LineStyle','none','FaceAlpha',0.1);hold on;
  38. if showleg
  39. legend([a_plot(1 ) a_plot(2) a_plot(20)],...
  40. legend1,legend2,legend3,'Location','NorthEast')
  41. legend('boxoff')
  42. end
  43. elseif simple==1
  44. legend([a_plot(1 ) a_plot(2)],...
  45. legend1,legend2,'Location','South')
  46. legend('boxoff')
  47. end
  48. title([titu])
  49. xlim(xl)
  50. ylim(yl)
  51. xlabel(xlab)
  52. ylabel('rho')
  53. set(gca, 'box', 'off')
  54. end