Mercurial > hg > camir-aes2014
comparison toolboxes/FullBNT-1.0.7/bnt/inference/dynamic/@pearl_dbn_inf_engine/private/init_pearl_dbn_ev_msgs.m @ 0:e9a9cd732c1e tip
first hg version after svn
author | wolffd |
---|---|
date | Tue, 10 Feb 2015 15:05:51 +0000 |
parents | |
children |
comparison
equal
deleted
inserted
replaced
-1:000000000000 | 0:e9a9cd732c1e |
---|---|
1 function msg = init_pearl_dbn_ev_msgs(bnet, evidence, engine) | |
2 | |
3 [ss T] = size(evidence); | |
4 pot_type = 'd'; | |
5 | |
6 % each hidden node absorbs lambda from its observed child (if any) | |
7 for t=1:T | |
8 for i=hnodes | |
9 c = engine.obschild(i); | |
10 if c > 0 | |
11 if t==1 | |
12 fam = family(bnet.dag, c); | |
13 e = bnet.equiv_class(c, 1); | |
14 CPDpot = CPD_to_pot(pot_type, bnet.CPD{e}, fam, bnet.node_sizes(:), bnet.cnodes(:), evidence(:,1)); | |
15 else | |
16 fam = family(bnet.dag, c, 2); % within 2 slice network | |
17 e = bnet.equiv_class(c, 2); | |
18 CPDpot = CPD_to_pot(pot_type, bnet.CPD{e}, fam, bnet.node_sizes(:), bnet.cnodes(:), evidence(:,t-1:t)); | |
19 end | |
20 temp = pot_to_marginal(CPDpot); | |
21 n = i + (t-1)*ss; | |
22 lam_msg = normalise(temp.T); | |
23 j = engine.child_index{n}(c+(t-1)*ss); | |
24 assert(j==1); | |
25 msg{n}.lambda_from_child{j} = lam_msg; | |
26 end | |
27 end | |
28 end |