blob: c0a79bdaf29e3ed76de4374efc96a00c6a06177c (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
function p = prob_CPD(CPD, domain, ns, cnodes, evidence)
% PROB_CPD Compute prob of a node given evidence on the parents (discrete)
% p = prob_CPD(CPD, domain, ns, cnodes, evidence)
%
% domain is the domain of CPD.
% node_sizes(i) is the size of node i.
% cnodes = all the cts nodes
% evidence{i} is the evidence on the i'th node.
ps = domain(1:end-1);
self = domain(end);
CPT = CPD_to_CPT(CPD);
if isempty(ps)
T = CPT;
else
assert(~any(isemptycell(evidence(ps))));
pvals = cat(1, evidence{ps});
i = subv2ind(ns(ps), pvals(:)');
T = reshape(CPT, [prod(ns(ps)) ns(self)]);
T = T(i,:);
end
p = T(evidence{self});
|
