Python＆amp; 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)