有没有人有使用scipy.stats.distributions的示例代码？

Question

我在努力弄清楚如何使用scipy.distributions包，并想知道如果有人能为我张贴一些示例代码。这似乎做我需要的一切，我只是无法弄清楚如何使用它。

我需要生成两个分布，一个对数正态和一个泊松。我知道每个的方差和λ

链接资源将很好的工作。

Answer 1

我假定你的意思是在scipy.stats分布。要创建的分布，生成随机变元，并计算pdf文件：

的Python 2.5.1（R251：54863，2008年2月4日，21时48分13秒） [4.0.1 GCC（苹果公司建立5465）]对达尔文 键入“help”，“版权”，“信用”或“许可”的详细信息。

>>> from scipy.stats import poisson, lognorm
>>> myShape = 5;myMu=10
>>> ln = lognorm(myShape)
>>> p = poisson(myMu)
>>> ln.rvs((10,)) #generate 10 RVs from ln
array([  2.09164812e+00,   3.29062874e-01,   1.22453941e-03,
         3.80101527e+02,   7.67464002e-02,   2.53530952e+01,
         1.41850880e+03,   8.36347923e+03,   8.69209870e+03,
         1.64317413e-01])
>>> p.rvs((10,)) #generate 10 RVs from p
array([ 8,  9,  7, 12,  6, 13, 11, 11, 10,  8])
>>> ln.pdf(3) #lognorm PDF at x=3
array(0.02596183475208955)

的其它方法（以及scipy.stats文档的其余部分）可以在新的 SciPy的文档站点。

Answer 2

下面是一些示例代码：概率分布在SciPy的