Mercurial > hg > camir-aes2014
comparison toolboxes/FullBNT-1.0.7/bnt/examples/static/Belprop/belprop_polytree_cg.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 % Inference on a conditional Gaussian model | |
2 | |
3 % Make the following polytree, where all arcs point down | |
4 | |
5 % 1 2 | |
6 % \ / | |
7 % 3 | |
8 % / \ | |
9 % 4 5 | |
10 | |
11 N = 5; | |
12 dag = zeros(N,N); | |
13 dag(1,3) = 1; | |
14 dag(2,3) = 1; | |
15 dag(3, [4 5]) = 1; | |
16 | |
17 ns = [2 1 2 1 2]; | |
18 | |
19 dnodes = 1; | |
20 %onodes = [1 5]; | |
21 bnet = mk_bnet(dag, ns, 'discrete', dnodes, 'observed', dnodes); | |
22 | |
23 bnet.CPD{1} = tabular_CPD(bnet, 1); | |
24 for i=2:N | |
25 bnet.CPD{i} = gaussian_CPD(bnet, i); | |
26 end | |
27 | |
28 engine = {}; | |
29 engine{end+1} = jtree_inf_engine(bnet); | |
30 engine{end+1} = pearl_inf_engine(bnet, 'protocol', 'parallel'); | |
31 | |
32 [time, engine] = cmp_inference_static(bnet, engine, 'maximize', 0, 'check_ll', 0, ... | |
33 'singletons_only', 0, 'observed', [1 3]); |