view toolboxes/FullBNT-1.0.7/bnt/general/mk_slice_and_half_dbn.m @ 0:e9a9cd732c1e tip

first hg version after svn

function bnet15 = mk_slice_and_half_dbn(bnet, int)
% function bnet = mk_slice_and_half_dbn(bnet, int)
% function bnet = mk_slice_and_half_dbn(bnet, int)
%
% Create a "1.5 slice" jtree, containing the interface nodes of slice 1
% and all the nodes of slice 2
% To keep the node numbering the same, we simply disconnect the non-interface nodes
% from slice 1, and set their size to 1.
% We do this to speed things up, and so that the likelihood is computed correctly.
% We do not need to do
% this if we just want to compute marginals (i.e., we can include nodes whose potentials will
% be left as all 1s).

intra15 = bnet.intra;
ss = length(bnet.intra);
nonint = mysetdiff(1:ss, int);
for i=nonint(:)'
  intra15(:,i) = 0;
  intra15(i,:) = 0;
  %assert(~any(bnet.inter(i,:)))
end
dag15 = [intra15      bnet.inter;
	 zeros(ss)    bnet.intra];
ns = bnet.node_sizes(:);
ns(nonint) = 1; % disconnected nodes get size 1
obs_nodes = [bnet.observed(:) bnet.observed(:)+ss];
bnet15 = mk_bnet(dag15, ns, 'discrete', bnet.dnodes, 'equiv_class', bnet.equiv_class(:), ...
		 'observed', obs_nodes(:));