Daniel@0: function bnet = mk_fhmm(N, Q, Y, discrete_obs)
Daniel@0: % MK_FHMM Make a factorial Hidden Markov Model
Daniel@0: %
Daniel@0: % There are N independent parallel hidden chains, each connected to the output
Daniel@0: %
Daniel@0: % e.g., N = 2 (vertical/diagonal edges point down)
Daniel@0: %
Daniel@0: % A1--->A2
Daniel@0: % | B1--|->B2
Daniel@0: % | /   |/
Daniel@0: % Y1    Y2
Daniel@0: %
Daniel@0: % [bnet, onode] = mk_chmm(n, q, y, discrete_obs)
Daniel@0: %
Daniel@0: % Each hidden node is discrete and has Q values.
Daniel@0: % If discrete_obs = 1, each observed node is discrete and has values 1..Y.
Daniel@0: % If discrete_obs = 0, each observed node is a Gaussian vector of length Y.
Daniel@0: 
Daniel@0: if nargin < 2, Q = 2; end
Daniel@0: if nargin < 3, Y = 2; end
Daniel@0: if nargin < 4, discrete_obs = 1; end
Daniel@0: 
Daniel@0: ss = N+1;
Daniel@0: hnodes = 1:N;
Daniel@0: onode = N+1;
Daniel@0: 
Daniel@0: intra = zeros(ss);
Daniel@0: intra(hnodes, onode) = 1;
Daniel@0: 
Daniel@0: inter = eye(ss);
Daniel@0: inter(onode,onode) = 0;
Daniel@0: 
Daniel@0: ns = [Q*ones(1,N) Y];
Daniel@0: 
Daniel@0: eclass1 = [hnodes onode];
Daniel@0: eclass2 = [hnodes+ss onode];
Daniel@0: if discrete_obs
Daniel@0:   dnodes = 1:ss;
Daniel@0: else
Daniel@0:   dnodes = hnodes;
Daniel@0: end
Daniel@0: bnet = mk_dbn(intra, inter, ns, 'discrete', dnodes, 'eclass1', eclass1, 'eclass2', eclass2, ...
Daniel@0: 	      'observed', onode);
Daniel@0: 
Daniel@0: for i=hnodes(:)'
Daniel@0:   bnet.CPD{i} = tabular_CPD(bnet, i);
Daniel@0: end
Daniel@0: i = onode;
Daniel@0: if discrete_obs
Daniel@0:   bnet.CPD{i} = tabular_CPD(bnet, i);
Daniel@0: else
Daniel@0:   bnet.CPD{i} = gaussian_CPD(bnet, i);
Daniel@0: end
Daniel@0: for i=hnodes(:)'+ss
Daniel@0:   bnet.CPD{i} = tabular_CPD(bnet, i);
Daniel@0: end
Daniel@0: 
Daniel@0: