Gerador Aleatório para Distribuição de Poisson

É uma distribuição de probabilidade discreta que expressa a probabilidade de um determinado número de eventos ocorrer em um intervalo fixo de tempo ou espaço se esses eventos ocorrerem com uma taxa média constante conhecida e independentemente do tempo desde o último evento.

Exemplo em Python

Exemplo com 1 milhão de sorteios

def pseudo_poisson(alpha, size=1):
    """
    Poisson generator from uniform generator
    """
    poisson = []
    for _ in range(size):
    
        # Sets seed based on the decimal portion of the current system clock
        t = time.perf_counter()
        seed = int(10**9*float(str(t-int(t))[0:]))
        
        U = pseudo_uniform(seed=seed, size=5*alpha)
        X,P,i = 0,1,0
        while P >= np.exp(-alpha):
            P = U[i]*P
            X += 1
            i += 1
        
        poisson.append(X)
    
    return np.array(poisson)