wolffd@0: function pi = CPD_to_pi(CPD, msg_type, n, ps, msg, evidence)
wolffd@0: % CPD_TO_PI Compute pi vector (noisyor)
wolffd@0: % pi = CPD_to_pi(CPD, msg_type, n, ps, msg)
wolffd@0: % Pearl p188 eqn 4.57
wolffd@0:   
wolffd@0: switch msg_type
wolffd@0:  case 'd',
wolffd@0:    pi = sum_prod_CPD_and_pi_msgs(CPD, n, ps, msg);
wolffd@0:    pi = [pi 1-pi]';
wolffd@0:  case 'g', 
wolffd@0:   error('can''t convert noisy-or CPD to Gaussian pi')
wolffd@0: end