benford.png

benford.py

from collections import Counter
import numpy as np
import matplotlib.pyplot as plt

def generer_nombres_avec_bornes(bornes, nb_valeurs_par_borne=5):
    """Renvoie un array qui rassemble des nombres entre 0 et une borne. Pour chaque borne, la fonction génère <len(bornes)> * nb_valeurs_par_borne (par défaut 5) valeurs"""
    return np.concatenate([np.random.randint(borne+1, size=nb_valeurs_par_borne) for borne in bornes])

def dessiner_avec_echantillon(echantillon):
    """À partir d'un échantillon de nombres, dessine le diagramme en bâtons correspondant au premier chiffre de chaque nombre"""
    counter = Counter(str(val)[0] for val in echantillon)
    taille = len(echantillon)
    plt.bar(range(1, 10), [counter[str(c)]/taille for c in range(1,10)])
    plt.show()
    
    

#dessiner_avec_echantillon(generer_nombres_avec_bornes(np.random.randint(1,10000000, size=100)))
#dessiner_avec_echantillon(generer_nombres_avec_bornes(generer_nombres_avec_bornes(np.random.randint(1,10000000, size=100))))
#dessiner_avec_echantillon(generer_nombres_avec_bornes(generer_nombres_avec_bornes(generer_nombres_avec_bornes(np.random.randint(1,10000000, size=100)))))
dessiner_avec_echantillon(generer_nombres_avec_bornes(generer_nombres_avec_bornes(generer_nombres_avec_bornes(generer_nombres_avec_bornes(np.random.randint(1,10000000, size=1000))))))