Skip to main content

Browse the archive

Revision fec1958a3f9502fd8e3054cb0b95f55fc77ba527 authored by Max Turner on 12 April 2016, 00:46:24 UTC, committed by Max Turner on 12 April 2016, 00:46:24 UTC 
More pacakge name changes for riekelab name convention
1 parent ce75ca2
Raw File
getRebounds.m 
function r = getRebounds(peaks_ind,trace,searchInterval)
%get rebound as fraction of peak amplitude

%trace = abs(trace);
peaks = trace(peaks_ind);
r = zeros(size(peaks));

for i=1:length(peaks)
   endPoint = min(peaks_ind(i)+searchInterval,length(trace));
   nextMin = getPeaks(trace(peaks_ind(i):endPoint),-1);
   if isempty(nextMin), nextMin = peaks(i); 
   else nextMin = nextMin(1); end
   nextMax = getPeaks(trace(peaks_ind(i):endPoint),1);
   if isempty(nextMax), nextMax = 0; 
   else nextMax = nextMax(1); end
   
   if nextMin<peaks(i) %not the real spike min
       r(i) = 0;
   else
       r(i) = nextMax; 
   end
end
back to top