Thank you @ricardoV94 for your proposed solution. I tried it but sadly still getting errors.
would the built in ‘shape’ args resolve this? like can I use it as below:
for control_var in [var for var in df.columns if var.startswith('cont_var')]:
x = df[control_var].values
control_beta = pm.Normal(f"{control_var}_control_coef",
mu=cont_beta_mu,
sigma=cont_beta_sigma,
dims= 'Region',
--> shape= (120,)
)
control_x = control_beta * x
print(control_x.shape)