wolffd@0: function CPD = update_params(CPD, ev, counts)
wolffd@0: % UPDATE_PARAMS Update the Dirichlet pseudo counts and compute the new MAP param estimates (tabular)
wolffd@0: %
wolffd@0: % CPD = update_params(CPD, ev) uses the evidence on the family from a single case.
wolffd@0: %
wolffd@0: % CPD = update_params(CPD, [], counts) does a batch update using the specified suff. stats.
wolffd@0: 
wolffd@0: if nargin < 3
wolffd@0:   n = length(ev);
wolffd@0:   data = cat(1, ev{:}); % convert to a vector of scalars
wolffd@0:   counts = compute_counts(data(:)', 1:n, mysize(CPD.CPT));
wolffd@0: end
wolffd@0:   
wolffd@0: CPD.prior = CPD.prior + counts;
wolffd@0: CPD.CPT = mk_stochastic(CPD.prior);