Mercurial > hg > camir-aes2014
annotate toolboxes/FullBNT-1.0.7/bnt/CPDs/@discrete_CPD/CPD_to_lambda_msg.m @ 0:e9a9cd732c1e tip
first hg version after svn
author | wolffd |
---|---|
date | Tue, 10 Feb 2015 15:05:51 +0000 |
parents | |
children |
rev | line source |
---|---|
wolffd@0 | 1 function lam_msg = CPD_to_lambda_msg(CPD, msg_type, n, ps, msg, p, evidence) |
wolffd@0 | 2 % CPD_TO_LAMBDA_MSG Compute lambda message (discrete) |
wolffd@0 | 3 % lam_msg = compute_lambda_msg(CPD, msg_type, n, ps, msg, p, evidence) |
wolffd@0 | 4 % Pearl p183 eq 4.52 |
wolffd@0 | 5 |
wolffd@0 | 6 switch msg_type |
wolffd@0 | 7 case 'd', |
wolffd@0 | 8 T = prod_CPT_and_pi_msgs(CPD, n, ps, msg, p); |
wolffd@0 | 9 mysize = length(msg{n}.lambda); |
wolffd@0 | 10 lambda = dpot(n, mysize, msg{n}.lambda); |
wolffd@0 | 11 T = multiply_by_pot(T, lambda); |
wolffd@0 | 12 lam_msg = pot_to_marginal(marginalize_pot(T, p)); |
wolffd@0 | 13 lam_msg = lam_msg.T; |
wolffd@0 | 14 case 'g', |
wolffd@0 | 15 error('discrete_CPD can''t create Gaussian msgs') |
wolffd@0 | 16 end |