@@ -93,7 +93,7 @@ def neg_post_phi_mono(phi, X, y):
# Beta(alpha, beta) prior on p → negative log-prior
npr_p = -((alpha - 1) * np.log(p + eps) + (beta - 1) * np.log(1 - p + eps))
- return nll + npr_r
+ return nll + npr_p
# Initialization (stable, simple)
def init_phi(X, y):
