Мне нужно нарисовать кривую плотности на гистограмме с фактической высотой столбцов (фактической частотой) в качестве оси Y.
Try1:
Я нашел соответствующий ответ здесь , но он нормализовал гистограмму к диапазону кривой.
Ниже мой код и вывод.
import numpy as np
import matplotlib.mlab as mlab
import matplotlib.pyplot as plt
from scipy.stats import norm
data = [125.36, 126.66, 130.28, 133.74, 126.92, 120.85, 119.42, 128.61, 123.53, 130.15, 126.02, 116.65, 125.24, 126.84,
125.95, 114.41, 138.62, 127.4, 127.59, 123.57, 133.76, 124.6, 113.48, 128.6, 121.04, 119.42, 120.83, 136.53, 120.4,
136.58, 121.73, 132.72, 109.25, 125.42, 117.67, 124.01, 118.74, 128.99, 131.11, 112.27, 118.76, 119.15, 122.42,
122.22, 134.71, 126.22, 130.33, 120.52, 126.88, 117.4]
(mu, sigma) = norm.fit(data)
x = np.linspace(min(data), max(data), 100)
plt.hist(data, bins=12, normed=True)
plt.plot(x, mlab.normpdf(x, mu, sigma))
plt.show()
Try2:
Там @DavidG предоставил опцию, определяемую пользователем функцию, даже если она точно не покрывает плотность гистограммы.
def gauss_function(x, a, x0, sigma):
return a * np.exp(-(x - x0) ** 2 / (2 * sigma ** 2))
test = gauss_function(x, max(data), mu, sigma)
plt.hist(data, bins=12)
plt.plot(x, test)
plt.show()
Результат для этого был,
Но фактическая гистограмма ниже, где ось Y колеблется от 0 до 8,
И я хочу нарисовать кривую плотности именно на этом. Будем благодарны за любую помощь в этом отношении.