You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
74 lines
2.0 KiB
Matlab
74 lines
2.0 KiB
Matlab
function plot_reverse_correl_freq(w,rcFreq,bin_labels,plotNodes,plotTimes,subBase)
|
|
colPlus=1;
|
|
if nargin<4
|
|
numNodes=length(rcFreq);
|
|
else
|
|
numNodes=length(plotNodes);
|
|
end
|
|
[poo,numTimes]=size(rcFreq{1});
|
|
maxNodes=size(w,1);
|
|
maxTime=floor(numTimes/2);
|
|
if nargin<5
|
|
plotTimes=[-maxTime:maxTime];
|
|
else
|
|
numTimes=length(plotTimes);
|
|
end
|
|
if nargin<6, subBase=0; end
|
|
|
|
numBins=size(rcFreq{1},1);
|
|
%binsToLabel=[1:floor(numBins/2):numBins+1];
|
|
binsToLabel=[floor(numBins/4),2+ceil(numBins*0.75)];
|
|
%binsToLabel=[1:numBins]
|
|
|
|
clf
|
|
nodeNum=0;
|
|
for node=plotNodes
|
|
nodeNum=nodeNum+1;
|
|
%for each node show the the reverse correlation at each time lag
|
|
toff=0;
|
|
for time=plotTimes
|
|
toff=toff+1;
|
|
maxsubplot(numNodes+1,colPlus+numTimes,(colPlus+numTimes)*(nodeNum)+toff+colPlus);
|
|
data=rcFreq{node}(:,time+maxTime+1);
|
|
if subBase %remove response of node with same phase but with orientation
|
|
%preference perpendicular to the stimulus
|
|
subNode=8*floor(node/8)+1;
|
|
data=data-rcFreq{subNode}(:,time+maxTime+1);
|
|
end
|
|
data=data./sum(data); %Ringbach et al 1997
|
|
plot(data,'LineWidth',2),
|
|
hold on, plot([0,numBins+5],[0,0],'k')
|
|
%if subBase, axis([0,numBins+3,-0.00001,0.000175]); else,
|
|
%axis([0,numBins+3,0,0.001]); end
|
|
if subBase, axis([0,numBins+3,-0.01,0.2]); else, axis([0,numBins+3,0,0.1]);
|
|
end
|
|
drawnow;
|
|
if nodeNum==1 & ~isempty(bin_labels)
|
|
set(gca,'XTick',binsToLabel,'XTickLabel',num2str(bin_labels(binsToLabel)',2));
|
|
else
|
|
set(gca,'XTick',[]);
|
|
end
|
|
if toff==1
|
|
set(gca,'YTick',[0:0.05:0.15]);
|
|
ylabel('p(\tau,f) / \Sigma_f p(\tau,f)','FontSize',12)
|
|
else
|
|
set(gca,'YTick',[]);
|
|
end
|
|
if nodeNum==1 & toff==ceil(numTimes/2)
|
|
xlabel('Spatial Frequency (cycles/pixel)','FontSize',12)
|
|
end
|
|
set(gca,'FontSize',10);
|
|
end
|
|
|
|
if node<=maxNodes
|
|
%show corresponding synaptic weigths
|
|
RF=w{node,1}-w{node,2};
|
|
lim=max(1e-9,0.9.*max(max(abs(RF))));
|
|
maxsubplot(numNodes+1,2+numTimes,(numTimes+2)*(nodeNum)+1);
|
|
imagesc(RF,[-lim,lim]);
|
|
axis('equal','tight'), set(gca,'XTick',[],'YTick',[]); drawnow;
|
|
end
|
|
end
|
|
|
|
set(gcf,'PaperPosition',2.*[1 1 numTimes 2*numNodes]);
|