To check out this repository please hg clone the following URL, or open the URL using EasyMercurial or your preferred Mercurial client.

Statistics Download as Zip
| Branch: | Revision:

root / _FullBNT / BNT / general / shrink_obs_dims_in_gaussian.m @ 8:b5b38998ef3b

History | View | Annotate | Download (438 Bytes)

1
function marg2 = shrink_obs_dims_in_gaussian(marg1, dom, evidence, ns)
2
% SHRINK_OBS_DIMS_IN_GAUSSIAN Remove observed dimensions from mu/Sigma
3
% function marg2 = shrink_obs_dims_in_gaussian(marg1, dom, evidence, ns)
4

    
5
% This is used by loopy
6

    
7
hdom = dom(isemptycell(evidence(dom)));
8
ndx = find_equiv_posns(hdom, dom);
9
b = block(ndx, ns(dom));
10
marg2.mu = marg1.mu(b);
11
marg2.Sigma = marg1.Sigma(b,b);
12
marg2.domain = marg1.domain;