2016-08-18 18 views
4

私はscipyのダウンロードから、フォン・ミーゼス分布に合うようにしようとしている(http://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.vonmises.htmlPythonとScipy:von misesの配布方法は?

だから、私は、データをフィットする場合それは、しかし

enter image description here

を返し

from scipy.stats import vonmises 
kappa = 3 
r = vonmises.rvs(kappa, size=1000) 
plt.hist(r, normed=True,alpha=0.2) 

を試してみましたその上に

vonmises.fit(r) 
# returns (1.2222011312461918, 0.024913780423670054, 2.4243546157480105e-30) 

vonmises.fit(r, loc=0, scale=1) 
# returns (1.549290021706847, 0.0013319431181202394, 7.1653626652619939e-29) 

しかし、返される値のいずれもVon Mises、kappaのパラメータではありません。

戻り値は何ですか? Ifeelは2番目の値はlo​​cまたは平均値です。しかし、最初に返される値が何であるか分かりません。

フォンミーゼスディストリビューションにはどのようにしたらいいですか?

+0

戻り値は、順番に 'kappa'、' loc'、 'scale'です。 – dbliss

答えて

7

戻り値は、κ、loc、およびscaleです。残念ながら、von Mises pdfはあまり適合しないようです。スケールを固定すると正しくフィットします。

>>> vonmises.fit(r, fscale=1) 
(2.994517240859579, -0.0080482378119089287, 1) 
+0

保存した私の一日!これはドキュメントに入れる必要があります! – cqcn1991

+0

@ cqcn1991その一部は[scipy documentation](http://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.rv_continuous.fit.html#scipy.stats.rv_continuous.fit)にあります。 :)。私はフィットのためにオプティマイザを変更しようとしましたが、それは役に立たないようです。 –

関連する問題