wolffd@0: function m = marginal_family(engine, b, i, t, add_ev)
wolffd@0: % MARGINAL_FAMILY Compute the marginal on the specified family (jtree_2TBN)
wolffd@0: % marginal = marginal_family(engine, b, i, t, add_ev)
wolffd@0: 
wolffd@0: bnet = bnet_from_engine(engine);
wolffd@0: if t==1
wolffd@0:   m = marginal_nodes(engine, b, family(bnet.dag, i), t, add_ev, 1);
wolffd@0: else
wolffd@0:   ss = length(bnet.intra);
wolffd@0:   fam = family(bnet.dag, i+ss);
wolffd@0:   m = marginal_nodes(engine, b, fam, t, add_ev, 1);
wolffd@0: end