Dimension error for Hierarchical model

The issue seems to be this line:

control_x = control_beta * x

because control_beta.shape==5 (one per region) and x.shape==(120,).