
function [theta_prior] = prior_theta(X)

theta_prior = sum(X)/(sum(sum(X)));
