1from scipy import stats
2ttest = stats.ttest_ind(mean_a, mean_b, equal_var = False)
3p_value = ttest[1]
4print('p_value: {:0.3f}'.format(p_value))
5
6# H0: mean_a - mean_b = 0
7# HA: mean_a - mean_b != 0
8# For one tail t-test divide the p-value by 2
9# For p-value <= 0.05: Reject H0