function [tests,pass,perf]=test_edr(varargin)
inputs={'verbose'};
verbose=0;
for n=1:nargin
if(~isempty(varargin{n}))
eval([inputs{n} '=varargin{n};']);
end
end
%Test the examples
% Note: This requires the Signal Processing toolbox
test_string={['signal=''fantasia/f1o02'';r_peaks=''ecg'';' ...
'data_type=1;channel=2;y=edr(data_type,signal,r_peaks,[],[],[],channel);' ...
'wfdb2mat(''f1o02'');[~,signal,Fs,~]=rdmat(''f1o02m'');' ...
'resp=signal(:,1);resp=resp-mean(resp);resp=resp*200;sec=length(resp)/Fs;' ...
'xax=[.25:.25:sec];r=interp1(y(:,1), y(:,2), xax,''spline'');' ...
'B=fir1(100,[.1 .5],''bandpass'');edr_filt=filtfilt(B,1,r);figure;plot(xax,edr_filt);' ...
'hold on;plot([1:length(resp)]/Fs,resp,''r'');legend(''edr'',''respiratory signal'');' ...
'xlabel(''time (s)'');close all']};
clean_up={['delete([pwd filesep ''f1o02m*'']);']};
[tests,pass,perf]=test_wrapper(test_string,clean_up,verbose);