R-alpha: New Incomplete Beta Function

[Kurt Hornik]

>>>>> Ross Ihaka writes:

> Here is a drop-in replacement for the R incomplete beta function.
> 	src/math/pbeta.c

> It is a slightly modified version of the cephes library one from
> Netlib.  In the few cases I tried it seems to give at least 14
> digit agreement with the one in S-PLUS (its hard to get more).

> I'm not sure what performance is like.  I'd like to know if it
> helps with some of the problems which have been reported for the
> t and F code in R.

Indeed, `qt(0.975, 3)' now also works on Debian Linux/GNU/ix86 when
compiled with `-O2 -g'.

Thanks!

(Of course, this is the only example I tried thus far ...)

-k
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-