Daniel@0: function bnet = mutilate_bnet(bnet, nodes, vals)
Daniel@0: % MUTILATE_BNET Clamp nodes to specific values (perform a surgical intervention)
Daniel@0: % bnet = mutilate_bnet(bnet, nodes, vals)
Daniel@0: %
Daniel@0: % We make all the clamped nodes roots.
Daniel@0: 
Daniel@0: ns = bnet.node_sizes;
Daniel@0: for i=1:length(nodes)
Daniel@0:   X = nodes(i);
Daniel@0:   x = vals(i);
Daniel@0:   bnet.dag(:,X) = 0;
Daniel@0:   bnet.CPD{X} = root_CPD(bnet, X, x);
Daniel@0: end