docs: def discrete_random_draw(data, nb=1):
      ''' Code from Steve Nguyen'''
      data = np.array(data)
      if not data.any():
          data = np.ones_like(data)
      data = data/data.sum()
      xk = np.arange(len(data))
      custm = stats.rv_discrete(name='custm', values=(xk, data))
      return custm.rvs(size=nb)