function [isi_R, isi_T, T_cens] = EM_load(full_filename)

% -------------------------------------------------
% Load the dataset

dataset = load(full_filename);


% -------------------------------------------------
% Split spike data in regular/truncated isi if necessary

if isnumeric(dataset) % then assume the dataset is a collection of (uncensored) isi
    
    isi_R = dataset;
    
    isi_T = [];
    
    T_cens = Inf;

        
else % if datafile is *.mat, assume it contains spike times cell array TS and trial_length
    

    T_cens = dataset.trial_length;
    
    [isi_R, isi_T] = EM_extract(dataset.TS, T_cens);
            
        
end

N_R = length(isi_R);