Help with Poisson model

Yegads, well, if you want to figure out how to fix that C++ there’s a very small chance that we’d see an interesting behavior :P. Long shot, but do you have the latest Rstan?

I’m honestly stumped on this. Without the specific data/model causing the problem, it’s hard to do anything other than make undirected guesses.

If it’s absolutely commenting and uncommenting this that causes the problem:

for (it in 1:I * T) {
  y_sim[it] = poisson_log_rng(y_hat[it]);
}

then I’m stumped! It doesn’t look like you’re using any fixed seeds or anything?