annotate toolboxes/FullBNT-1.0.7/bnt/CPDs/@discrete_CPD/Old/convert_to_table.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 T = convert_to_table(CPD, domain, local_ev, obs_bitv)
Daniel@0 2 % CONVERT_TO_TABLE Convert a discrete CPD to a table
Daniel@0 3 % function T = convert_to_table(CPD, domain, local_ev, obs_bitv)
Daniel@0 4 %
Daniel@0 5 % We convert the CPD to a CPT, and then lookup the evidence on the discrete parents.
Daniel@0 6 % The resulting table can easily be converted to a potential.
Daniel@0 7
Daniel@0 8
Daniel@0 9 CPT = CPD_to_CPT(CPD);
Daniel@0 10 obs_child_only = ~any(obs_bitv(1:end-1)) & obs_bitv(end);
Daniel@0 11
Daniel@0 12 if obs_child_only
Daniel@0 13 sz = size(CPT);
Daniel@0 14 CPT = reshape(CPT, prod(sz(1:end-1)), sz(end));
Daniel@0 15 o = local_ev{end};
Daniel@0 16 T = CPT(:, o);
Daniel@0 17 else
Daniel@0 18 odom = domain(obs_bitv);
Daniel@0 19 vals = cat(1, local_ev{find(obs_bitv)}); % undo cell array
Daniel@0 20 map = find_equiv_posns(odom, domain);
Daniel@0 21 index = mk_multi_index(length(domain), map, vals);
Daniel@0 22 T = CPT(index{:});
Daniel@0 23 end