[SciPy-dev] scipy.stats: sf for symmetric distributions.
Wed Sep 16 16:01:47 CDT 2009
On 16-Sep-09, at 3:53 PM, David Warde-Farley wrote:
> incidentally, this breaks two
> ttest tests because they differ in the 17th decimal place, how sure
> are you of those hard-coded numbers?).
I just did the integral with Sage and requested 100 bits of precision,
more than we need, and got
1/4*sqrt(2)/(1/2*x^2 + 1)^(3/2)
sage: n(2 * integral(thepdf, x, -infinity, -1), prec=100)
And in IPython:
In : np.float64(0.42264973081037423549085121950)
So it seems that ...21 is the right answer and the hardcoded ...27 is
actually the result of roundoff, it would seem.
I've fixed the tests.
More information about the Scipy-dev