0001
0002
0003 load spikes2
0004
0005
0006 spikes = ss_dejitter(spikes);
0007 spikes = ss_outliers(spikes);
0008 spikes = ss_kmeans(spikes);
0009 spikes = ss_energy(spikes);
0010 spikes = ss_aggregate(spikes);
0011
0012
0013
0014
0015
0016
0017
0018 projection = pcasvd(spikes.waveforms);
0019 [assigns,log] = kkm(spikes.waveforms(:,1:3),0);
0020
0021
0022
0023
0024
0025
0026 ssg_databrowse2d(spikes);
0027 ssg_databrowse3d(spikes);
0028
0029
0030
0031 figure; plot(spikes.waveforms'); axis tight; title('Spike Waveforms');
0032 figure; histxt(spikes.waveforms); colormap jet(256);
0033
0034
0035 figure; plot(spikes.outliers.waveforms'); axis tight;
0036
0037
0038 figure; set(gcf, 'Renderer', 'OpenGL'); clusterXT(spikes, spikes.overcluster.assigns); title('Local Clusters');
0039
0040
0041 figure; set(gcf, 'Renderer', 'OpenGL'); clusterXT(spikes, spikes.hierarchy.assigns); title('Final Clusters');
0042 figure; showclust(spikes, spikes.hierarchy.assigns);
0043 figure; correlations(spikes, spikes.hierarchy.assigns);
0044 figure; aggtree(spikes); title('Aggregation Tree');