I think the issue is that \theta \not\perp n in my problem. When I introduce n as a regressor using Bernoulli likelihood, I find a huge coefficient.
If I have time, I may go back and confirm whether the difference persists with Binomial likelihood when including n as a regressor. But I think I have a grip on this now; thanks for your feedback!