[SciPy-dev] Warning about remaining issues in stats.distributions ?
Yaroslav Halchenko
lists at onerussian.com
Tue Dec 9 13:32:14 EST 2008
> * distributions that have problems for some range of parameters
so a good (imho) piece to add to unittests for the 'issues' to be fixed:
scipy.stats.rdist(1.32, 0, 1).cdf(-1.0+numpy.finfo(float).eps)
(tried on the SVN trunk to verify that it fails... discover the reason on
your own ;-))
For myself I resolved it with
__eps = N.sqrt(N.finfo(float).eps)
rdist = rdist_gen(a=-1.0+__eps, b=1.0-__eps, ....
but I am not sure if that is the cleanest way... and may be some other
distributions would need such tweakery to make them more stable.
--
.-.
=------------------------------ /v\ ----------------------------=
Keep in touch // \\ (yoh@|www.)onerussian.com
Yaroslav Halchenko /( )\ ICQ#: 60653192
Linux User ^^-^^ [175555]
More information about the SciPy-Dev
mailing list