K-fold CV with observation-level random effects

[avehtari]

@mhollanders

Or do you predict a new random effect for each observation with a _rng() call or something?

Yes. You assume that you are predicting for a site from which you did not have observations, and thus you draw the site-specific parameters (aka “random effects”) from the prior using _rng(), which corresponds to integrating over them.