annotate toolboxes/FullBNT-1.0.7/bnt/inference/dynamic/@hmm_inf_engine/find_mpe.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 mpe = find_mpe(engine, evidence)
|
|
Daniel@0
|
2 % FIND_MPE Find the most probable explanation (Viterbi)
|
|
Daniel@0
|
3 % mpe = enter_evidence(engine, evidence, ...)
|
|
Daniel@0
|
4 %
|
|
Daniel@0
|
5 % evidence{i,t} = [] if if X(i,t) is hidden, and otherwise contains its observed value (scalar or column vector)
|
|
Daniel@0
|
6 %
|
|
Daniel@0
|
7
|
|
Daniel@0
|
8 obslik = mk_hmm_obs_lik_matrix(engine, evidence);
|
|
Daniel@0
|
9 path = viterbi_path(engine.startprob, engine.transprob, obslik);
|
|
Daniel@0
|
10 bnet = bnet_from_engine(engine);
|
|
Daniel@0
|
11 ns = bnet.node_sizes_slice;
|
|
Daniel@0
|
12 ns(bnet.observed) = 1;
|
|
Daniel@0
|
13 ass = ind2subv(ns, path);
|
|
Daniel@0
|
14 mpe = num2cell(ass');
|
|
Daniel@0
|
15 mpe(bnet.observed,:) = evidence(bnet.observed,:);
|
|
Daniel@0
|
16
|
