Mercurial > hg > camir-aes2014
annotate toolboxes/FullBNT-1.0.7/bnt/inference/dynamic/@ff_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, nodes, t) |
wolffd@0 | 2 % MARGINAL_NODES Compute the marginal on the specified query nodes (ff) |
wolffd@0 | 3 % marginal = marginal_nodes(engine, i, t) |
wolffd@0 | 4 % returns Pr(X(i,t) | Y(1:T)), where X(i,t) is the i'th node in the t'th slice. |
wolffd@0 | 5 % If enter_evidence used filtering instead of smoothing, this will return Pr(X(i,t) | Y(1:t)). |
wolffd@0 | 6 |
wolffd@0 | 7 if nargin < 3, t = 1; end |
wolffd@0 | 8 assert(length(nodes)==1); |
wolffd@0 | 9 i = nodes(end); |
wolffd@0 | 10 if myismember(i, engine.hnodes) |
wolffd@0 | 11 marginal = pot_to_marginal(engine.marginals{i,t}); |
wolffd@0 | 12 else |
wolffd@0 | 13 marginal = pot_to_marginal(dpot(i, 1, 1)); % observed |
wolffd@0 | 14 end |
wolffd@0 | 15 |
wolffd@0 | 16 bnet = bnet_from_engine(engine); |
wolffd@0 | 17 ss = length(bnet.intra); |
wolffd@0 | 18 % we convert the domain to the unrolled numbering system |
wolffd@0 | 19 % so that update_ess extracts the right evidence. |
wolffd@0 | 20 marginal.domain = nodes+(t-1)*ss; |