wolffd@0: function marginal = marginal_family(engine, i, t)
wolffd@0: % MARGINAL_FAMILY Compute the marginal on the specified family (jtree_unrolled_dbn)
wolffd@0: % marginal = marginal_family(engine, i, t)
wolffd@0: 
wolffd@0: if nargin < 3, t = 1; end
wolffd@0: 
wolffd@0: bnet = bnet_from_engine(engine);
wolffd@0: ss = length(bnet.intra);
wolffd@0: marginal = marginal_family(engine.sub_engine, i + (t-1)*ss);
wolffd@0: