wolffd@0: function [ll, p] = prob_path(prior, transmat, obsmat, qs)
wolffd@0: % PROB_PATH Compute the prob. of a specific path (state sequence) through an HMM.
wolffd@0: % [ll, p] = prob_path(prior, transmat, obsmat, states)
wolffd@0: %
wolffd@0: % ll = log prob path
wolffd@0: % p(t) = Pr(O(t)) * Pr(Q(t) -> Q(t+1)) for 1<=t<T, p(T) = Pr(O(T))
wolffd@0: 
wolffd@0: T = size(obsmat, 2);
wolffd@0: p = zeros(1,T);
wolffd@0: p(1) = prior(qs(1)) * obsmat(qs(1),1);
wolffd@0: for t=2:T
wolffd@0:   p(t) = transmat(qs(t-1), qs(t)) * obsmat(qs(t),t);
wolffd@0: end
wolffd@0: 
wolffd@0: ll = sum(log(p));