wolffd@0:  Example due to Wang Hee Lin" <engp1622@nus.edu.sg
wolffd@0: 
wolffd@0: 
wolffd@0: intra = zeros(2);
wolffd@0: intra(1,2) = 1; 
wolffd@0: inter = zeros(2);
wolffd@0: inter(1,1) = 1; 
wolffd@0: 
wolffd@0: Q = 2; % num hidden states
wolffd@0: O = 2; % num observable symbols
wolffd@0: ns = [Q O];%number of states
wolffd@0: dnodes = 1:2;
wolffd@0: %onodes = [1:2]; % only possible with jtree, not hmm
wolffd@0: onodes = [2]; 
wolffd@0: bnet = mk_dbn(intra, inter, ns, 'discrete', dnodes, 'observed', onodes);
wolffd@0: for i=1:4
wolffd@0:   bnet.CPD{i} = tabular_CPD(bnet, i);
wolffd@0: end
wolffd@0: 
wolffd@0: prior0 = normalise(rand(Q,1));
wolffd@0: transmat0 = mk_stochastic(rand(Q,Q));
wolffd@0: obsmat0 = mk_stochastic(rand(Q,O));
wolffd@0: 
wolffd@0: %engine = smoother_engine(hmm_2TBN_inf_engine(bnet));
wolffd@0: engine = smoother_engine(jtree_2TBN_inf_engine(bnet));
wolffd@0: 
wolffd@0: ss = 2;%slice size(ss)
wolffd@0: ncases = 10;%number of examples
wolffd@0: T=10;
wolffd@0: max_iter=2;%iterations for EM
wolffd@0: cases = cell(1, ncases);
wolffd@0: for i=1:ncases
wolffd@0:   ev = sample_dbn(bnet, T);
wolffd@0:   cases{i} = cell(ss,T);
wolffd@0:   cases{i}(onodes,:) = ev(onodes, :);
wolffd@0: end
wolffd@0: [bnet2, LLtrace] = learn_params_dbn_em(engine, cases, 'max_iter', 4);