Try something like mu = pm.Normal('mu', alpha + beta_1*x['feature_1_2'].values, sd=1.)
mu = pm.Normal('mu', alpha + beta_1*x['feature_1_2'].values, sd=1.)