wolffd@0: function m = marginal_nodes(engine, query)
wolffd@0: % MARGINAL_NODES Compute the marginal on the specified query nodes (quickscore)
wolffd@0: % marginal = marginal_nodes(engine, query)
wolffd@0: %
wolffd@0: % 'query' must be a single disease (root) node.
wolffd@0: 
wolffd@0: assert(length(query)==1);
wolffd@0: p = engine.post(query);
wolffd@0: m.T = [1-p p]';
wolffd@0: m.domain = query;
wolffd@0: