# Set lower and upper bounds of ribbon
plt.fill_between(x='day', y1='lower', y2='upper', data=cinci_so2)

# Add point-estimate reference line
plt.plot('day', 'mean', 'w-', alpha=0.5, data=data)