t_stat, p_t = stats.ttest_ind(grupo_A, grupo_B)
data = df['total_bill']
mean = np.mean(data)
sem = stats.sem(data) # standard error of mean
ci = stats.t.interval(0.95, len(data)-1, loc=mean, scale=sem)
print(f"95% CI: ci")
stats.ttest_rel(before_measurement, after_measurement)