% You can assume n_plus and n_minus are never 0 together
%
function score = log_bayesian_model_score(n_plus, n_minus)

% total counts for each node
n = n_plus + n_minus;

%
% YOU MUST FILL IN THE SCORE COMPUTATION
% 
% You may use log_fact(n) below that is equal to log n! for n >= 0
% 

score = 


% -----------------------------
function [logfac] = log_fact(n) 

logfac = gammaln(n+1); 

