annotate toolboxes/FullBNT-1.0.7/bnt/examples/static/SCG/scg1.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 % Same as cg1, except we call stab_cond_gauss_inf_engine
Daniel@0 2
Daniel@0 3 bnet = mk_incinerator_bnet;
Daniel@0 4
Daniel@0 5 engines = {};
Daniel@0 6 engines{end+1} = stab_cond_gauss_inf_engine(bnet);
Daniel@0 7 engines{end+1} = jtree_inf_engine(bnet);
Daniel@0 8 engines{end+1} = cond_gauss_inf_engine(bnet);
Daniel@0 9 nengines = length(engines);
Daniel@0 10
Daniel@0 11 F = 1; W = 2; E = 3; B = 4; C = 5; D = 6; Min = 7; Mout = 8; L = 9;
Daniel@0 12 n = 9;
Daniel@0 13 dnodes = [B F W];
Daniel@0 14 cnodes = mysetdiff(1:n, dnodes);
Daniel@0 15
Daniel@0 16 evidence = cell(1,n); % no evidence
Daniel@0 17 ll = zeros(1, nengines);
Daniel@0 18 for e=1:nengines
Daniel@0 19 [engines{e}, ll(e)] = enter_evidence(engines{e}, evidence);
Daniel@0 20 end
Daniel@0 21 %assert(approxeq(ll(1), ll)))
Daniel@0 22 ll
Daniel@0 23
Daniel@0 24 % Compare to the results in table on p1107.
Daniel@0 25 % These results are printed to 3dp in Cowell p150
Daniel@0 26
Daniel@0 27 mu = zeros(1,n);
Daniel@0 28 sigma = zeros(1,n);
Daniel@0 29 dprob = zeros(1,n);
Daniel@0 30 addev = 1;
Daniel@0 31 tol = 1e-2;
Daniel@0 32 for e=1:nengines
Daniel@0 33 for i=cnodes(:)'
Daniel@0 34 m = marginal_nodes(engines{e}, i, addev);
Daniel@0 35 mu(i) = m.mu;
Daniel@0 36 sigma(i) = sqrt(m.Sigma);
Daniel@0 37 end
Daniel@0 38 for i=dnodes(:)'
Daniel@0 39 m = marginal_nodes(engines{e}, i, addev);
Daniel@0 40 dprob(i) = m.T(1);
Daniel@0 41 end
Daniel@0 42 assert(approxeq(mu([E D C L Min Mout]), [-3.25 3.04 -1.85 1.48 -0.214 2.83], tol))
Daniel@0 43 assert(approxeq(sigma([E D C L Min Mout]), [0.709 0.770 0.507 0.631 0.459 0.860], tol))
Daniel@0 44 assert(approxeq(dprob([B F W]), [0.85 0.95 0.29], tol))
Daniel@0 45 %m = marginal_nodes(engines{e}, bnet.names('E'), addev);
Daniel@0 46 %assert(approxeq(m.mu, -3.25, tol))
Daniel@0 47 %assert(approxeq(sqrt(m.Sigma), 0.709, tol))
Daniel@0 48 end
Daniel@0 49
Daniel@0 50 % Add evidence (p 1105, top right)
Daniel@0 51 evidence = cell(1,n);
Daniel@0 52 evidence{W} = 1; % industrial
Daniel@0 53 evidence{L} = 1.1;
Daniel@0 54 evidence{C} = -0.9;
Daniel@0 55
Daniel@0 56 ll = zeros(1, nengines);
Daniel@0 57 for e=1:nengines
Daniel@0 58 [engines{e}, ll(e)] = enter_evidence(engines{e}, evidence);
Daniel@0 59 end
Daniel@0 60 %assert(all(approxeq(ll(1), ll)))
Daniel@0 61 ll
Daniel@0 62
Daniel@0 63 for e=1:nengines
Daniel@0 64 for i=cnodes(:)'
Daniel@0 65 m = marginal_nodes(engines{e}, i, addev);
Daniel@0 66 mu(i) = m.mu;
Daniel@0 67 sigma(i) = sqrt(m.Sigma);
Daniel@0 68 end
Daniel@0 69 for i=dnodes(:)'
Daniel@0 70 m = marginal_nodes(engines{e}, i, addev);
Daniel@0 71 dprob(i) = m.T(1);
Daniel@0 72 end
Daniel@0 73 assert(approxeq(mu([E D C L Min Mout]), [-3.90 3.61 -0.9 1.1 0.5 4.11], tol))
Daniel@0 74 assert(approxeq(sigma([E D C L Min Mout]), [0.076 0.326 0 0 0.1 0.344], tol))
Daniel@0 75 assert(approxeq(dprob([B F W]), [0.0122 0.9995 1], tol))
Daniel@0 76 end
Daniel@0 77