How to account for multiple betas in example documentation?

Thanks, super helpful! Could you recommend a procedure if the betas were all chained together, say:

with pm.Model() as model_1:
    a = pm.Normal("a", 0.0, 10.0)
    betas = pm.Normal('betas', mu=[40, 25, 55], sigma=5, dims='features')
    ...

Or should we stick to making discrete entries for each parameter (to make our lives easier later?)

Also in your example, where does idata come into play; it seems as an unused variable.