I encounter a job that involve a very large hierarchical model (around 40k parameters and 1m samples). We used to code them in rJAGS and it took long to converge. We recently recoded it in rstan, it seems that it is doing worse in stan.
Is it ture that stan preferably to fit a comparatively small model (less p) and not work well in a very big one?
Thanks,
TJ
No, that is not true in general. Post your code and you might get some pointers as to what might be awry.
Exactly the opposite if you look at relative performance compared to other tools.
You want to compare effective sample size adjusting for convergence, not just iteration speed.