wolffd@0: function marginal = marginal_family(engine, i, add_ev)
wolffd@0: % MARGINAL_FAMILY Compute the marginal on the specified family (jtree)
wolffd@0: % marginal = marginal_family(engine, i)
wolffd@0: 
wolffd@0: if nargin < 3, add_ev = 0; end
wolffd@0: assert(~add_ev);
wolffd@0: 
wolffd@0: bnet = bnet_from_engine(engine);
wolffd@0: fam = family(bnet.dag, i);
wolffd@0: c = engine.clq_ass_to_node(i);
wolffd@0: marginal = pot_to_marginal(marginalize_pot(engine.clpot{c}, fam));