Prediction using sample_ppc in Hierarchical model

Unfortunately it might be an issue with Multinomial, see eg this issue: https://github.com/pymc-devs/pymc3/issues/2614