Error occurred during calling the sampler; sampling not done

[increasechief]

Is it typical to have the rng inside of your .stan file this way?

generated quantities{
  vector[N] sim_Y;
  for (n in 1:N)  sim_Y[n] = normal_rng(mean(xb), sig)';
}

Specifying the number of samples for rng

But this brings up a more general question of where to do posterior summarization. It’s a mix. Do it wherever it is easiest/makes the most sense. But usually if you already wrote the model in Stan it’s easy to just copy-paste a few things and compute whatever you want in the generated quantities. generated quantities{ vector[10] y_pred; real E_y_pred; // Expectation of posterior prediction for (i in 1:10) y_pred[i] = normal_rng(beta * x_new[i], sigma); E_y_pred = mean(y_pred);…