Mercurial > hg > camir-aes2014
comparison toolboxes/FullBNT-1.0.7/bnt/inference/static/@gibbs_sampling_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) | |
2 % ENTER_EVIDENCE Add the specified evidence to the network (gibbs_sampling_inf_engine) | |
3 % [engine, loglik] = enter_evidence(engine, evidence) | |
4 % | |
5 % evidence{i} = [] if if X(i) is hidden, and otherwise contains its observed value | |
6 % | |
7 % loglik is not computed... we just return a 0 value | |
8 | |
9 bnet = bnet_from_engine(engine); | |
10 | |
11 engine.hnodes = find(isemptycell(evidence)); | |
12 engine.onodes = mysetdiff(1:length(evidence), engine.hnodes); | |
13 | |
14 engine.evidence = zeros(engine.slice_size, 1); | |
15 | |
16 % Reset all counts since they are no longer valid | |
17 engine.marginal_counts = {}; | |
18 %engine.state = sample_bnet (bnet, 1, 0); | |
19 engine.state = cell2num(sample_bnet(bnet)); | |
20 | |
21 % For speed, we use a normal (not cell) array. We're making use of | |
22 % the current restriction to discrete nodes. | |
23 for i = engine.onodes | |
24 engine.evidence(i) = evidence{i}; | |
25 end | |
26 | |
27 loglik = 0; | |
28 | |
29 |