Mercurial > hg > camir-aes2014
annotate toolboxes/FullBNT-1.0.7/bnt/potentials/@cgpot/normalize_pot.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 [pot, loglik] = normalize_pot(pot) |
| wolffd@0 | 2 % NORMALIZE_POT Convert the CG potential Pr(X,E) into Pr(X|E) and return log Pr(E). |
| wolffd@0 | 3 % [pot, loglik] = normalize_pot(pot) |
| wolffd@0 | 4 |
| wolffd@0 | 5 % Marginalize down to [], so that the normalizing constant becomes Pr(E) |
| wolffd@0 | 6 temp = cg_can_to_mom(marginalize_pot(pot, [])); |
| wolffd@0 | 7 %loglik = temp.mom{1}.logp; |
| wolffd@0 | 8 [temp2, loglik] = normalize_pot(temp.mom{1}); |
| wolffd@0 | 9 |
| wolffd@0 | 10 % Adjust scale factor to reflect the fact that the pot now represents Pr(X | E) instead of Pr(X,E). |
| wolffd@0 | 11 |
| wolffd@0 | 12 scale = -loglik; |
| wolffd@0 | 13 if 1 |
| wolffd@0 | 14 switch pot.subtype |
| wolffd@0 | 15 case 'm' |
| wolffd@0 | 16 for i=1:pot.dsize |
| wolffd@0 | 17 pot.mom{i} = rescale_pot(pot.mom{i}, scale); |
| wolffd@0 | 18 end |
| wolffd@0 | 19 case 'c' |
| wolffd@0 | 20 for i=1:pot.dsize |
| wolffd@0 | 21 pot.can{i} = rescale_pot(pot.can{i}, scale); |
| wolffd@0 | 22 end |
| wolffd@0 | 23 end |
| wolffd@0 | 24 end |