wolffd@0: function logp = enumerate_HMM_loglik(prior, transmat, obsmat)
wolffd@0: % ENUMERATE_HMM_LOGLIK Compute the log likelihood of a sequence by exhaustive (O(Q^T)) enumeration.
wolffd@0: % logp = enumerate_HMM_loglik(prior, transmat, obsmat)
wolffd@0: %
wolffd@0: % Inputs:
wolffd@0: % prior(i) = Pr(Q(1) = i)
wolffd@0: % transmat(i,j) = Pr(Q(t+1)=j | Q(t)=i)
wolffd@0: % obsmat(i,t) = Pr(y(t) | Q(t)=i)
wolffd@0: 
wolffd@0: Q = length(prior);
wolffd@0: T = size(obsmat, 2);
wolffd@0: sizes = repmat(Q, 1, T);
wolffd@0: 
wolffd@0: psum = 0;
wolffd@0: for i=1:Q^T
wolffd@0:   qs = ind2subv(sizes, i); % make the state sequence 
wolffd@0:   psum = psum + prob_path(prior, transmat, obsmat, qs);
wolffd@0: end
wolffd@0: logp = log(psum)
wolffd@0: 
wolffd@0: