Bambi vs pymc model specification

Questions
2

I think it’s a broadcasting issue with deer_y. Can you make the modification

# Gaussian Likelihood
    y_pred = pm.Normal('y_pred',
                       mu=μ,
                       sd=ϵ,
                       observed=deer_y.squeeze())

and tell me if you get results that line up with Bambi?

1 Like