annotate toolboxes/FullBNT-1.0.7/bnt/inference/dynamic/@ff_inf_engine/marginal_nodes.m @ 0:cc4b1211e677
tip
initial commit to HG from
Changeset:
646 (e263d8a21543) added further path and more save "camirversion.m"
| author |
Daniel Wolff |
| date |
Fri, 19 Aug 2016 13:07:06 +0200 |
| parents |
|
| children |
|
| rev |
line source |
|
Daniel@0
|
1 function marginal = marginal_nodes(engine, nodes, t)
|
|
Daniel@0
|
2 % MARGINAL_NODES Compute the marginal on the specified query nodes (ff)
|
|
Daniel@0
|
3 % marginal = marginal_nodes(engine, i, t)
|
|
Daniel@0
|
4 % returns Pr(X(i,t) | Y(1:T)), where X(i,t) is the i'th node in the t'th slice.
|
|
Daniel@0
|
5 % If enter_evidence used filtering instead of smoothing, this will return Pr(X(i,t) | Y(1:t)).
|
|
Daniel@0
|
6
|
|
Daniel@0
|
7 if nargin < 3, t = 1; end
|
|
Daniel@0
|
8 assert(length(nodes)==1);
|
|
Daniel@0
|
9 i = nodes(end);
|
|
Daniel@0
|
10 if myismember(i, engine.hnodes)
|
|
Daniel@0
|
11 marginal = pot_to_marginal(engine.marginals{i,t});
|
|
Daniel@0
|
12 else
|
|
Daniel@0
|
13 marginal = pot_to_marginal(dpot(i, 1, 1)); % observed
|
|
Daniel@0
|
14 end
|
|
Daniel@0
|
15
|
|
Daniel@0
|
16 bnet = bnet_from_engine(engine);
|
|
Daniel@0
|
17 ss = length(bnet.intra);
|
|
Daniel@0
|
18 % we convert the domain to the unrolled numbering system
|
|
Daniel@0
|
19 % so that update_ess extracts the right evidence.
|
|
Daniel@0
|
20 marginal.domain = nodes+(t-1)*ss;
|
