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
26
27
28
29
30
31
32
33
34
35
36
37
38
|
function T = convert_to_table(CPD, domain, evidence)
% CONVERT_TO_TABLE Convert a Gaussian CPD to a table
% T = convert_to_table(CPD, domain, evidence)
sz = CPD.sizes;
ns = zeros(1, max(domain));
ns(domain) = sz;
odom = domain(~isemptycell(evidence(domain)));
ps = domain(1:end-1);
cps = ps(CPD.cps);
dps = ps(CPD.dps);
self = domain(end);
cdom = [cps(:)' self];
ddom = dps;
cnodes = cdom;
[m, C, W] = gaussian_CPD_params_given_dps(CPD, domain, evidence);
ns(odom) = 1;
dpsize = prod(ns(dps));
self = domain(end);
assert(myismember(self, odom));
self_val = evidence{self};
T = zeros(dpsize, 1);
if length(cps) > 0
assert(~any(isemptycell(evidence(cps))));
cps_vals = cat(1, evidence{cps});
for i=1:dpsize
T(i) = gaussian_prob(self_val, m(:,i) + W(:,:,i)*cps_vals, C(:,:,i));
end
else
for i=1:dpsize
T(i) = gaussian_prob(self_val, m(:,i), C(:,:,i));
end
end
|
