Use to rstan and rethinking pacakge. trying to switch to pymc3...getting an error

Questions
7

No - actually I should write:

mu = pm.Normal('mu', alpha + beta_1*x['feature_1_2'].values, sd=1., model_output.shape)

Because it should have the same shape as the observed (it represents the latent mean of the observation)