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)

デモ用に入手した 2 行のスクリプト (このパッケージの作成者から提供されたもの) を使用していますが、この問題に対処する方法がわかりません。

スクリプトは次のとおりです。

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