annotate toolboxes/FullBNT-1.0.7/bnt/general/convert_dbn_CPDs_to_pots.m @ 0:cc4b1211e677 tip

initial commit to HG from Changeset: 646 (e263d8a21543) added further path and more save "camirversion.m"
author Daniel Wolff
date Fri, 19 Aug 2016 13:07:06 +0200
parents
children
rev   line source
Daniel@0 1 function CPDpot = convert_dbn_CPDs_to_pots(bnet, evidence, pot_type, softCPDpot)
Daniel@0 2 % CONVERT_DBN_CPDS_TO_POTS Convert CPDs of (possibly instantiated) DBN nodes to potentials
Daniel@0 3 % CPDpot = convert_dbn_CPDs_to_pots(bnet, evidence, pot_type, softCPDpot)
Daniel@0 4 %
Daniel@0 5 % CPDpot{n,t} is a potential containing P(n,t|pa(n,t), ev)
Daniel@0 6 % softCPDpot{n,t} is a potential containing P(n,t|pa(n,t), ev) insted of using n's CPD
Daniel@0 7
Daniel@0 8 [ss T] = size(evidence);
Daniel@0 9
Daniel@0 10 if nargin < 4, softCPDpot = cell(ss,T); end
Daniel@0 11 CPDpot = softCPDpot;
Daniel@0 12
Daniel@0 13 % Convert CPDs of instantiated nodes to potential form
Daniel@0 14 t = 1;
Daniel@0 15 for n=1:ss
Daniel@0 16 fam = family(bnet.dag, n);
Daniel@0 17 e = bnet.equiv_class(n, 1);
Daniel@0 18 if isempty(softCPDpot{n,t})
Daniel@0 19 CPDpot{n,t} = convert_to_pot(bnet.CPD{e}, pot_type, fam(:), evidence(:,1));
Daniel@0 20 end
Daniel@0 21 end
Daniel@0 22 for n=1:ss
Daniel@0 23 fam = family(bnet.dag, n, 2);
Daniel@0 24 e = bnet.equiv_class(n, 2);
Daniel@0 25 for t=2:T
Daniel@0 26 if isempty(softCPDpot{n,t})
Daniel@0 27 CPDpot{n,t} = convert_to_pot(bnet.CPD{e}, pot_type, fam(:), evidence(:,t-1:t));
Daniel@0 28 end
Daniel@0 29 end
Daniel@0 30 end