kruskal_wallis_h_test.py

from scipy.stats import kruskal

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

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

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