Hi all -

As you may know, I have an R package that does time-varying latent variable models with Stan (GitHub - saudiwin/idealstan: idealstan offers item-response theory (IRT) ideal-point estimation for binary, ordinal, counts and continuous responses with time-varying and missing-data inference. Latent space model also included. Full and approximate Bayesian sampling with 'Stan' (www.mc-stan.org).). So I’ve seen a lot of the problems with identifiability. However, there is something new that I’m picking up that I hadn’t noticed before, and that is convergence to modes with different log-posterior values. For example, I just ran a model with two chains that ended up with the step sizes/treedepth values in the attached files.

As can be seen, one chain converged to a more likely mode with a *lower* treedepth and *higher* stepsize. The other chain found a less likely mode that required longer transitions.

Here’s the question–when dealing with this kind of modality (i.e., the model has an invariance to rotations problem), does it make sense to always select the chain with the higher log-likelihood, especially if the treedepth is lower? It would seem that in this case one chain is clearly superior to the other.