How to properly do out-of-sample prediction for hierarchical model

[peeeffchang]

Update: I think I’ve seen someone else making the same observation, that it turns out sample_posterior_predictive works only if the number of out-of-sample forecast is the same as training set.