Daniel@0: function marginal = marginal_family(engine, i, t, add_ev)
Daniel@0: % MARGINAL_FAMILY Compute the joint distribution on a set of family (smoother_engine)
Daniel@0: % function marginal = marginal_family(engine, i, t, add_ev)
Daniel@0: 
Daniel@0: if nargin < 4, add_ev = 0; end
Daniel@0: marginal = marginal_family(engine.tbn_engine, engine.b{t}, i, t, add_ev);