2017-02-10 13 views
0

アルファ= 29(スケール)とベータ= 3(サイズ)のガンマ分布をプロットします。言い換えれば、ガンマ(29,3)のpdfをプロットしたいと思います。 documentationによれば、pythonガンマ関数はパラメータaとxのみを持ち、sizeパラメータが存在しない場合はどうすればいいですか?pythonでalphaとbetaパラメータを使ってガンマ分布をプロットする方法

私はドキュメントによるとlocはベータ版だったが、私はそれが実際にオフセットだと思うので、以下のコードが間違っている...

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

x = np.linspace (0, 100, 200) 
y1 = stats.gamma.pdf(x, a=29, loc=3) #a is alpha, loc is beta??? 
plt.plot(x, y1, "y-", label=(r'$\alpha=29, \beta=3$')) 


plt.ylim([0,0.08]) 
plt.xlim([0,150]) 
plt.show() 

答えて

3

思った、あなたはスケールパラメータ(シータ)を使用したいのですが、シータの逆であるベータを定義しているので、1 /βの値としてスケールを渡します。この例では、1/3または0.33333になります。

ので、試してみてください。

y1 = stats.gamma.pdf(x, a=29, scale=0.33333) 
関連する問題