construct_equiv_hier_hmm.m

The BNL toolbox is a set of Matlab functions for defining and estimating the parameters of a Bayesi

%equiv_class: vector. (ordered) nodes having the same CPT have same value
%equiv_class_time: vector.  (ordered) nodes having a CPT governed by the
    %same set of parameters have the same value. CPTs may be different (eg when because
    %different values on covariates)
%lin_pred_struct: cell each element gives the structure of the linear
%predictor for each equiv_class of nodes
equiv_class=zeros(nr_nodes,1);
equiv_class(startdaynode)=1;
equiv_class(mysetdiff(hid_day_nodes,startdaynode))=2;
equiv_class(startsignalnodes)=3;
equiv_class(terminal_merged_nodes.nodenrs)=4;
equiv_class(mysetdiff(hid_signal_nodes,startsignalnodes))=5;

%%define equiv_class_time
equiv_class_time=equiv_class;
f5=find(equiv_class==5);
equiv_class(f5)=5:length(f5)+4;