annotate toolboxes/FullBNT-1.0.7/bnt/general/compute_joint_pot.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 [jpot, loglik] = compute_joint_pot(bnet, nodes, evidence, domain)
|
|
Daniel@0
|
2 % COMPUTE_JOINT_POT Compute the global joint potential of a Bayes net
|
|
Daniel@0
|
3 % function jpot = compute_joint_pot(bnet, nodes, evidence, domain)
|
|
Daniel@0
|
4
|
|
Daniel@0
|
5 if nargin < 4, domain = nodes; end
|
|
Daniel@0
|
6
|
|
Daniel@0
|
7 onodes = find(~isemptycell(evidence));
|
|
Daniel@0
|
8 pot_type = determine_pot_type(bnet, onodes, domain);
|
|
Daniel@0
|
9
|
|
Daniel@0
|
10 jpot = mk_initial_pot(pot_type, domain, bnet.node_sizes, bnet.cnodes, onodes);
|
|
Daniel@0
|
11 for i=nodes(:)'
|
|
Daniel@0
|
12 e = bnet.equiv_class(i);
|
|
Daniel@0
|
13 fam = family(bnet.dag, i);
|
|
Daniel@0
|
14 pot = convert_to_pot(bnet.CPD{e}, pot_type, fam(:), evidence);
|
|
Daniel@0
|
15 jpot = multiply_by_pot(jpot, pot);
|
|
Daniel@0
|
16 end
|
|
Daniel@0
|
17 %[jpot, loglik] = normalize_pot(jpot); % causes errors in asia_dt1 etc
|
