Python 和 Scipy : How to fit a von mises distribution?

Question

我正在尝试拟合 von Mises 分布，来自 scipy ( http://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.vonmises.html )

所以我试过了

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)