Code Python pour vérifier et représenter graphiquement s'il suit la loi de Benford
main.py
a=""" """
#Pendant ce temps, copiez la séquence de nombres d'Excel, etc.
import collections
import matplotlib.pyplot as plt
a=a.split("\n")
a = [x for x in a if x != '']
#Supprimer null
l_head=list(map(lambda x: str(x)[0], a))
#Obtenez le premier personnage
l_count=collections.Counter(l_head).most_common(9)
#9e à partir du haut car il peut contenir 0
#l_count.sort(key=lambda x: x[0])
#Trier
#Lorsqu'il est représenté graphiquement, il est trié dans l'ordre de 123456789 ou la fréquence d'apparition. Tous les 123 une fois retiré.
l_count
l_rate=list(map(lambda x: x*100 /sum(list(zip(*l_count))[1]), list(zip(*l_count))[1]))
#Calcul du ratio
l_index=list(zip(*l_count))[0]
l_rate
plt.xticks([0.15, 1.15, 2.15, 3.15, 4.15, 5.15, 6.15, 7.15, 8.15], l_index)
plt.bar(range(9), l_rate, color='g', width=0.3, label='Selected_data', align="center")
ben={1:30.1, 2:17.6, 3:12.5, 4:9.7, 5:7.9, 6:6.7, 7:5.8, 8:5.1, 9:4.6}
ben_label= list(map(lambda x: ben[int(x)], l_index))
plt.bar([0.3,1.3,2.3,3.3,4.3,5.3,6.3,7.3,8.3], ben_label, color='b', width=0.3, label='Benford_law', align="center")
plt.legend(bbox_to_anchor=(0.5, 1),loc=2)
plt.show()
Exemple:
Population par unité administrative au Japon
Superficie pour chaque unité administrative au Japon
Densité de la population par unité administrative au Japon
Nombre de ménages par unité administrative au Japon
Tous les chiffres reflétés ici
Les statistiques proviennent des chiffres du Bureau des statistiques https://www.stat.go.jp/data/index.html
Recommended Posts