有没有人有使用scipy.stats.distributions的示例代码?
题
我在努力弄清楚如何使用scipy.distributions包,并想知道如果有人能为我张贴一些示例代码。这似乎做我需要的一切,我只是无法弄清楚如何使用它。
我需要生成两个分布,一个对数正态和一个泊松。我知道每个的方差和λ
链接资源将很好的工作。
解决方案
我假定你的意思是在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的文档站点。
其他提示
下面是一些示例代码:概率分布在SciPy的
不隶属于 StackOverflow