wolffd@0: function B = eval_pdf_cond_multinomial(data, obsmat)
wolffd@0: % EVAL_PDF_COND_MULTINOMIAL Evaluate pdf of conditional multinomial 
wolffd@0: % function B = eval_pdf_cond_multinomial(data, obsmat)
wolffd@0: %
wolffd@0: % Notation: Y = observation (O values), Q = conditioning variable (K values)
wolffd@0: %
wolffd@0: % Inputs:
wolffd@0: % data(t) = t'th observation - must be an integer in {1,2,...,K}: cannot be 0!
wolffd@0: % obsmat(i,o) = Pr(Y(t)=o | Q(t)=i)
wolffd@0: %
wolffd@0: % Output:
wolffd@0: % B(i,t) = Pr(y(t) | Q(t)=i)
wolffd@0: 
wolffd@0: [Q O] = size(obsmat);
wolffd@0: T = prod(size(data)); % length(data);
wolffd@0: B = zeros(Q,T);
wolffd@0: 
wolffd@0: for t=1:T
wolffd@0:   B(:,t) = obsmat(:, data(t));
wolffd@0: end