utall-Strand MVAR estimation algorithm</span>
criterion = 1;      <span class="comment">% 1: AIC, Akaike Information Criteria</span>

disp(<span class="string">'Running MVAR estimation routine.'</span>)

[IP,pf,A,pb,B,ef,eb,vaic,Vaicv] = mvar(u,maxIP,alg,criterion);

disp([<span class="string">'Number of channels = '</span> int2str(nChannels) <span class="string">' with '</span> <span class="keyword">...</span>
    int2str(nSegLength) <span class="string">' data points; MAR model order = '</span> int2str(IP