https://github.com/PerezOrtegaJ/Neural_Ensemble_Analysis
Tip revision: 9d37fd031dfbdb4eb69faa449d0a6416267a7d4f authored by Jesús Pérez on 28 July 2020, 20:36:58 UTC
Update README.md
Update README.md
Tip revision: 9d37fd0
Get_Ensembles_Length.m
function [widths,nPeaks] = Get_Ensembles_Length(indices,sequence)
% Identify the legnth of the ensembles and plot the distribution of each one
%
% widths = Plot_Ensembles_Length(indices,sequence)
%
% By Jesus Perez-Ortega, Apr 2020
% Find vectors
id = find(indices>0);
% Get the ensemble id
ensembles = unique(sequence)';
nEns = length(ensembles);
% Get colors
for i = 1:nEns
seqEns = [];
% Create binary signal to identify the lenght of each activation
seqEns(id(sequence==ensembles(i))) = 1;
% find peaks
[~,w] = Find_Peaks_Or_Valleys(seqEns);
% get number of peaks
nPeaks(i) = length(w);
% assign widths
widths{ensembles(i)} = w;
end