Mercurial > hg > camir-ismir2012
annotate toolboxes/FullBNT-1.0.7/bnt/CPDs/@tabular_CPD/log_prior.m @ 0:cc4b1211e677 tip
initial commit to HG from
Changeset:
646 (e263d8a21543) added further path and more save "camirversion.m"
| author | Daniel Wolff |
|---|---|
| date | Fri, 19 Aug 2016 13:07:06 +0200 |
| parents | |
| children |
| rev | line source |
|---|---|
| Daniel@0 | 1 function L = log_prior(CPD) |
| Daniel@0 | 2 % LOG_PRIOR Return log P(theta) for a tabular CPD |
| Daniel@0 | 3 % L = log_prior(CPD) |
| Daniel@0 | 4 |
| Daniel@0 | 5 switch CPD.prior_type |
| Daniel@0 | 6 case 'none', |
| Daniel@0 | 7 L = 0; |
| Daniel@0 | 8 case 'dirichlet', |
| Daniel@0 | 9 D = CPD.dirichlet(:); |
| Daniel@0 | 10 L = sum(log(D + (D==0))); |
| Daniel@0 | 11 case 'entropic', |
| Daniel@0 | 12 % log-prior = log exp(-H(theta)) = sum_i theta_i log (theta_i) |
| Daniel@0 | 13 fam_sz = CPD.sizes; |
| Daniel@0 | 14 psz = prod(fam_sz(1:end-1)); |
| Daniel@0 | 15 ssz = fam_sz(end); |
| Daniel@0 | 16 C = reshape(CPD.CPT, psz, ssz); |
| Daniel@0 | 17 L = sum(sum(C .* log(C + (C==0)))); |
| Daniel@0 | 18 end |