wolffd@0: % Make the HHMM in Figure 1 of the NIPS'01 paper
wolffd@0: 
wolffd@0: Qsize = [2 3 2];
wolffd@0: D = 3;
wolffd@0: 
wolffd@0: % transprob{d}(i,k,j), transprob{1}(i,j)
wolffd@0: % termprob{d}(k,j), termprob{1}(1,j)
wolffd@0: % startprob{d}(k,j), startprob{1}(1,j)
wolffd@0: % obsprob(k, o) for discrete outputs
wolffd@0:     
wolffd@0: % LEVEL 1
wolffd@0: %       1 2 e
wolffd@0: A{1} = [0 0 1;
wolffd@0: 	0 0 1];
wolffd@0: [transprob{1}, termprob{1}] = remove_hhmm_end_state(A{1});
wolffd@0: startprob{1} = [0.5 0.5];
wolffd@0: Q1args = {'startprob', startprob{1}, 'transprob', transprob{1}};
wolffd@0: 
wolffd@0: % LEVEL 2
wolffd@0: A{2} = zeros(Qsize(2), Qsize(1), Qsize(2)+1);
wolffd@0: 
wolffd@0: %              1 2 3 e
wolffd@0: A{2}(:,1,:) = [0 1 0 0
wolffd@0: 	       0 0 1 0
wolffd@0: 	       0 0 0 1];
wolffd@0: 
wolffd@0: %              1 2 3 e
wolffd@0: A{2}(:,2,:) = [0 1 0 0
wolffd@0: 	       0 0 1 0
wolffd@0: 	       0 0 0 1];
wolffd@0: 
wolffd@0: [transprob{2}, termprob{2}] = remove_hhmm_end_state(A{2});	       
wolffd@0: 
wolffd@0: % always enter level 2 in state 1
wolffd@0: startprob{2} = [1 0 0
wolffd@0: 		1 0 0];
wolffd@0: 
wolffd@0: Q2args = {'startprob', startprob{2}, 'transprob', transprob{2}};
wolffd@0: F2args = {'CPT', termprob{2}};
wolffd@0: 
wolffd@0: 
wolffd@0: % LEVEL 3
wolffd@0: 
wolffd@0: A{3} = zeros([Qsize(3) Qsize(1:2) Qsize(3)+1]);
wolffd@0: endstate = Qsize(3)+1;
wolffd@0: %    Qt-1(3) Qt(1) Qt(2) Qt(3)
wolffd@0: %                               1   2   e
wolffd@0: A{3}(1,      1,    1,    endstate) = 1.0;
wolffd@0: A{3}(:,      1,    2,    :) = [0.0 1.0 0.0
wolffd@0:             	               0.5 0.0 0.5];
wolffd@0: A{3}(1,      1,    3,    endstate) = 1.0;
wolffd@0: 
wolffd@0: A{3}(1,      2,    1,    endstate) = 1.0;
wolffd@0: A{3}(:,      2,    2,    :) = [0.0 1.0 0.0
wolffd@0:             	               0.5 0.0 0.5];
wolffd@0: A{3}(1,      2,    3,    endstate) = 1.0;
wolffd@0: 
wolffd@0: A{3} = reshape(A{3}, [Qsize(3) prod(Qsize(1:2)) Qsize(3)+1]);
wolffd@0: [transprob{3}, termprob{3}] = remove_hhmm_end_state(A{3});	       
wolffd@0: 
wolffd@0: % define the vertical entry points to level 3
wolffd@0: startprob{3} = zeros(Qsize);
wolffd@0: %            Q1 Q2 Q3
wolffd@0: startprob{3}(1, 1, 1) = 1.0;
wolffd@0: startprob{3}(1, 2, 1) = 1.0;
wolffd@0: startprob{3}(1, 3, 1) = 1.0;
wolffd@0: 
wolffd@0: startprob{3}(2, 1, 1) = 1.0;
wolffd@0: startprob{3}(2, 2, 1) = 1.0;
wolffd@0: startprob{3}(2, 3, 1) = 1.0;
wolffd@0: 
wolffd@0: startprob{3} = reshape(startprob{3}, prod(Qsize(1:2)), Qsize(3));
wolffd@0: 
wolffd@0: chars = ['a', 'b', 'c', 'd', 'x', 'y'];
wolffd@0: Osize = length(chars);
wolffd@0: 
wolffd@0: obsprob = zeros([Qsize Osize]);
wolffd@0: %       1 2 3 O
wolffd@0: obsprob(1,1,1,find(chars == 'a')) =  1.0;
wolffd@0: 
wolffd@0: obsprob(1,2,1,find(chars == 'x')) =  1.0;
wolffd@0: obsprob(1,2,2,find(chars == 'y')) =  1.0;
wolffd@0: 
wolffd@0: obsprob(1,3,1,find(chars == 'b')) =  1.0;
wolffd@0: 
wolffd@0: obsprob(2,1,1,find(chars == 'c')) =  1.0;
wolffd@0: 
wolffd@0: obsprob(2,2,1,find(chars == 'x')) =  1.0;
wolffd@0: obsprob(2,2,2,find(chars == 'y')) =  1.0;
wolffd@0: 
wolffd@0: obsprob(2,3,1,find(chars == 'd')) =  1.0;
wolffd@0: 
wolffd@0: obsprob = reshape(obsprob, prod(Qsize), Osize);
wolffd@0: 
wolffd@0: [intra, inter, Qnodes, Fnodes, Onode] = mk_hhmm_topo(D);
wolffd@0: 
wolffd@0: hhmm.Qnodes = Qnodes;
wolffd@0: hhmm.Fnodes = Fnodes;
wolffd@0: hhmm.Onode = Onode;
wolffd@0: hhmm.D = D;
wolffd@0: hhmm.Qsize = Qsize;
wolffd@0: hhmm.Osize = Osize;
wolffd@0: hhmm.startprob = startprob;
wolffd@0: hhmm.transprob = transprob;
wolffd@0: hhmm.termprob = termprob;
wolffd@0: hhmm.obsprob = obsprob;
wolffd@0: hhmm.A = A;
wolffd@0: 
wolffd@0: