Mauch MIREX 2010

function p = prob_node(CPD, self_ev, pev)

% PROB_NODE Compute P(y|pa(y), theta) (tabular)

% p = prob_node(CPD, self_ev, pev)

%

% self_ev{m} is the evidence on this node in case m

% pev{i,m} is the evidence on the i'th parent in case m

% If there is a single case, self_ev can be a scalar instead of a cell array

ncases = size(pev, 2);

%assert(~any(isemptycell(pev))); % slow

%assert(~any(isemptycell(self_ev))); % slow

CPT = CPD_to_CPT(CPD);  

sz = mysize(CPT);

nparents = length(sz)-1;

assert(nparents == size(pev, 1));

if ncases==1 

  x = cat(1, pev{:});

  if iscell(y)

    y = self_ev{1};

  else

    y = self_ev;

  end

  switch nparents

   case 0, p = CPT(y);

   case 1, p = CPT(x(1), y);

   case 2, p = CPT(x(1), x(2), y);

   case 3, p = CPT(x(1), x(2), x(3), y);

   otherwise,

    ind = subv2ind(CPD.sizes, [x y]);

    p = CPT(ind);

  end

else

  x = num2cell(pev)'; % each row is a case

  y = cat(1, self_ev{:})';

  ind = subv2ind(CPD.sizes, [x y]);

  p = CPT(ind);

end