find_carryover_significance_zerolag.m

function Otot=find_carryover_significance_zerolag(Otot)
% this checks whether the increment of the O information going one order up
% is significant with respect to the value at the lower order

maxsize=length(Otot);
for isize=maxsize:-1:4
    higher=isize;
    nlower_red=size(Otot(higher-1).index_var_red,1);
    nlower_syn=size(Otot(higher-1).index_var_syn,1);
    nhigher_red=size(Otot(higher).index_var_red,1);
    nhigher_syn=size(Otot(higher).index_var_syn,1);
    sig_red=zeros(nhigher_red,1);
    sig_syn=zeros(nhigher_syn,1);
    for ilower=1:nlower_red
        for ihigher=1:nhigher_red
            if sum(ismember(Otot(higher-1).index_var_red(ilower,:),Otot(higher).index_var_red(ihigher,:)))==higher-1
                v1=Otot(higher-1).bootsigCI_red(ilower,:);
                v2=Otot(higher).bootsigCI_red(ihigher,:);
                l1=v1(1);l2=v2(1);r1=v1(2);r2=v2(2);
                r = bsxfun(@min, r1, r2);
                l = bsxfun(@max, l1, l2);
                sig_red(ihigher)= l >= r;
            end
        end
    end
    Otot(isize).inc_sig_red=sig_red;
    for ilower=1:nlower_syn
        for ihigher=1:nhigher_syn
            if sum(ismember(Otot(higher-1).index_var_syn(ilower,:),Otot(higher).index_var_syn(ihigher,:)))==higher-1
                v1=Otot(higher-1).bootsigCI_syn(ilower,:);
                v2=Otot(higher).bootsigCI_syn(ihigher,:);
                l1=v1(1);l2=v2(1);r1=v1(2);r2=v2(2);
                r = bsxfun(@min, r1, r2);
                l = bsxfun(@max, l1, l2);
                sig_syn(ihigher)= l >= r;
            end
        end
    end
    Otot(isize).inc_sig_syn=sig_syn;
end