Sample_ppc shape error with multi-index model

No problem! I filed a bug on Github: https://github.com/pymc-devs/pymc3/issues/3147

I like the 1d vector set up better, because it is also easier if you have missing data.