Python＆Scipy：如何适应von mises发行版？

Question

我正试图从scipy（ http://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.vonmises.html ）中修改von Mises发行版。

所以我试过了

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的参数。

返回值是多少？ 我觉得第二个是loc，或者是平均值。 但不知道第一个返回值是什么。

我应该如何适应von mises发行版？

Answer 1

返回的值是kappa，loc和scale。 不幸的是，von Mises pdf似乎并没有让自己屈服。 如果您修复比例，它确实适合：

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