私は scipy.distributions パッケージの使用方法を理解するのに苦労しており、誰かが私のためにサンプルコードを投稿できるかどうか疑問に思いました. 必要なものはすべて揃っているようですが、使い方がわかりません。
対数正規分布とポアソン分布の 2 つの分布を生成する必要があります。それぞれの分散とラムダを知っています。
リソースへのリンクも同様に機能します。
私はあなたがの分布を意味すると仮定しますscipy.stats
。分布を作成するには、ランダムな変量を生成し、pdfを計算します。
Python 2.5.1(r251:54863、2008年2月4日、21:48:13)[GCC 4.0.1(Apple Inc.ビルド5465)]ダーウィンタイプ「help」、「copyright」、「credits」、または「license」詳細については。
>>> 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の確率分布