Участок нормальной кривой не ... нормальный

Question

Я строю график распределения результатов теста и подобранной кривой к этим оценкам:

h = sorted(data['Baseline'])  #sorted
fit = stats.norm.pdf(h, np.mean(h), np.std(h))
plt.plot(h,fit,'-o')
plt.hist(h,normed=True)      #use this to draw histogram of your data
plt.show()

Тем не менее, график pdf не выглядит нормальным (см. Излом на кривой около x = 60). Смотрите вывод:

Я не уверен, что здесь происходит ... любая помощь приветствуется. Это потому, что между поставленными наблюдениями проводится нормальная линия? Может предоставить вам реальные данные, если это необходимо, есть только 60 наблюдений.

Answer 1

Да, вы оцениваете норма-pdf по сохранности. Вместо этого вы хотели бы создать некоторые другие данные, такие как

h = sorted(data['Baseline'])  #sorted
x = np.linspace(h.min(), h.max(), 151)

fit = stats.norm.pdf(x, np.mean(h), np.std(h))

plt.plot(x,fit,'-')
plt.hist(h,normed=True) 
plt.show()

Обратите внимание, что данные не выглядят нормально распределенными вообще. Поэтому потенциально вы бы предпочли использовать другое распределение или, возможно, выполнить оценку плотности ядра.