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.