Mathematics.
compared to those on the lower end. from scipy.stats import skewnorm a = 3 mean, std = 70, 140 fig, ax = plt.subplots(1, 1) x = np.linspace(0,600, 200) ax.plot(x, skewnorm.pdf(x, a, mean, std),'r-', lw=5) r = skewnorm.rvs(a, mean, std, size=5000) ax.hist(r, density=True, histtype='stepfilled',