In addition, you can make it work by supplying a value to the graph:
In [13]: with pm.Model() as model:
...: beta_soc_dem_edu_spend = pm.Normal('β(soc_dem, edu spend)', shape=1
...: )
...: print('β(soc_dem, edu spend) shape is: %s'%beta_soc_dem_edu_spend.s
...: hape)
...: print('β(soc_dem, edu spend) shape is: %s'%beta_soc_dem_edu_spend.s
...: hape.eval({beta_soc_dem_edu_spend: np.zeros(1)}))