function figureOut=PBMultiPanelPlot(out) figure; Con=out.connect; blurredImage=out.blurredImage; P=out.P; molaspass=interp1([1 31 82 133 184 256],[0 0 0; 0 0 .75; .5 0 .8; 1 .1 0; 1 .9 0; 1 1 1],1:256); blanding=interp1([1 128 129 256],[0 1 1; .248 .2559 .2559; .2559 .248 .248; 1 0 0],1:256); blanding=flipud(blanding); subplot(3,3,1); Con2=(Con-min(Con(:)))/(max(Con(:))-min(Con(:))); Con2=(Con2-mode(Con2(:)))*256+128; image(Con2) colormap(blanding) freezeColors subplot(3,3,4:9); imagesc(blurredImage'); colormap(molaspass); hold on classBreaks=[8.5 16.5 24.5 32.5 41.5 50.5]; for i=1:length(classBreaks) plot([0 P.N], classBreaks(i)*ones(1,2),'w-'); end xlim([0 P.N]); subplot(3,3,1:3); hold on; tick=0; for i=1:size(out.inputPSCs,2); temp=out.inputPSCs(:,i); % temp=temp(2:end)-temp(1:end-1); temp=find(temp==1); scatter(temp,zeros(size(temp))-i,'k.'); end % temp=out.Iin'>6*10^-8; % blurWidth=30*(1e-3/out.P.dt); % temp=imfilter(double(temp([33:40 43:50],:)),fspecial('gaussian',[1,blurWidth],blurWidth/5)); % imagesc(temp); set(gca,'Xtick',[]); set(gca,'Ytick',[]);