Mercurial > hg > camir-aes2014
view toolboxes/FullBNT-1.0.7/bnt/CPDs/@gaussian_CPD/convert_to_table.m @ 0:e9a9cd732c1e tip
first hg version after svn
author | wolffd |
---|---|
date | Tue, 10 Feb 2015 15:05:51 +0000 |
parents | |
children |
line wrap: on
line source
function T = convert_to_table(CPD, domain, evidence) % CONVERT_TO_TABLE Convert a Gaussian CPD to a table % T = convert_to_table(CPD, domain, evidence) sz = CPD.sizes; ns = zeros(1, max(domain)); ns(domain) = sz; odom = domain(~isemptycell(evidence(domain))); ps = domain(1:end-1); cps = ps(CPD.cps); dps = ps(CPD.dps); self = domain(end); cdom = [cps(:)' self]; ddom = dps; cnodes = cdom; [m, C, W] = gaussian_CPD_params_given_dps(CPD, domain, evidence); ns(odom) = 1; dpsize = prod(ns(dps)); self = domain(end); assert(myismember(self, odom)); self_val = evidence{self}; T = zeros(dpsize, 1); if length(cps) > 0 assert(~any(isemptycell(evidence(cps)))); cps_vals = cat(1, evidence{cps}); for i=1:dpsize T(i) = gaussian_prob(self_val, m(:,i) + W(:,:,i)*cps_vals, C(:,:,i)); end else for i=1:dpsize T(i) = gaussian_prob(self_val, m(:,i), C(:,:,i)); end end