如何解决拟合自定义分布 scipy.stats 会导致溢出
我正在尝试将广义误差分布拟合到我拥有的一些数据中。分布的形式为
我尝试了以下实现
void l(int a[45])
然而这给了
import numpy as np
import scipy.stats as st
from scipy.special import gamma
class ged(st.rv_continuous):
def _pdf(self,x,mu,sigma,kappa):
term1 = gamma(3*kappa)/gamma(kappa)
exponent = (term1 * ((x - mu)/sigma)**2)**(1/(2*kappa))
term2 = np.exp(-exponent)
term3 = 2*sigma*gamma(kappa+1)
fx = term1**0.5 * term2/term3
return fx
ged_inst = ged(name='ged')
data = np.random.normal(size=1000)
ged_inst.fit(data,0.01,1)
我如何正确实施此分发?我试图适应真实数据(不是问题中生成的玩具普通数据)
解决方法
正如评论中发布的那样,要使其正常工作,我需要覆盖默认的 _argcheck
函数。以下工作:
class ged(st.rv_continuous):
def _pdf(self,x,mu,sigma,kappa):
term1 = gamma(3*kappa)/gamma(kappa)
exponent = (term1 * ((x - mu)/sigma)**2)**(1/(2*kappa))
term2 = np.exp(-exponent)
term3 = 2*sigma*gamma(kappa+1)
fx = term1**0.5 * term2/term3
return fx
def _argcheck(self,kappa):
s = sigma > 0
k = kappa < 1
return s and k
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。