我正在尝试拟合来自 scipy 的冯·米塞斯分布(http://docs.scipy.org/doc/scipy/reference/ generated/scipy.stats.vonmises.html 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,或平均值。但不知道第一个返回值是什么。
我应该如何拟合冯·米塞斯分布?
返回值是 kappa、loc 和scale。不幸的是,冯·米塞斯的 pdf 似乎并不适合拟合。如果你修正了比例,它确实适合:
>>> vonmises.fit(r, fscale=1)
(2.994517240859579, -0.0080482378119089287, 1)
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)