Home > chronux_1_50 > spikesort > correlations.m

correlations

PURPOSE ^

CORRELATIONS another temporary utility to show clusters

SYNOPSIS ^

function correlations(spikes, useassigns, show);

DESCRIPTION ^

    CORRELATIONS  another temporary utility to show clusters
       CORRELATIONS(SPIKES, [USEASSIGNS], [SHOW]);

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function correlations(spikes, useassigns, show);
0002 %    CORRELATIONS  another temporary utility to show clusters
0003 %       CORRELATIONS(SPIKES, [USEASSIGNS], [SHOW]);
0004 
0005 %   Last Modified By: sbm on Fri Sep 16 07:03:16 2005
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;  % (msec);
0020 
0021 for r = 1:K
0022     for c = r:K  % upper right triangle only (b/c symmetric)
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;  % blank out autocorr peak
0029             bar(lags,cross,1.0);  shading flat;
0030             set(gca, 'XLim', [-maxlag, maxlag]);
0031         else
0032             cla;  % show blank if <= 1 pts
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

Generated on Mon 09-Oct-2006 00:54:52 by m2html © 2003