Shape error when making out-of-sample predictions

Try to add shape=output.shape to the observed variable. You have to tell PyMC how the shape depends on the parameters if you want to only change the predictors.