Mercurial > hg > camir-aes2014
view toolboxes/FullBNT-1.0.7/bnt/CPDs/@tabular_CPD/Old/BIC_score_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 score = BIC_score_CPD(CPD, fam, data, ns, cnodes) % BIC_score_CPD Compute the BIC score of a tabular CPD % score = BIC_score_CPD(CPD, fam, data, ns, cnodes) if iscell(data) local_data = cell2num(data(fam,:)); else local_data = data(fam, :); end counts = compute_counts(local_data, CPD.sizes); CPT = mk_stochastic(counts); % MLE tiny = exp(-700); CPT = CPT + (CPT==0)*tiny; % replace 0s by tiny LL = sum(log(CPT(:)) .* counts(:)); N = size(data, 2); score = LL - 0.5*CPD.nparams*log(N);