blob: a651e1922861238349f18e615e58646866127014 (
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
26
27
|
function m = pot_to_marginal(pot)
% POT_TO_MARGINAL Convert a cgpot to a marginal structure.
% m = pot_to_marginal(pot)
pot = cg_can_to_mom(pot);
m.domain = pot.domain;
n = pot.csize;
d = length(pot.mom);
if n==0
m.mu = [];
m.Sigma = [];
else
m.mu = zeros(n, d);
m.Sigma = zeros(n, n, d);
end
m.T = 0*myones(pot.dsizes);
for i=1:pot.dsize
s = struct(pot.mom{i}); % violate privacy of object
if n > 0
m.mu(:,i) = s.mu;
m.Sigma(:,:,i) = s.Sigma;
end
m.T(i) = exp(s.logp);
end
if isvectorBNT(m.T)
m.T = m.T(:)';
end
|
