tlm21@cam.ac.uk writes:
> Full_Name: Tim Massingham
> Version: 1.2.2
> OS: Debian/Linux
> Submission from: (NULL) (131.111.8.68)
>
>
> In 1.2.2 sources (also in 0.90.1. I haven't been able to check other
versions)
> Line 103 in qbeta.c should read:
> w = y * sqrt(h + r) / h - (t - s) * (r + 5. / 6. - 2 / (3 * h));
>
> since otherwise the 5 / 6 will evaluate to zero (I think).
> The Statlib fortran code uses five / six instead.
>
> Cheers, thanks for great software.
Whoops! Fortunately it sits in code that calculates an initial
approximation, so one would hope it only affects speed of convergence.
Thanks for pointing that out.
--
O__ ---- Peter Dalgaard Blegdamsvej 3
c/ /'_ --- Dept. of Biostatistics 2200 Cph. N
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk) FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._