tips_grouped = tips_all.groupby(['smoker', 'size'])
tips = tips_grouped.mean()
tips['CI'] = tips_grouped.total_bill.apply(stats.sem) * 1.96
tips.reset_index(inplace=True)
Plotting errors bars from dataframe using Seaborn FacetGrid