"Less than 15 spikes detected" error message

[Hugues1273]

Dear all,

I tried to make spike sorting using wave_clus. I have my continuous data in .mat files. When I load a file with the data, I get this message error:

Error using wave_clus>load_data_button_Callback (line 230)
Less than 15 spikes detected

Error in gui_mainfcn (line 95)
        feval(varargin{:});

Error in wave_clus (line 63)
    gui_mainfcn(gui_State, varargin{:});

Error while evaluating UIControl Callback.

I have already find an issue with a problem close to the mine. The advice was to modify the parameter file. I did it but I still have this error message. Here is my file parameter:

% LOAD PARAMS                                                                                                                               
par.segments_length = 5;             % length (in minutes) of segments in which the data is cutted (default 5min).                          
par.sr = 48000;                      % sampling rate (in Hz). This parameter will be only used if the data file don't have a sr.            

% PLOTTING PARAMETERS                                                                                                                       
par.cont_segment = true;                                                                                                                    
par.max_spikes_plot = 1000;          % max. number of spikes to be plotted                                                                  
par.print2file = true;               % If is not true, print the figure (only for batch scripts).                                           
par.cont_plot_samples = 100000;      % number of samples used in the one-minute (maximum) sample of continuous data to plot.                
par.to_plot_std = 1;                 % # of std from mean to plot                                                                           
par.all_classes_ax = 'mean';         % 'mean'/'all'. If it's 'mean' only the mean waveforms will be ploted in the axes with all the classes 
par.plot_feature_stats = false;                                                                                                             

% SPC PARAMETERS                                                                                                                            
par.mintemp = 0.00;                  % minimum temperature for SPC                                                                          
par.maxtemp = 0.251;                 % maximum temperature for SPC                                                                          
par.tempstep = 0.01;                 % temperature steps                                                                                    
par.SWCycles = 100;                  % SPC iterations for each temperature (default 100)                                                    
par.KNearNeighb = 11;                % number of nearest neighbors for SPC                                                                  
par.min_clus = 20;                   % minimum size of a cluster (default 60)                                                               
par.randomseed = 0;                  % if 0, random seed is taken as the clock value (default 0)                                            
%par.randomseed = 147;               % If not 0, random seed                                                                                
%par.temp_plot = 'lin';              % temperature plot in linear scale                                                                     
par.temp_plot = 'log';               % temperature plot in log scale                                                                        

par.c_ov = 0.7;                      % Overlapping coefficient to use for the inclusion criterion.                                          
par.elbow_min  = 0.4;                %Thr_border parameter for regime border detection.                                                     

% DETECTION PARAMETERS                                                                                                                      
par.tmax = 'all';                    % maximum time to load                                                                                 
%par.tmax= 180;                      % maximum time to load (in sec)                                                                        
par.tmin= 0;                         % starting time for loading (in sec)                                                                   
par.w_pre = 20;                      % number of pre-event data points stored (default 20)                                                  
par.w_post = 44;                     % number of post-event data points stored (default 44))                                                
par.alignment_window = 10;           % number of points around the sample expected to be the maximum                                        
par.stdmin = 1;                      % minimum threshold for detection                                                                      
par.stdmax = 500;                     % maximum threshold for detection                                                                     
par.detect_fmin = 30;               % high pass filter for detection                                                                        
par.detect_fmax = 3000;              % low pass filter for detection (default 1000)                                                         
par.detect_order = 5;                % filter order for detection. 0 to disable the filter.                                                 
par.sort_fmin = 30;                 % high pass filter for sorting                                                                          
par.sort_fmax = 3000;                % low pass filter for sorting (default 3000)                                                           
par.sort_order = 2;                  % filter order for sorting. 0 to disable the filter.                                                   
par.ref_ms = 1.5;                    % detector dead time, minimum refractory period (in ms)                                                
par.detection = 'both';               % type of threshold ('pos','neg','both')                                                              
% par.detection = 'neg';                                                                                                                    
% par.detection = 'pos';                                                                                                                    

% INTERPOLATION PARAMETERS                                                                                                                  
par.int_factor = 5;                  % interpolation factor                                                                                 
par.interpolation = 'y';             % interpolation with cubic splines (default)                                                           
% par.interpolation = 'n';                                                                                                                  

% FEATURES PARAMETERS                                                                                                                       
par.min_inputs = 10;         % number of inputs to the clustering                                                                           
par.max_inputs = 0.75;       % number of inputs to the clustering. if < 1 it will the that proportion of the maximum.                       
par.scales = 4;                        % number of scales for the wavelet decomposition                                                     
par.features = 'wav';                % type of feature ('wav' or 'pca')                                                                     
%par.features = 'pca'                                                                                                                       

% FORCE MEMBERSHIP PARAMETERS                                                                                                               
par.template_sdnum = 3;             % max radius of cluster in std devs.                                                                    
par.template_k = 10;                % # of nearest neighbors                                                                                
par.template_k_min = 10;            % min # of nn for vote                                                                                  
%par.template_type = 'mahal';       % nn, center, ml, mahal                                                                                 
par.template_type = 'center';       % nn, center, ml, mahal                                                                                 
par.force_feature = 'spk';          % feature use for forcing (whole spike shape)                                                           
%par.force_feature = 'wav';         % feature use for forcing (wavelet coefficients).                                                       
par.force_auto = true;              %automatically force membership (only for batch scripts).                                               

% TEMPLATE MATCHING                                                                                                                         
par.match = 'y';                    % for template matching                                                                                 
%par.match = 'n';                   % for no template matching                                                                              
par.max_spk = 40000;                % max. # of spikes before starting templ. match.                                                        
par.permut = 'y';                   % for selection of random 'par.max_spk' spikes before starting templ. match.                            
% par.permut = 'n';                 % for selection of the first 'par.max_spk' spikes before starting templ. match.                         

% HISTOGRAM PARAMETERS                                                                                                                      
par.nbins = 100;                    % # of bins for the ISI histograms                                                                      
par.bin_step = 1;                   % percentage number of bins to plot    

I know there is more than 15 spikes in my data (much more!) so I am really desperated... :(

If someone has a solution I would be so grateful!!

Thank you so much

[ferchaure]

Hi, par.stdmin=1 is too low (It can detect a lot of spikes but discard them because of the refractory period). Try with par.stdmin=3.5 Maybe an issue with the input file? Can you plot a segment of the data?

[Hugues1273]

Thank you so much for your answer! You are very reactive! thank you!

I put 3.5 but it did not change anything. I still have the same error message :(

I put you a 1 second segment of my record get with spyking-circus. is that a way to do that with wave clus?

[Hugues1273]

Maybe I could send you a file to try?

[ferchaure]

The data looks fine,

If you want, sent me the data to : fchaure at fi.uba.ar