Upload files to 'code'

code/code/Fig_1/Fig_1.m

@@ -0,0 +1,39 @@
+load('data_fig_1.db','-mat');
+figure;
+
+%% theta frequency
+frwin=3;
+for a=1:2
+    rx=size(data.unwrapped_phase_right{a},1);
+    freq_right{a}=([data.unwrapped_phase_right{a}(:,1+frwin:end) nan(rx,frwin)]-[nan(rx,frwin) data.unwrapped_phase_right{a}(:,1:end-frwin)])/2/pi/0.02/frwin/2;
+    rx=size(data.unwrapped_phase_left{a},1);
+    freq_left{a}=([data.unwrapped_phase_left{a}(:,1+frwin:end) nan(rx,frwin)]-[nan(rx,frwin) data.unwrapped_phase_left{a}(:,1:end-frwin)])/2/pi/0.02/frwin/2;    
+end
+
+%% acceleration
+win=10;    %gfr=fr;
+
+
+    rx=size(data.x_right,1);
+    acc_right=-([nan(rx,win) data.x_right(:,1:end-win)]+[data.x_right(:,1+win:end) nan(rx,win)]-2*data.x_right)/0.02/0.02/win/win;
+    rx=size(data.x_left,1);    
+    acc_left=([nan(rx,win) data.x_left(:,1:end-win)]+[data.x_left(:,1+win:end) nan(rx,win)]-2*data.x_left)/0.02/0.02/win/win;    
+
+
+ fct=0; 
+th=8; axlim=[-120  220]; 
+subplot(2,1,1)
+
+gacc=acc_right;
+gfr=freq_right;
+do_plots;
+
+set(AX(2),'XTick',[]);
+
+subplot(2,1,2)
+    gfr=freq_left; 
+    gacc=acc_left;    
+
+do_plots;
+set(AX(2),'XTick',[]);
+

code/code/Fig_1/Fig_1.m~

@@ -0,0 +1,48 @@
+load('data_fig_1.db','-mat');
+
+
+subplot(2,1,1)
+
+%% theta frequency
+frwin=3;
+for a=1:3
+    rx=size(data.unwrapped_phase_right{a},1);
+    freq_right{a}=([data.unwrapped_phase_right{a}(:,1+frwin:end) nan(rx,frwin)]-[nan(rx,frwin) data.unwrapped_phase_right{a}(:,1:end-frwin)])/2/pi/0.02/frwin/2;
+    rx=size(data.unwrapped_phase_left{a},1);
+    freq_left{a}=([data.unwrapped_phase_left{a}(:,1+frwin:end) nan(rx,frwin)]-[nan(rx,frwin) data.unwrapped_phase_left{a}(:,1:end-frwin)])/2/pi/0.02/frwin/2;    
+end
+
+%% acceleration
+win=10;    %gfr=fr;
+
+
+    rx=size(x,1);
+    acc_right=-([nan(rx,win) x(:,1:end-win)]+[x(:,1+win:end) nan(rx,win)]-2*x)/0.02/0.02/win/win;
+
+
+axlim=0; fct=0; th=0; 
+th=ths(r); axlim=[-120  220]; %fct=0.0127;
+new_figfig
+% set(AX(2),'XTickLabel',[]);
+% set(AX(1),'XTickLabel',[]);
+
+set(AX(2),'XTick',[]);
+% L=legend('acceleration','CA1','MEC','Location', 'NorthEast');
+% set(L,'Color','none')
+% set(L,'box','off')
+% set(L,'FontSize',10,'FontName','Arial')
+subplot(2,1,2)
+% win=5;
+if r==3
+    gfr=data.unwrapped_phase_left; gacc=lacc;
+    rx=size(lx,1);
+    gacc=([nan(rx,win) lx(:,1:end-win)]+[lx(:,1+win:end) nan(rx,win)]-2*lx)/0.02/0.02/win/win;    
+else
+    gfr=fr; gacc=acc; 
+    rx=size(x,1);
+    gacc=-([nan(rx,win) x(:,1:end-win)]+[x(:,1+win:end) nan(rx,win)]-2*x)/0.02/0.02/win/win;    
+end
+% axlim=[-160 260];
+new_figfig
+set(AX(2),'XTick',[]);
+

code/code/Fig_1/do_plots.m

@@ -0,0 +1,62 @@
+
+% selx=423:624;
+selx=451:601;
+bx=(-500:500)/50;
+col1=.1*[1 1 1];%[.0 .6 .6];
+col2=[255 175 3]/255;%[77 255 255]/255;%[255 195 77]/255;%[225 54 255]/255;%[255 110 238]/255;[0 .63 .54];
+col3=[255 0 0]/255;%[106 77 255]/255;%[77 255 106]/255;%[.60 .65 .0];%[77 255 106]/255;%[0 1 1]*.9;
+xxx=gacc(:,selx);
+mxxx=nanmean(xxx);
+sxxx=nanstd(xxx)/2;
+hold on;
+
+
+[AX,H1,H2]=plotyy(bx(selx),nanmean(gacc(:,selx)),bx(selx),[nanmean(gfr{1}(:,selx))' nanmean(gfr{2}(:,selx))']);
+uistack(AX(1), 'top')
+
+set(AX(1),'Color','none')
+set(AX(1),'FontName','Arial');
+set(AX(2),'FontName','Arial');
+if th==0
+    th=nanmean(nanmean(gfr{2}(:,[423:450 600:624])));
+    th=.5*round(th/.5);
+end
+if fct==0
+    fct= (nanmean([max(nanmean(gfr{2}(:,selx))) max(nanmean(gfr{1}(:,selx)))])-th)/max(mxxx);
+end
+set(H1,'LineWidth',1);set(H2,'LineWidth',2);
+set(H1,'LineStyle','--','Color',col1);
+ set(H2(1),'Color',col2)
+set(H2(2),'Color',col3)
+if axlim==0
+    axlim=1.01*[1.15*min(mxxx-sxxx) max(mxxx+sxxx)];
+end
+axs=-1000:100:1000;
+set(AX(1),'YLim',axlim);
+set(AX(1),'YLim',axlim); %axlim=get(AX(2),'YLim');
+set(AX(1),'YTick',axs);
+set(AX(1),'XTick',-3:1:3);
+set(AX(2),'YLim',axlim*fct+th-0.0);%[6.332 10]);
+set(AX(2),'YTick',6:11);
+set(AX(1),'XLim',[-1 2]);
+set(AX(2),'XLim',[-1 2]);
+set(AX(1),'box','off');
+set(AX(1),'YColor','k');
+set(AX(2),'YColor',col3);
+AX(1)=reg_plots(AX(1));
+AX(2)=reg_plots(AX(2));
+
+
+
+AX_annotations = copyobj(AX(1),get(AX(1),'Parent'));
+A=area(AX_annotations,bx(selx),[mxxx-sxxx; 2*sxxx]'); 
+set(A(1),'Visible','off');
+set(A(2),'FaceColor',col1+(1-col1)*.8);
+set(A(2),'EdgeColor','none');
+uistack(AX_annotations, 'bottom');
+set(AX_annotations,'Color',[1 1 1]);
+set(AX_annotations,'XTick',[],'YTick',[]);
+dsp=diff(axlim)/200*120;
+ set(AX(1),'dataaspectratio',[1 dsp 1]);
+ set(AX(2),'dataaspectratio',[1 dsp*fct 1]);
+ set(AX_annotations,'dataaspectratio',[1 dsp 1]);

code/code/Fig_1/reg_plots.m

@@ -0,0 +1,3 @@
+function ax=reg_plots(ax)
+set(ax,'FontSize',13);
+set(ax,'FontName','Arial');