function [avgPI,stdPI,avgPE,stdPE] = calculatePauses(newP,newT);
% newP contains the markers for the start and stop of each pause around
% each good peak - DECIMATED. newP(1) is the start of the pause for
% peak #1, newP(2) is the stop for the pause around peak #1, etc.
% newT is the same for the troughs.
% the average and standard deviation of the post-inspiratory pauses (PI) and
% post-expiratory pauses (PE) are returned.
allDiffs = diff(newP); % finds the diffs between ALL adjacent points in newP
Pdiffs = []; % initialize empty array to hold the true differences
n = max(size(allDiffs));
for i = 1:n
if mod(i,2) > 0 % save all ODD indexed values (1,3,5, ...) into Pdiffs
Pdiffs = [Pdiffs, (allDiffs(i) * 5)]; % multiply by 5 for undecimated size
end;
end;
Pdiffs = abs(Pdiffs); % so all values are positive
allTDiffs = diff(newT); % finds the diffs between ALL adjacent points in newT
Tdiffs = []; % initialize empty array to hold the true differences
n = max(size(allTDiffs));
for i = 1:n
if mod(i,2) > 0 % save all ODD indexed values (1,3,5, ...) into Pdiffs
Tdiffs = [Tdiffs, (allTDiffs(i) * 5)]; % multiply by 5 for undecimated size
end;
end;
Tdiffs = abs(Tdiffs); % so all values are positive
avgPI = mean(Pdiffs); avgPE = mean(Tdiffs); % averages
stdPI = std(Pdiffs); stdPE = std(Tdiffs); % standard deviations