wolffd@0: % Make the following HHMM
wolffd@0: %
wolffd@0: %     LH                  RH
wolffd@0: %    /                      \
wolffd@0: %   /                        \
wolffd@0: %  LR -> UD -> RL -> DU       RL -> UD -> LR -> DU
wolffd@0: %   \
wolffd@0: %    \
wolffd@0: %     Q1 -> Q2
wolffd@0: %
wolffd@0: % where level 1 is fully interconnected (not shown)
wolffd@0: % level 2 is left-right
wolffd@0: % and each model at level 3 is a 2 state LR shared HMM 
wolffd@0: 
wolffd@0: Qsizes = [2 4 2];
wolffd@0: D = 3;
wolffd@0: 
wolffd@0: % LEVEL 1
wolffd@0: 
wolffd@0: startprob1 = 'ergodic';
wolffd@0: transprob1 = 'ergodic';
wolffd@0: 
wolffd@0: 
wolffd@0: % LEVEL 2
wolffd@0: 
wolffd@0: startprob = zeros(2, 4);
wolffd@0: %        Q1  Q2
wolffd@0: startprob(1, 1) = 1;
wolffd@0: startprob(2, 3) = 1;
wolffd@0: 
wolffd@0: transprob = zeros(2, 4, 4);
wolffd@0: transprob(1,:,:) = [0 1 0 0
wolffd@0: 		    0 0 1 0
wolffd@0: 		    0 0 0 1
wolffd@0: 		    0 0 0 1];
wolffd@0: transprob(2,:,:) = [0 0 0 1
wolffd@0: 		    1 0 0 0
wolffd@0: 		    0 1 0 0
wolffd@0: 		    0 0 0 1];
wolffd@0: 
wolffd@0: Q2args = {'startprob', startprob, 'transprob', transprob};
wolffd@0: 
wolffd@0: % always terminate in state 4 (default)
wolffd@0: % F2args
wolffd@0: 
wolffd@0: % LEVEL 3
wolffd@0: 
wolffd@0: % Defaults are fine: always start in state 1, left-right model, finish in state 2
wolffd@0: 
wolffd@0: 
wolffd@0: % OBS LEVEl
wolffd@0: 
wolffd@0: chars = ['L', 'l', 'U', 'u', 'R', 'r', 'D', 'd'];
wolffd@0: Osize = length(chars);
wolffd@0: 
wolffd@0: obsprob = zeros([4 2 Osize]);
wolffd@0: %       Q2 Q3 O
wolffd@0: obsprob(1, 1, find(chars == 'L')) =  1.0;
wolffd@0: obsprob(1, 2, find(chars == 'l')) =  1.0;
wolffd@0: 
wolffd@0: obsprob(2, 1, find(chars == 'U')) =  1.0;
wolffd@0: obsprob(2, 2, find(chars == 'u')) =  1.0;
wolffd@0: 
wolffd@0: obsprob(3, 1, find(chars == 'R')) =  1.0;
wolffd@0: obsprob(3, 2, find(chars == 'r')) =  1.0;
wolffd@0: 
wolffd@0: obsprob(4, 1, find(chars == 'D')) =  1.0;
wolffd@0: obsprob(4, 2, find(chars == 'd')) =  1.0;
wolffd@0: 
wolffd@0: Oargs = {'CPT', obsprob};
wolffd@0: 
wolffd@0: 
wolffd@0: bnet = mk_hhmm3('Qsizes', Qsizes, 'Osize', Osize', 'discrete_obs', 1, 'Oargs', Oargs, 'Q1args', Q1args, 'Q2args', Q2args);
wolffd@0: 
wolffd@0: T = 20;
wolffd@0: usecell = 0;
wolffd@0: evidence = sample_dbn(bnet, T, usecell);      
wolffd@0: %chars(evidence(end,:))
wolffd@0: 
wolffd@0: Q1 = 1; Q2 = 2; Q3 = 3; F3 = 4; F2 = 5; obs = 6;
wolffd@0: Qnodes = [Q1 Q2 Q3]; Fnodes = [F2 F3];
wolffd@0: 
wolffd@0: pretty_print_hhmm_parse(evidence, Qnodes, Fnodes, obs, chars);
wolffd@0: 
wolffd@0: eclass = bnet.equiv_class;
wolffd@0: S=struct(bnet.CPD{eclass(Q2,2)})