paired_students_t_test.py

from scipy.stats import ttest_rel

# Randomly generate the data
x1 = rng.normal(loc=0.00, scale=1, size=100)
x2 = x1 + rng.normal(loc=0.25, scale=1, size=100)

# Calculate test statistic and p-value
stat, p = ttest_rel(x1, x2)

# Interpreation
print('stat=%.3f, p=%.3f' % (stat, p))
if p > 0.05:
    print('Do not reject the null hypothesis and conclude the means of the paired samples are the same.')
else:
    print('Reject the null hypothesis and conclude the means of the paired samples are not the same.')