Mercurial > hg > camir-aes2014
view toolboxes/FullBNT-1.0.7/bnt/CPDs/@discrete_CPD/Old/prob_CPD.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 p = prob_CPD(CPD, domain, ns, cnodes, evidence) % PROB_CPD Compute prob of a node given evidence on the parents (discrete) % p = prob_CPD(CPD, domain, ns, cnodes, evidence) % % domain is the domain of CPD. % node_sizes(i) is the size of node i. % cnodes = all the cts nodes % evidence{i} is the evidence on the i'th node. ps = domain(1:end-1); self = domain(end); CPT = CPD_to_CPT(CPD); if isempty(ps) T = CPT; else assert(~any(isemptycell(evidence(ps)))); pvals = cat(1, evidence{ps}); i = subv2ind(ns(ps), pvals(:)'); T = reshape(CPT, [prod(ns(ps)) ns(self)]); T = T(i,:); end p = T(evidence{self});