Kanban lead time weibull distribution estimator

weibull_estimate.py

import numpy as np

def weibull_estimate(min_, max_, shape=1.25):
    """Draw a random number between min_ and max_ using a Weibull
    distribution with shape `shape`.

    Thanks to Troy Magennis for the inspiration!
    """

    return min(max_, ((max_ - min_) / 4.0) * np.random.weibull(shape) + min_)

# e.g. to plot a histogram:

import matplotlib.pyplot as plt

plt.hist([weibull_estimate(5, 20) for x in range(1000)])