---

import re, numpy as np
from scipy import stats
import matplotlib.pyplot as plt
from collections import Counter

s = """сюда вставь весь список"""

x = re.findall(r'NA|Jan|Mar|Apr|May|Jun|Jul|Aug|Sept|Oct|Nov', s)
x = [i for i in x if i != 'NA']

cnt = Counter(x)
obs = np.array(list(cnt.values()))
k = len(cnt)
n = len(x)
exp = n / k

chi = sum((obs - exp) ** 2 / exp)
crit = stats.chi2.ppf(1 - 0.03, k - 1)

print(k - 1)
print(f"{chi:.8f}")
print(f"{crit:.8f}")
print("Да" if chi > crit else "Нет")

plt.bar(cnt.keys(), cnt.values())
plt.show()