wolffd@0: % make the structure of an embedded HMM with 2 rows and 3 columns
wolffd@0: 
wolffd@0: % 1------------>2
wolffd@0: % |\   \        | \  \
wolffd@0: % 3->4->5       6->7->8
wolffd@0: 
wolffd@0: n = 8;
wolffd@0: dag = zeros(n);
wolffd@0: dag(1,[2 3 4 5])=1;
wolffd@0: dag(2,[6 7 8])=1;
wolffd@0: for i=3:4
wolffd@0:   dag(i,i+1)=1;
wolffd@0: end
wolffd@0: for i=6:7
wolffd@0:   dag(i,i+1)=1;
wolffd@0: end
wolffd@0: ns = 2*ones(1,n);
wolffd@0: bnet = mk_bnet(dag,ns);
wolffd@0: for i=1:n
wolffd@0:   bnet.CPD{i}=tabular_CPD(bnet,i);
wolffd@0: end
wolffd@0: [jtree, root, cliques] =  graph_to_jtree(moralize(bnet.dag), ones(1,n), {}, {});
wolffd@0: %[jtree, root, cliques, B, w, elim_order, moral_edges, fill_in_edges] = dag_to_jtree(bnet);
wolffd@0: