%% plot correlations of N/c model (Fig. 6) clear all; close all % load data for all four trials load Data_Fig6 Cavg=squeeze(nanmean(CC,2)); Csem=squeeze(nanstd(CC,[],2)./sqrt(size(CC,2))); ll=size(Cavg,1); set(0,'DefaultAxesFontSize',12) fig=figure; [ha, pos] = tight_subplot(3,5,.03,.1,.15); for pcount=1:ll axes(ha(pcount)); %h = fill([lags,fliplr(lags)], [Cavg(pcount,:)+Csem(pcount,:),fliplr(Cavg(pcount,:)-Csem(pcount,:))], 'k','LineStyle','none'); %set(h,'facealpha',.3) ylim([-.25 .25]); hold on; plot(lags,Cavg(pcount,:),'k') axis=([-100 100 -0.2 0.2]); set(gca,'Clipping','off','Box','off','TickDir','out','Color','none'); set(gca, 'XTick', [-50,0,50]) if pcount<6 str=['c = ' num2str((pcount-1)*.1)]; text(.4,1.2,str,'Units', 'Normalized', 'VerticalAlignment', 'Top','FontSize',12) end if pcount==1 text(-1.2,.5,'N=5','Units', 'Normalized', 'VerticalAlignment', 'Top','FontSize',12) end if pcount==6 text(-1.2,.5,'N=20','Units', 'Normalized', 'VerticalAlignment', 'Top','FontSize',12) end if pcount==11 text(-1.2,.5,'N=100','Units', 'Normalized', 'VerticalAlignment', 'Top','FontSize',12) end if pcount==13 text(.1,-.25,'time lag (ms)','Units', 'Normalized', 'VerticalAlignment', 'Top','FontSize',12) end if pcount==6 h=text(-.6,0,'GPi-VLa correlation','Units', 'Normalized', 'VerticalAlignment', 'Top','FontSize',12) set(h,'Rotation',90); end end set(ha(1:10),'XTickLabel',''); set(ha([2:5,7:10,12:15]),'YTickLabel','') set(ha(11:15),'XTickLabel',{'-50','0','50'}); %ylabel('GPi-VLa correlation') %xlabel('time lag (ms)') % average effect (0-10ms) CCdip=nanmean(CC(:,:,100:110),3); CCdipmean=nanmean(CCdip,2); % CCdipsem=nanstd(CCdip,[],2); CCdipsem=nanstd(CCdip,[],2)./sqrt(size(CCdip,2)); % average 5 CCFs of each trial CCdiptrial1=nanmean(CCdip(:,1:5),2); CCdiptrial2=nanmean(CCdip(:,6:10),2); CCdiptrial3=nanmean(CCdip(:,11:15),2); CCdiptrial4=nanmean(CCdip(:,16:20),2); CCdipN5=[CCdiptrial1(1:5),CCdiptrial2(1:5),CCdiptrial3(1:5),CCdiptrial4(1:5)]; CCdipN20=[CCdiptrial1(6:10),CCdiptrial2(6:10),CCdiptrial3(6:10),CCdiptrial4(6:10)]; CCdipN100=[CCdiptrial1(11:15),CCdiptrial2(11:15),CCdiptrial3(11:15),CCdiptrial4(11:15)]; ct=[0:.1:.4]'; % approximation -(c(N-1)/N-1/N)*k; % determine k by least squares fit k=[.01:.001:0.4]; N5=CCdipmean(1:5); N20=CCdipmean(6:10); N100=CCdipmean(11:15); yN5=@(k) -k*4/5*ct - k/5; yN20=@(k) -k*19/20*ct - k/20; yN100=@(k) -k*99/100*ct - k/100; for i=1:length(k) dist5(i,:)=(N5-yN5(k(i))).^2; dist20(i,:)=(N20-yN20(k(i))).^2; dist100(i,:)=(N100-yN100(k(i))).^2; end dist=[dist5,dist20,dist100]; distsum=sum(dist,2); %plot(k,distsum) [min,minind]=min(distsum); kmin=k(minind); y5=-kmin*4/5*ct - kmin/5; y20=-kmin*19/20*ct - kmin/20; y100=-kmin*99/100*ct - kmin/100; set(0,'DefaultAxesFontSize',12) fig=figure; subplot(2,2,3:4) hold on plot(ct,y5,'--','Color',[.7 .2 .2]) plot(ct,y20,'--','Color',[.2 .7 .2]) plot(ct,y100,'--','Color',[.2 .2 .7]) errorbar(ct,CCdipmean(1:5),CCdipsem(1:5),'Color',[.7 .2 .2],'LineWidth',1) errorbar(ct,CCdipmean(6:10),CCdipsem(6:10),'Color',[.2 .7 .2],'LineWidth',1) errorbar(ct,CCdipmean(11:15),CCdipsem(11:15),'Color',[.2 .2 .7],'LineWidth',1) text(.8,.9,'N=5','Units', 'Normalized', 'VerticalAlignment', 'Top','Color',[.7 .2 .2],'FontSize',12) text(.8,.8,'N=20','Units', 'Normalized', 'VerticalAlignment', 'Top','Color',[.2 .7 .2],'FontSize',12) text(.8,.7,'N=100','Units', 'Normalized', 'VerticalAlignment', 'Top','Color',[.2 .2 .7],'FontSize',12) xlim([-.01 .41]) xlabel('GPi-GPi correlation (c)'); ylabel({'GPi-VLa correlation' '(average [0,10]ms)'}) set(gca,'Box','off','TickDir','out','Color','none');