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