blob: 460ec782ea6da7de86b834ab7d13a2aa10472ece (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
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
|
