def _vector_of_log_likelihoods(theta_1, theta_2, x): lik = np.log(0.5 * norm.pdf(x, theta_1, SIGMA_x) + 0.5 * norm.pdf(x, theta_1 + theta_2, SIGMA_x)) return lik