comparison toolboxes/FullBNT-1.0.7/bnt/inference/dynamic/@stable_ho_inf_engine/enter_evidence.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 [engine, loglik] = enter_evidence(engine, evidence, varargin)
2 % ENTER_EVIDENCE Add the specified evidence to the network (jtree_unrolled_dbn)
3 % [engine, loglik] = enter_evidence(engine, evidence, ...)
4 %
5 % evidence{i,t} = [] if if X(i,t) is hidden, and otherwise contains its observed value (scalar or column vector)
6 %
7 % The following optional arguments can be specified in the form of name/value pairs:
8 % [default value in brackets]
9 %
10 % maximize - if 1, does max-product instead of sum-product [0]
11 % filter - if 1, does filtering (not supported), else smoothing [0]
12 %
13 % e.g., engine = enter_evidence(engine, ev, 'maximize', 1)
14
15 maximize = 0;
16 filter = 0;
17
18 % parse optional params
19 args = varargin;
20 nargs = length(args);
21 if nargs > 0
22 for i=1:2:nargs
23 switch args{i},
24 case 'maximize', maximize = args{i+1};
25 case 'filter', filter = args{i+1};
26 otherwise,
27 error(['invalid argument name ' args{i}]);
28 end
29 end
30 end
31
32 if filter
33 error('jtree_unrolled_dbn does not support filtering')
34 end
35
36 if size(evidence,2) ~= engine.nslices
37 error(['engine was created assuming there are ' num2str(engine.nslices) ...
38 ' slices, but evidence has ' num2str(size(evidence,2))])
39 end
40
41 [engine.unrolled_engine, loglik] = enter_evidence(engine.unrolled_engine, evidence, 'maximize', maximize);
42
43