Non-repetitive random number in numpy

by

numpy.random.Generator.choice offers a replace argument to sample without replacement:

from numpy.random import default_rng

rng = default_rng()
numbers = rng.choice(20, size=10, replace=False)

If you’re on a pre-1.17 NumPy, without the Generator API, you can use random.sample() from the standard library:

print(random.sample(range(20), 10))

You can also use numpy.random.shuffle() and slicing, but this will be less efficient:

a = numpy.arange(20)
numpy.random.shuffle(a)
print a[:10]

There’s also a replace argument in the legacy numpy.random.choice function, but this argument was implemented inefficiently and then left inefficient due to random number stream stability guarantees, so its use isn’t recommended. (It basically does the shuffle-and-slice thing internally.)

Some timings:

import timeit
print("when output size/k is large, np.random.default_rng().choice() is far far quicker, even when including time taken to create np.random.default_rng()")
print(1, timeit.timeit("rng.choice(a=10**5, size=10**4, replace=False, shuffle=False)", setup="import numpy as np; rng=np.random.default_rng()", number=10**3)) #0.16003450006246567
print(2, timeit.timeit("np.random.default_rng().choice(a=10**5, size=10**4, replace=False, shuffle=False)", setup="import numpy as np", number=10**3)) #0.19915290002245456

print(3, timeit.timeit("random.sample( population=range(10**5), k=10**4)", setup="import random", number=10**3))   #5.115292700007558

print("when output size/k is very small, random.sample() is quicker")
print(4, timeit.timeit("rng.choice(a=10**5, size=10**1, replace=False, shuffle=False)", setup="import numpy as np; rng=np.random.default_rng()", number=10**3))  #0.01609779999125749
print(5, timeit.timeit("random.sample( population=range(10**5), k=10**1)", setup="import random", number=10**3))  #0.008387799956835806

So numpy.random.Generator.choice is what you usually want to go for, except for very small output size/k.

More Related Contents:

  1. how to generate integer random number in fortran 90 in the range [0,5]?
  2. How to test random numbers?
  3. Java program that ask the user to enter the number of random values they wish to generate
  4. How to generate random binary numbers 0 or 1 with length of N and with option to control the probability of having 0 or 1?
  5. Is there a way to ordinate numbers in a list?
  6. Random / noise functions for GLSL
  7. Java Generate Random Number Between Two Given Values [duplicate]
  8. How to generate a random integer in the range [0,n] from a stream of random bits without wasting bits?
  9. How many double numbers are there between 0.0 and 1.0?
  10. How does a random number generator work?
  11. Generate a random number in a certain range in MATLAB
  12. C# Random Numbers aren’t being “random”
  13. Why are my random numbers always the same?
  14. How to specify upper and lower limits when using numpy.random.normal
  15. How to get a normal distribution within a range in numpy? [duplicate]
  16. Generate a random sample of points distributed on the surface of a unit sphere
  17. Generating uniform random numbers in Lua
  18. What is the best way to generate random numbers in C++?
  19. Creating same random number sequence in Python, NumPy and R
  20. Shuffling NumPy array along a given axis
  21. How to generate a random normal distribution of integers
  22. Generating non-repeating random numbers in Python
  23. Weighted random sample without replacement in python
  24. How do I create a random String by sampling from alphanumeric characters?
  25. Generate N random numbers whose sum is a constant K – Excel
  26. how does random() actually work?
  27. Vectorizing `numpy.random.choice` for given 2D array of probabilities along an axis
  28. Working with floating point NumPy arrays for comparison and related operations
  29. Generate random locations within a triangular domain
  30. Generate random array of floats between a range

Leave a Comment