numpy 1.9.0: ValueError: probabilities do not sum to 1
I think 1.7e-6 is a large enough relative error to be worth complaining about. You can renormalize easily enough, though, if you're confident the error is negligible:
>>> probs = np.array(probs)>>> probs /= probs.sum()>>> probs.sum()1.0>>> samples = np.random.choice(arr, size=1000, replace=True, p=probs)>>> samples[:5]array([ 1.37635054, 1.1287515 , 1.7229892 , 19.8967587 , 2.07953181])