Then, the only thing that it is left are the values that you are supplying as arguments to the RV, such as b_lower, b_upper, etc. You might be supplying a data type that messes with shape, such as a list. You can’t do that. It should be something like:
b_lower = 0.
low= pm.Normal('low', mu= b_lower, sigma= 0.001)