Noninvasive Fetal ECG: The PhysioNet/Computing in Cardiology Challenge 2013 1.0.0

function [Tq std_rr iqr_rr med_rr ampl] = MyBuildSeq_3(S, Fs, bs)
%
% function: attempts to build sequence of annotatios
%

v_mu    = 350:10:480;
v_sigma = 10:5:50;
rr_iqr  = zeros(numel(v_mu)*numel(v_sigma), 1); 
rr_med  = zeros(numel(v_mu)*numel(v_sigma), 1); 
err     = zeros(numel(v_mu)*numel(v_sigma), 1); 
ampl    = zeros(numel(v_mu)*numel(v_sigma), 1); 
v       = zeros(numel(v_mu)*numel(v_sigma), 2); 
j = 1;
for mj=1:numel(v_mu)
    for sj=1:numel(v_sigma)
        Tq = BuildChain(S, Fs, v_mu(mj), v_sigma(sj));
        if( isempty(Tq)==0 )
            rr        = diff(Tq);
            rr_iqr(j) = iqr(rr);
            rr_med(j) = median(rr);
            
            err(j)    = std(rr); 
            ampl(j)   = std(S(Tq));
            v(j,1)    = mj;
            v(j,2)    = sj;
            j = j+1; 
        end
    end
end

rr_iqr(j:end) = [];
rr_med(j:end) = [];
err(j:end)    = [];
ampl(j:end)   = [];
v(j:end, :)   = [];

[~, I]  = min(err);
sigma       = v_sigma( v(I,2) );
mu          = v_mu( v(I,1) );
Tq          = BuildChain(S, Fs, mu, sigma);

std_rr = err(I);
iqr_rr = rr_iqr(I);
med_rr = rr_med(I);
ampl   = ampl(I);