Thanks for the quick reply. Here is an example that gives the error:
y=1
mu=0.000519101
phi=1.70972
I get the error about half the times I fit my data (less often with the simulated data above). As soon as I get the error, that chain crashes.
Interestingly, I don’t get it with BRMS, which can also fit this model. I guess the priors used there avoid it.