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

[Q O] = size(obsmat);
T = prod(size(data)); % length(data);
B = zeros(Q,T);

for t=1:T
  B(:,t) = obsmat(:, data(t));
end