python - matplotlib.mlab.normpdf() 的正确用法是什么？

Question

我打算让我正在编写的程序的一部分在多个原始文本源上自动生成各种统计数据的高斯分布，但是我在按照以下指南生成图表时遇到了一些问题:

python pylab plot normal distribution

剧情代码的大致要点如下。

import numpy as np
import matplotlib.mlab as mlab
import matplotlib.pyplot as pyplot

meanAverage = 222.89219487179491    # typical value calculated beforehand
standardDeviation = 3.8857889432054091    # typical value calculated beforehand

x = np.linspace(-3,3,100)
pyplot.plot(x,mlab.normpdf(x,meanAverage,standardDeviation))
pyplot.show()

它所做的只是生成一条看起来相当平坦且无用的 y = 0 行!谁能看出这里的问题是什么？

干杯。

Answer 1

如果您阅读 matplotlib.mlab.normpdf 的文档，此功能已取消，您应该使用 scipy.stats.norm.pdf相反。

Deprecated since version 2.2: scipy.stats.norm.pdf

因为您的分布均值约为 222，所以您应该使用 np.linspace(200, 220, 100)。

因此您的代码将如下所示:

import numpy as np
from scipy.stats import norm
import matplotlib.pyplot as pyplot

meanAverage = 222.89219487179491    # typical value calculated beforehand
standardDeviation = 3.8857889432054091    # typical value calculated beforehand

x = np.linspace(200, 220, 100)
pyplot.plot(x, norm.pdf(x, meanAverage, standardDeviation))
pyplot.show()