annotate toolboxes/FullBNT-1.0.7/bnt/inference/static/@gaussian_inf_engine/marginal_nodes.m @ 0:e9a9cd732c1e tip

first hg version after svn
author wolffd
date Tue, 10 Feb 2015 15:05:51 +0000
parents
children
rev   line source
wolffd@0 1 function marginal = marginal_nodes(engine, query)
wolffd@0 2 % MARGINAL_NODES Compute the marginal on the specified query nodes (gaussian)
wolffd@0 3 % marginal = marginal_nodes(engine, query)
wolffd@0 4
wolffd@0 5 % Compute sum_{Hsum} Pr(Hkeep, Hsum | o)
wolffd@0 6 H = engine.hnodes;
wolffd@0 7 bnet = bnet_from_engine(engine);
wolffd@0 8 ns = bnet.node_sizes;
wolffd@0 9 Hkeep = myintersect(H, query);
wolffd@0 10 Hsum = mysetdiff(H, Hkeep);
wolffd@0 11
wolffd@0 12 [marginal.mu, marginal.Sigma] = marginalize_gaussian(engine.Hmu, engine.HSigma, Hkeep, Hsum, ns);
wolffd@0 13 marginal.domain = query;
wolffd@0 14 marginal.T = 1;
wolffd@0 15