1
0

plot3lines_models2nd2.m 2.3 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273
  1. function plot3lines_models(time,a,colores,xl,yl,legend1,legend2,titu,simple,xlab,addingcluster,timesignclust,coloresclus,rdline,rdlinecolor,legend3,showleg)
  2. for pl=1:2;
  3. dia_acc=a{pl}
  4. timetoplot=time(150:end);
  5. mean_da=nanmean(dia_acc(150:end,:),2);
  6. sem_da=nanstd(dia_acc(150:end,:),0,2)/sqrt(sum(~isnan(dia_acc(1,:))));
  7. line([timetoplot(1),timetoplot(end)], [0,0], 'Color', 'k','LineStyle','-','LineWidth',2);hold on;
  8. % line([timetoplot(1),timetoplot(end)], [chancelevel,chancelevel], 'Color', 'k','LineStyle','--');hold on;
  9. a_plot(pl)=plot(timetoplot,mean_da,'Color',colores(pl,:),'LineWidth',3);hold on;
  10. 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);
  11. a_plot(20)=plot(time,rdline,'--','LineWidth',2,'Color',rdlinecolor);hold on;
  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','-.');hold on;
  34. % a_plot(7)=line([2000,2000],[-10 10], 'Color', 'k','LineStyle','--');hold on;
  35. a_plot(8)=fill([0 350 0 350],[20 20 0 0],'k','LineStyle','none','FaceAlpha',0.2);
  36. a_plot(9)=fill([0 350 0 350],[20 20 0 0],'k','LineStyle','none','FaceAlpha',0.0);
  37. a_plot(10)=fill([0 350 0 350],[20 20 0 0],'k','LineStyle','none','FaceAlpha',0.0);
  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