Hi @linas! The model that I posted (zoib.stan) also has \mu varying for each observation.
If that doesn’t solve your problem, could you please clarify what your intended model is (e.g., using math notation), so we can get a better sense for whether your Stan model matches your intended model?