0001 function correlations(spikes, useassigns, show);
0002
0003
0004
0005
0006
0007 if (nargin < 2),
0008 if (isfield(spikes.hierarchy, 'assigns')), useassigns = spikes.hierarchy.assigns;
0009 elseif (isfield(spikes.overcluster, 'assigns')), useassigns = spikes.overcluster.assigns;
0010 else useassigns = ones(size(spikes.waveforms,1),1);
0011 end
0012 end
0013 if (nargin < 3)
0014 list = unique(useassigns);
0015 show = list(1:min(5,length(list)));
0016 end
0017
0018 K = length(show);
0019 maxlag = 0.050;
0020
0021 for r = 1:K
0022 for c = r:K
0023 subplot(K,K,c+(r-1)*K);
0024 selectrow = find(useassigns == show(r));
0025 selectcol = find(useassigns == show(c));
0026 if ((length(selectrow) > 1) && (length(selectcol) > 1))
0027 [cross,lags] = pxcorr(spikes.spiketimes(selectrow), spikes.spiketimes(selectcol), 1000, maxlag);
0028 if (r == c), cross(find(lags == 0)) = 0; end;
0029 bar(lags,cross,1.0); shading flat;
0030 set(gca, 'XLim', [-maxlag, maxlag]);
0031 else
0032 cla;
0033 end
0034 if (r == c), ylabel(sprintf('Cluster #%d', show(r))); end;
0035 if (r == 1), title(sprintf('Cluster #%d', show(c))); end;
0036 end
0037 end