fig = plt.figure(figsize=(10, 7))
bins = plt.hist(stats[:, 0], bins=20)
centers = (bins[1][1:]+bins[1][:-1])/2
plt.xlabel('Frequencey of Max Amplitude', fontsize=20)