2

我最近开始探索 Python 并遇到了一个名为PaCal的包的问题

一切看起来都很好，除了我在任何时候想打印一些数据时都会出现这个错误（比如在 中print A.mean()），错误行是：

Traceback (most recent call last):
  File "C:\Users\rmobenta\Desktop\tt.py", line 12, in <module>
    print A.interval(0.95)
  File "C:\Python27\lib\site-packages\pacal\distr.py", line 229, in interval
    return self.quantile(p_lim), self.quantile(1.0 - p_lim)
  File "C:\Python27\lib\site-packages\pacal\distr.py", line 215, in quantile
    return self.get_piecewise_cdf().inverse(y)
  File "C:\Python27\lib\site-packages\pacal\segments.py", line 1721, in inverse
    x = findinv(segi.f,  a = segi.a, b = segi.b, c = y, rtol = params.segments.cumint.reltol, maxiter = params.segments.cumint.maxiter) # TODO PInd, MInf
  File "C:\Python27\lib\site-packages\pacal\utils.py", line 384, in findinv
    return brentq(lambda x : fun(x) - c, a, b, **kwargs)
  File "C:\Python27\lib\site-packages\scipy\optimize\zeros.py", line 414, in brentq
    raise ValueError("rtol too small (%g < %g)" % (rtol, _rtol))
ValueError: rtol too small (1e-16 < 4.44089e-16)

我正在使用我为演示获得的两行脚本（由该软件包的作者提供）并且不知道如何解决这个问题。

这是脚本：

from pacal import *
Y = UniformDistr(1, 2)
X = UniformDistr(3, 4)
A = atan(Y / X)
A.plot()

print A.mean()
print A.interval(0.95)
4

1 回答 1

0

问题来自于 params.py : 中定义l.141的PaCalsegments.vumint.reltol = 1e-16。

这是rtol在segments.py中传递给SciPy 函数的值brentq()。

最后将其与numpy.finfo(float).eps * 2（l.413和scipy/optimize/zeros.pyl.10的）进行比较，不幸的是较小。

所以这可能是 PaCal 实现的问题，与您的代码无关。

请注意，您提供的值interval()对应于默认值（l.222distr.py ）。

我认为您应该联系 PaCal 开发人员以获取更多信息并可能打开一个问题。

于 2014-06-23T14:18:46.923 回答