With time differences like that there are presumably dramatic differences in what the chains are doing. I’m guessing they don’t mix?
From the code you posted, two things:
Do a non-centered parameterization for indRE, see here.
Tighten the prior on tau. Maybe a half-normal or an exponential.
But since you’ve already run the code a few times you should look at the output and see if the differences in the inferences made by the fast/slow chains show any sort of application-specific behavior (like maybe the results are interpretable and they give you clues about what is going on).
Also you can do code blocks on discourse by surrounding code on each side by three backticks (or there’s a little code formatting block in the toolbar of the post). Something like this: ```write code here```
I was guessing that should be something related to mixing too.
I had run two chains, and as this happens, the first chain finished usually fast, and the second chain takes longer.
For instance, if I fixed 10000 iterations with 2000 of burnin and 8 of thin, both chains arrived at 3000 iterations at the same time, but then the second chain stopped for a while to update (sometimes until the first chain get 8000 iterations). After that, the time between 4000 to 5000, 5000 to 6000 and so on, increased. Because of that, I supposed that could be something with mixing but I was not sure about that.
Thank You also for the hints about the post. This was my first one, and I was not quite sure how I could do this.