Statistical Thinking in Python (Part 1)
Justin Bois
Teaching Professor at the California Institute of Technology

samples = rng.poisson(6, size=10000)
x, y = ecdf(samples)
_ = plt.plot(x, y, marker='.', linestyle='none')
plt.margins(0.02)
_ = plt.xlabel('number of successes')
_ = plt.ylabel('CDF')
plt.show()

Statistical Thinking in Python (Part 1)