plt.hlines(xmin='lower', xmax='upper', y='y', data=data, 
           linewidth=5, color='steelblue', alpha=0.5)

# Point-estimate for reference
plt.plot('est', 'y', 'k|', data = data)