Mercurial > hg > camir-aes2014

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

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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 (pearl_dbn)
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.onodes)
wolffd@0 11 marginal.T = engine.marginal{i,t};
wolffd@0 12 else
wolffd@0 13 marginal.T = 1; % observed
wolffd@0 14 end
wolffd@0 15
wolffd@0 16 % we convert the domain to the unrolled numbering system
wolffd@0 17 % so that update_ess extracts the right evidence.
wolffd@0 18 marginal.domain = nodes+(t-1)*engine.ss;