Pylab，绘制概率分布

Question

我正在使用pylab进行科学计算。我用hist方法制作了一个直方图。不过，我也想画出一个概率分布。有没有人知道怎么做？

Answer 1

Scipy里有很多概率分布的功能，这些功能都是现成的，可以直接使用。这里有一些文档里的例子：

norm: 正态分布的连续随机变量。

alpha: 阿尔法分布的连续随机变量。

beta: 贝塔分布的连续随机变量。

chi: 卡方分布的连续随机变量。

erlang: 埃尔朗分布的连续随机变量。

gamma: 伽马分布的连续随机变量。

...

接下来是一些关于如何从这些分布中抽样并绘制 norm, beta, chi 和 alpha 的例子：

from scipy.stats import *
from pylab import hist, show,subplot

N = 2000

subplot(221)
X = norm.rvs(size=N)
hist(X,normed=True,bins=30)

subplot(222)
X = beta.rvs(2,1,size=N)
hist(X,normed=True,bins=30)

subplot(223)
X = chi.rvs(2.0, size=N)
hist(X,normed=True,bins=30)

subplot(224)
X = gamma.rvs(1.0, size=N)
hist(X,normed=True,bins=30)
show()