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)