dossier-stats/mediane.py at main · simjnd/dossier-stats · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
import matplotlib.pyplot as plt
import numpy as np

t = [ 162.5, 167.5, 172.5, 172.5, 172.5, 172.5, 172.5, 172.5, 172.5, 172.5, 172.5, 172.5, 172.5, 177.5, 177.5, 177.5, 177.5, 177.5, 177.5, 177.5, 177.5, 177.5, 177.5, 177.5, 177.5, 177.5, 182.5, 182.5, 182.5, 182.5, 182.5]
p = [ 47.5, 57.5, 52.5, 52.5, 52.5, 52.5, 57.5, 62.5, 62.5, 62.5, 67.5, 72.5, 72.5, 52.5, 57.5, 62.5, 62.5, 67.5, 67.5, 67.5, 72.5, 72.5, 72.5, 72.5, 77.5, 92.5, 67.5, 72.5, 77.5, 77.5, 97.5]

print("-------------TAILLE--------------------------------------------------------------")
e = t
avg = np.average(e)
print(f"Moyenne: {round(avg, 4)} cm")
q2 = np.median(e)
print(f"Médiane: {q2} cm")
q1 = np.quantile(e, .25)
print(f"Premier Quartile: {q1} cm")
q3 = np.quantile(e, .75)
print(f"Troisième Quartile: {q3} cm")
print(f"Variance: {round(np.var(e), 4)}")
print(f"Ecart-type: {round(np.std(e), 4)} cm")
print("")
eu = np.unique(e)
print("Fréquence")
for ei in eu:
    occ = np.count_nonzero(e == ei)
    print(f"Taille: {ei}, Occ: {occ}, Freq: {round((occ / len(e)) * 100, 4)}%")
print("")
print(f"Etendue: {np.sort(e)[-1] - np.sort(e)[0]} cm")
print(f"Intervalle Inter-Quartile: {q3 - q1}")
ecart_absolu_moyen = 0
for ei in e:
    ecart_absolu_moyen += np.abs(avg - ei)
ecart_absolu_moyen = round(ecart_absolu_moyen / len(e), 4)
print(f"Ecart absolu moyen: {ecart_absolu_moyen}")
ecart_absolu_median = 0
for ei in e:
    ecart_absolu_median += np.abs(np.median(e) - ei)
ecart_absolu_median = round(ecart_absolu_median / len(e), 4)
print(f"Ecart absolu median: {ecart_absolu_median}")

print(f"Coefficient de variation: {round(np.std(e) / np.average(e), 4)}")
print(f"Ecart moyen relatif: {round(ecart_absolu_moyen / np.average(e), 4)}")
print(f"Coefficient Interquartile Relatif: {round((q3 - q1) / q2, 4)}")

print("-------------POIDS--------------------------------------------------------------")
e = p
avg = np.average(e)
print(f"Moyenne: {round(avg, 4)} kg")
q2 = np.median(e)
print(f"Médiane: {q2} kg")
q1 = np.quantile(e, .25)
print(f"Premier Quartile: {q1} kg")
q3 = np.quantile(e, .75)
print(f"Troisième Quartile: {q3} kg")
print(f"Variance: {round(np.var(e), 4)}")
print(f"Ecart-type: {round(np.std(e), 4)} kg")
print("")
eu = np.unique(e)
print("Fréquence")
for ei in eu:
    occ = np.count_nonzero(e == ei)
    print(f"Taille: {ei}, Occ: {occ}, Freq: {round((occ / len(e)) * 100, 4)}%")
print("")
print(f"Etendue: {np.sort(e)[-1] - np.sort(e)[0]} kg")
print(f"Intervalle Inter-Quartile: {q3 - q1}")
ecart_absolu_moyen = 0
for ei in e:
    ecart_absolu_moyen += np.abs(avg - ei)
ecart_absolu_moyen = round(ecart_absolu_moyen / len(e), 4)
print(f"Ecart absolu moyen: {ecart_absolu_moyen}")
ecart_absolu_median = 0
for ei in e:
    ecart_absolu_median += np.abs(np.median(e) - ei)
ecart_absolu_median = round(ecart_absolu_median / len(e), 4)
print(f"Ecart absolu median: {ecart_absolu_median}")

print(f"Coefficient de variation: {round(np.std(e) / np.average(e), 4)}")
print(f"Ecart moyen relatif: {round(ecart_absolu_moyen / np.average(e), 4)}")
print(f"Coefficient Interquartile Relatif: {round((q3 - q1) / q2, 4)}")